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SECTION  1 
INTRODUCTION 


This  section  presents  Information  that  Is  Introductory  to  the 
following  sections.  However,  It  goes  beyond  a  simple  introduction  by 
incorporating  a  summary  of  this  report  together  with  conclusions  drawn 
from  the  work  performed,  and  recommendations  regarding  future  related 
work. 


1 . 1  BACKGROUND 

The  purpose  of  developing  a  model  of  urban  mass  fires  is  to  pre¬ 
dict  the  spread,  development,  and  ultimate  effects  of  uncontrolled  urban 
fires.  The  model  predictions  need  to  be  in  sufficient  detail  and  accuracy 
to  provide  decision  makers  with  information  regarding  the  mitigation  and/ 
or  encouragement  of  selected  mass  fire  effects.  As  always,  the  ideal  is 
precise  and  accurate  predictions.  However,  the  detail  and  precision  that 
can  be  useful  to  decision  makers  may  be  far  less  than  that  required  by 
fire  phenomenology  researchers. 

For  these  reasons,  and  in  the  interest  of  providing  decision 
makers  with  a  (minimally  useful)  model  as  early  as  possible,  we  have 
defined  a  mass  fire  model  concept  and  developed  a  demonstration  implemen¬ 
tation  of  the  concept.  It  is  anticipated  that  the  phenomenology  models 
employed  may  be  criticized  as  being  premature,  imprecise,  and/or  insuffi¬ 
ciently  substantiated  by  empirical  data.  We  agree,  beforehand,  with  these 
allegations.  The  evolutionary  concept  of  simulation  development  demands 
that  the  physical  models  be  continuously  subject  to  improvement. 
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This  work  Is  concerned  with  simulating  fires  Initiated  In  urban 
areas  by  nuclear  detonations.  Such  fires  are  characterized  by  a  large 
number  of  more  or  less  simultaneous  Ignitions  distributed  over  an  area 
perhaps  hundreds  of  square  kilometers  In  extent.  The  only  fire  behavior 
experience,  which  approaches  these  circumstances.  Is  limited  to  the  mass 
fires  of  World  War  II.  These  mass  fires  resulted  from  conventional 
Incendiary  bombing  of  German  and  Japanese  cities  and  from  the  nuclear 
detonations  over  Hiroshima  and  Nagasaki.  In  addition,  some  large  wildland 
fires,  burning  in  heavy  fuels  and  rapidly  propagating  by  branding  far 
ahead  of  the  main  fire  front,  may  have  produced  sufficient  densities  of 
discrete  Ignitions  to  have  approximated  mass  fire  conditions. 

Observer  reports  and  ex  post  facto  reconstructions  of  these 
events  establish  some  qualitative  and  semiqualltatlve  features  of  mass 
fires,  which  a  valid  model  must  replicate,  but  which  alone  cannot  furnish 
a  sufficient  basis  for  developing  a  model  of  mass  fire  behavior.  A  model 
of  urban  mass  fire  behavior,  thus,  cannot  be  totally  based  upon  the 
results  of  full  scale  empirical  observations.  The  data  are  either  Inade¬ 
quate  or  nonexistent,  and  are  likely  to  remain  so.  Enough  documented  fire 
behavior  experience  does  exist,  however,  to  provide  critical  checks  upon 
the  results  of  any  mass  fire  behavior  model.  For  example,  a  useful  fire 
behavior  model  should  replicate  the  development  of  the  Hamburg  fire  storm, 
as  well  as  the  failure  of  fire  storms  to  develop  In  Berlin. 

In  the  development  of  a  computerized  simulation  of  physical 
phenomena,  it  Is  Important  to  address  the  computational  basis  of  the  simu¬ 
lation  as  well  as  its  physical  basis.  Complete  knowledge  of  the  relevant 
physical  processes  will  not  assure  development  of  a  useful  simulation, 
unless  suitably  timely,  detailed,  accurate,  and  economical  computations 
are  possible. 


The  work  described  In  this  report  addresses  both  of  these 
aspects  of  the  problem  by  describing  a  concept  for  a  computerized 
simulation  of  mass  fire  behavior,  which  Incorporates  physically  based  fire 
phenomenology  models. 

1.2  CONCLUSIONS  AND  RECOMMENDATIONS 

The  most  Important  results  of  this  work  are  the  model  concept 
and  the  demonstration  model.  The  demonstration  model  is  an  implementation 
of  the  model  concept  In  simplified  form.  The  demonstration  model  shows 
the  feasibility  of  developing  a  computerized  simulation  of  mass  fires 
based  on  the  concept  described  in  Section  4.  Additionally,  the  demonstra¬ 
tion  model  has  potential  for  both  short-term  and  long-term  utility. 

In  the  short  term,  It  is  recommended  that  the  physical  realism 
of  the  fire  phenomenology  models  be  improved,  that  the  simulation  be 
calibrated  and  that  parametric  results  be  generated.  In  the  construction 
of  the  demonstration  model,  it  was  not  possible  to  address  the  fire 
phenomenology  models  in  the  required  detail.  Consequently,  a  great  deal 
of  Improvement  in  overall  model  realism  can  be  accomplished  with  a  modest 
Investment  to  Improve  the  phenomenology  models. 

1.2.1  Hydrodynamic  Modeling 

The  total  atmospheric  perturbation,  and  the  part  of  it  describ¬ 
ing  street  level  winds,  is  essential  not  only  to  describe  mass  fire 
conditions  but  also  to  model  the  spread  and  development  of  an  urban  mass 
fire.  The  coarse  resolution  global  atmospheric  hydrodynamics  of  a  large 
area  fire  Is  a  useful,  but  limited,  approach.  Both  experimental  and  theo¬ 
retical  evidence  Indicates  that  the  actual  distribution  of  separated 
fires,  and  the  spontaneous  breakup  of  area  fires  into  separate  columns 
exercise  major  Influence  upon  fire-induced  winds.  Not  necessarily  to  the 
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exclusion  of  other  approaches,  our  model  concept  opts  for  a  hydrocalcu¬ 
lation  based  upon  parameterized  multiple  plume  properties,  which  can 
Incorporate  the  requisite  resolution.  The  feasibility,  In  principle,  of 
this  approach  is  shown  by  the  demonstration  model,  but  further  investi¬ 
gation  and  theoretical  study  is  required  to  develop  these  ideas  Into 
reliable  computational  procedures. 

1.2.2  Branding  and  Contagion 

Limits  of  spatial  resolution,  uncertainties  of  fuel  properties, 
and  insufficient  understanding  of  phenomenology  are  difficulties  that 
models  of  spread  mechanisms  In  urban  fires  must  confront,  even  when  hydro- 
dynamic  flow  fields  are  specified.  Our  demonstration  model  Includes  the 
two  propagation  mechanisms  of  branding  and  contagion  (flaming  contact). 
Radiative  Ignition  Is  omitted,  but  could  be  Included  as  well.  Some 
aspects  of  these  spread  mechanisms  are  amenable  to  simple  physical  estima¬ 
tion.  An  Important  goal  is  to  improve  the  estimates  and  extend  the  areas 
In  which  estimates  can  be  made  beyond  the  level  presently  employed  in  the 
demonstration  model.  Where  possible,  physical  estimates  Inspire  confi¬ 
dence  not  otherwise  obtainable.  In  view  of  the  difficulties  mentioned 
above,  there  are  definite  limits  upon  these  physical  refinements  and  the 
degree  to  which  they  could  be  used  to  model  propagation  mechanisms. 
Recognition  of  these  limits  and  their  implications  for  modeling  procedures 
is  an  important  research  effort  in  model  development. 

1.2.3  Simulation  Grid  and  Subgrid  Characteristics 

From  consideration  of  the  tradeoffs  between  the  desire  for 
detail  and  accuracy,  and  the  need  for  computational  economy  and  storage 
efficiency,  it  was  decided  that  the  demonstration  model  would  employ  a 
spatial  grid  of  city  block  sized  elements  (i.e.,  cells).  Over  these 
cells,  uniform  values  of  variables  and  effects  would  be  assumed.  For 


example,  a  cell  would  be  characterized  by  a  fuel  type,  the  wind  velocity 
would  be  constant  over  a  cell,  and  combustion  would  be  (spatially)  uniform 
over  a  cell.  City  block  sized  grid  cells  are  physically  appropriate  in 
that  each  cell  is,  thus,  surrounded  by  a  line  fuel  break  (i.e.,  street)  of 
specifiable  width.  In  this  grid  configuration,  the  fire,  in  order  to 
spread,  must  be  able  to  propagate  over  the  streets.  In  the  process  of 
developing  and  debugging  the  demonstration  simulation,  it  became  clear 
that  the  city  block  sized  resolution  was  Insufficient  to  describe  effects 
to  a  city  block  sized  level  of  detail.  For  example,  ignition  of  a  non- 
burning  cell  by  a  burning  cell  depends  on  how  and  where  the  fuel  is 
burning  within  the  Ignited  cell.  Indeed,  it  was  necessary  to  incorporate 
some  subgrid  detail  In  the  demonstration  model.  The  grid  structure 
employed  and  the  nature  and  detail  of  the  subgrid  calculations  can  signif¬ 
icantly  impact  the  computational  econoity  and  storage  efficiency  of  any 
mass  fire  model.  For  these  reasons,  it  is  recommended  that  a  more 
detailed  study  be  made  of  the  simulation  grid  structure  and  the  nature  and 
detail  of  the  subgrid  calculations  to  be  employed.  This  study  should 
recognize  that  these  two  subjects  are  not  independent  and  that  they  both 
relate  to  the  problem  of  fuel  specification. 

The  simulation  can  be  calibrated  against  historical  data  on 
large  fires  and  should  be  able  to  replicate  at  least  gross  details  of 
known  situations.  A  model,  calibrated  to  even  this  degree,  is  of  poten¬ 
tial  use  to  decision  makers.  Parametric  studies  of  historical  fires  in 
which  selected  parameters  are  varied  to  produce  fire  effects  as  a  function 
of  the  varied  parameter  will  provide  additional,  sensitivity  information 
for  use  by  decision  makers. 

In  the  long  term,  the  structure  of  the  demonstration  model  will 
provide  the  vehicle  for  an  evolving  mass  fire  simulation.  An  Intimate 
understanding  of  the  need  for,  and  the  nature  of  evolving  simulations  is 
explicit  in  the  design  of  the  demonstration  model  structure.  It  readily 


allows  for  modification  of  existing,  or  Incorporation  of  new,  computation 
processes  and  data  types. 


1 . 3  SUMMARY 

Section  2  presents  a  succinct  description  of  the  state-of-the- 
art  in  urban  mass  fire  modeling.  An  exhaustive  review  of  the  literature 
was  not  attempted;  however,  the  important  sources  are  discussed  and  refer¬ 
enced.  It  was  found  that,  although  work  is  in  progress  on  the  constituent 
parts  of  a  mass  fire  behavior  model,  no  satisfactory  model  exists. 

Section  3  describes  the  USOA  Forest  Service-developed  Rothermel 
model  that  has  been  used  to  predict  the  behavior  of  fires  in  lower  level 
(grass,  brush,  etc.)  wildland  fuels.  The  Rothermel  model  is  described  in 
considerable  detail  in  order  to  clearly  point  out  why  it  is  considered 
inappropriate  for  describing  fires  in  an  urban  environment. 

Section  4  presents  the  physical  basis  for  a  computerized  simu¬ 
lation  of  the  initialization  and  spread  of  uncontrolled  urban  fires, 
including  the  onset  and  development  of  mass  fires.  Such  a  simulation  must 
Incorporate  physically  based  models  of  all  the  pertinent  phenomena.  These 
models  must  parametrically  define  cause-effect  relationships  for  all  of 
the  Important  variables  and  variable  interactions  involved  in  the  effects 
to  be  simulated.  The  required  physically-based  phenomenology  models  are 
described  and  discussed  in  detail.  Particular  emphasis  is  placed  on  the 
development  of  useful  (as  opposed  to  ideal)  phenomenological  models  for 
timely  incorporation  in  a  mass  fire  simulation. 

Section  5  presents  the  computational  basis  for  a  computerized 
simulation  of  urban  mass  fires.  A  demonstration  version  of  the  simulation 
is  described  in  terms  of  the  overall  simulation  structure,  the  mathemati¬ 
cal  and  logic  structure  of  each  major  model,  and  data  entry,  management, 
and  output. 
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Section  6  presents  sample  results  for  eleven  example  cases. 
These  cases  Involve  different  initial  ignition  densities  and  patterns, 
together  with  different  Initial  ambient  wind  conditions.  These  results 
show  the  initiation  and  development  of  uncontrolled  urban  fires  In  terms 
of  fire  geometry,  heat  production  rates,  and  wind  velocities. 


SECTION  2 

MASS  FIRE  MODELING  (STATE-OF-THE-ART) 


2.1  SUMMARY 

There  presently  exists  no  satisfactory  model  by  which  the 
development  of  nuclear- induced  urban  fires  may  be  described.  Research 
into  the  constituent  mechanisms  required  as  components  of  such  a  model  is 
in  progress,  but  model  development  need  not  wait  for  this  research  to  be 
completed.  Modeling  is  possible  at  any  state  of  knowledge,  and,  with 
proper  structuring,  can  be  done  in  an  open-ended  fashion  to  Incorporate 
new  knowledge  when  it  develops. 

Significant  advances  of  understanding  have  been  acquired  since 
the  existing  models  were  devised.  The  need  to  incorporate  fire  spread  and 
mass  fire  phenomena  as  related  parts  of  a  unified  simulation  Is  the  major 
requirement,  which  can  be  met  at  this  time.  The  potentially  lengthy  and 
costly  effort  to  determine  fuel  properties  and  inventories  can  be  guided 
by  the  degree  to  wh  ich  any  practical  model  can  make  use  of  them. 

2.2  BACKGROUND 

The  initiation,  development,  and  effects  of  mass  fires  resulting 
from  nuclear  attack  on  urban  areas  has  been  a  problem  of  concern  since  the 
advent  of  nuclear  weapons.  Investigations  motivated  by  this  problem  have 
amassed  a  volume  of  material,  dealing  with  its  various  aspects,  which  is 
summarized  in  the  review  of  Martin  (Reference  1),  supplemented  by  the 
recent  report  of  Small  and  Brode  (Reference  2).  Both  of  these  documents 
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contain  extensive  bibliographies.  Since  1978,  the  Defense  Civil  Prepared¬ 
ness  Agency,  now  part  of  the  Federal  Emergency  Management  Agency,  has 
sponsored  an  Integrated  research  program  aimed  at  achieving  a  reliable 
method  for  predicting  the  Incendiary  outcome  of  nuclear  attack.  A 
periodic  report  and  assessment  of  this  program  has  been  published  yearly 
(References  3  and  4).  The  following  discussion  here  does  not  attempt  to 
repeat  these  general  surveys,  but  focuses  specifically  upon  certain  model¬ 
ing  problems. 

2.3  INITIAL  CONDITIONS 

Much  of  the  previous  and  presently  ongoing  research  Is  devoted 
to  establishing  the  initial  fire  conditions  resulting  from  nuclear 
attack.  The  effect  of  the  thermal  pulse,  as  modified  by  distance  and 
atmospheric  conditions,  in  producing  ignition  of  various  inflammable 
materials  exposed  to  it;  the  potential  of  these  ignitions  to  spread  and 
result  in  full  building  involvement;  and  the  effect  of  blast  to  (I)  alter 
substantially  the  fuel  properties  and  distribution;  (2)  to  extinguish  or 
enhance  primary  Ignitions;  and  (3)  to  produce  secondary  Ignitions  by 
breakage  of  gas  and  power  lines  are  all  important  aspects  of  the  initia¬ 
tion  problem  (Reference  4).  These  very  important  processes,  which  others 
are  investigating,  are  not  our  direct  concern  here.  Given  the  initial 
conditions  (of  both  ignitions  and  fuel),  we  are  Interested  in  modeling  the 
fire  development. 

Though  not  our  immediate  concern,  the  faithful  incorporation  of 
these  very  detailed  mechanisms  to  initialize  a  model  is  an  outstanding 
problem.  Given  uncertainties  of  structural  behavior  and  the  impossibility 
of  resolving  small  detail  in  a  city-wide  simulation,  some  judicious  model¬ 
ing  choices  are  required,  even  assuming  adequate  understanding  of  basic 
mechanisms. 


21 


2.4 


SPREAD  MODELS 


Existing  fire  spread  models  are  reviewed  by  Weisbecker  and  Lee 
(Reference  5)  and  Alger  and  Martin  (Reference  3).  Many  models  of  fire 
spread  within  compartments  and  buildings  have  been  built,  but  attempts  at 
spread  models  on  a  cit/-wide  scale  are  limited  to  three  models.  These  are 
the  SSI  (Scientific  Services,  Inc.),  the  URS  (URS  Research  Company),  and 
the  IITRI  (Illinois  Institute  of  Technology  Research  Institute)  models  and 
their  modifications.  Relying  upon  the  conclusions  of  the  papers  cited 
above,  It  appears  that  all  of  these  models  (which  were  developed  in  the 
late  1960's)  may  suffer  serious  defects. 

All  of  these  models  stress  fuel  distributions,  which  are  incor¬ 
porated  at  a  level  of  detail  ranging  from  high  detail  (actual  building  and 
content  surveys)  in  the  SSI  model  to  use-factor  classification  of  areas  in 
the  URS  and  IITRI  models.  Fuel  representation  does  not  account  for  blast 
damage,  which  is  viewed  as  a  serious  restriction.  The  basic  spread  mech¬ 
anism  of  all  these  models  is  radiative  contagion,  with  some  simplified 
spread  by  firebrands  in  the  IITRI  and  URS  models.  The  failure  to  ade¬ 
quately  model  fire  parameters,  upon  which  these  processes  depend,  is 
considered  a  shortcoming.  In  practice  these  models  spread  the  fire  in  a 
stochastic  manner,  based  upon  fuel  type.  A  final,  and  perhaps  decisive, 
criticism  of  these  pioneering  models  is  their  lack  of  any  interactive 
features  of  mass  fire  phenomena,  fire- induced  winds,  their  effect  upon 
spread,  and  plume  coalescence. 

Consistent  with  the  Project  Flambeau  results,  we  do  not  agree 
with  the  emphasis  upon  radiation  as  a  primary  spread  mechanism,  but 
believe  that  both  branding  and  flame  contact  are  of  equal,  if  not  greater, 
importance.  Mass  fire  phenomena  are  essential  to  determine  spread.  With¬ 
out  minimizing  the  true  role  of  fuel  properties,  it  may  be  that  this 
problem  aspect  has  been  overemphasized  in  the  past.  Resolution  limita¬ 
tions  will,  in  any  model,  limit  the  detail  of  fuel  descriptions. 
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In  any  case,  the  accumulation  of  understanding  permits  substan¬ 
tial  Improvement  upon  these  early  spread  models, 

2.5  FIRESTORM  MODELS 

The  presence  of  strong,  hurricane-force  winds  directed  radially 
inward  toward  the  large  urban  fires  of  WWII  has  prompted  diverse  attempts 
to  understand  and  model  this  phenomenon,  called  a  "firestorm." 

An  analysis  of  historical  experience  (Reference  6)  established 
some  preliminary  conditions  thought  necessary  for  firestorm  phenomena.  By 
employing  historical  analysis  in  conjunction  with  modeling,  Lommanssen,  et 
al  (Reference  7)  developed  buildup  and  existence  criteria  for  firestorms. 
Firestorm  criteria  were  also  developed  (Reference  8)  on  the  basis  of  WWII 
experience  and  empirical  data  on  flame  merging. 

The  need  for  experimental  data  motivated  two  programs  to  test 
the  behavior  of  large  fires:  the  Flambeau  program  (References  9  and  10), 
and  the  Euroka  fire  in  Australia  (Reference  11),  a  fire  of  260  meter 
radius.  The  data  from  these  fires,  though  none  were  as  large  as  an  urban 
mass  fire,  fills  an  important  gap  between  historical  mass  fires  and  labo¬ 
ratory  fires.  We  have  found  the  insights  reported  by  Countryman  from  the 
Flambeau  program  especially  valuable. 

Aside  from  fire  plume  modeling,  done  by  several  Investigators, 
the  theoretical  work  of  Nielsen  (Reference  12)  and  Nielsen  and  Tao  (Refer¬ 
ence  13)  specifically  models  the  global  hydrodynamics  of  a  large  urban 
fire.  Recently,  a  new  approach  to  atmospheric  disturbances  and  physical 
parameters  of  large  area  fires  has  been  initiated  (Reference  2).  Some 
detailed,  two-dimensional  analytic  models  of  fire  plumes  In  the  atmosphere 
have  been  published  (References  14  through  17). 


We  believe  the  hydrodynamic  properties  of  urban  mass  fires  to  be 
of  central  importance.  The  extent  to  which  global  (lew-resolution) 
hydrodynamics  Is  sufficient  needs  further  examination.  The  firestorm 
phenomenon,  however,  Is  only  one  aspect  of  the  problem.  The  hydrodynamic 
modeling  must  be  broadly  enough  structured  to  encompass  any  developmental 
behavior  which  occurs.  A  good  part  of  the  problem  Is  to  understand  what 
can  happen. 

2.6  WILDLAND  FIRES 

Our  task  excludes  wildland  fires,  but  any  discussion  of  the 
state-of-the-art  for  modeling  large  fires  must  take  cognizance  of  the 
Rothermel  Model  (Reference  18),  which  occupies  a  unique  status  as  the  only 
well-verified  spread  model  in  existence.  This  model,  its  limitations,  and 
possible  applicability  to  the  urban  fire  problem,  are  examined  In  detail 
in  Section  3. 


SECTION  3 

ANALYSIS  OF  ROTHERMEL  RATE  OF  SPREAD  EQUATION 


3.1  SUMMARY  AND  CONCLUSIONS 

The  Rothermel  rate  of  spread  (RROS)  equation  (Reference  18)  for 
estimating  localized  rate  of  fire  spread  in  a  homogeneous  fuel  bed  is 
described  in  terms  of  a  qualitative  theoretical  basis  and  a  quantitative 
empirical  definition.  Generalizations  of  the  basic  rate  of  spread  equa- 
tion  are  described  that  allow  estimation  of  fire  spread  rates  over  an  area 
of  interest  in  heterogeneous,  natural  understory  fuels.  It  is  shown  that 
local  scalar  spread  rate  is  necessary  but  insufficient  to  allow  fire 
growth  estimates  to  be  made.  A  fire  front  propagation  algorithm  and  a 
local  directional  spread  rate  algorithm  are  described,  together  with  local 
scalar  spread  rate,  that  provide  a  computational  procedure  for  estimating 
the  growth  of  wildland  fires.  Difficulties  associated  with  input  data 
acquisition  and  data  management  resulting  from  the  fuel  characteristics 
required  by  the  rate  of  spread  equation  are  discussed  and  their  resolution 
is  described. 

The  Rothermel  rate  of  spread  equation  is  demonstrably  suitable 
(Reference  19)  for  estimating  the  growth  of  free  burning  fires  in  lower 
story,  natural  fuels.  From  the  analyses  described  here,  and  in  the 
referenced  report.  It  is  concluded  that  the  Rothermel  equation  Is  of 
limited  utility  in  estimating  the  growth  of  urban  fires  for  the  following 
reasons: 


1.  The  quantitative  coefficients  of  the  RROS 
equation  are  valid  only  for  understory 
natural  fuels  ( 1 .e. ,  grass,  brush,  etc.). 
They  were  experimentally  determined  from 


laboratory  fires  In  continuous  fuel  beds 
where  the  fuel  was  simulated  by  uniform, 
homogeneous  constructs  of  uniform  fuel 
particles  and  the  environment  was  simu¬ 
lated  by  a  wind  tunnel. 

2.  The  fire  front  propagation  and  local 
directional  spread  rate  algorithms, 
required  for  wildland  fire  growth  esti¬ 
mates  were  developed  specifically  for 
understory  natural  fuels  and  are  not 
applicable  to  estimating  the  growth  of 
fires  In  an  urban  fire  environment. 

3.  It  Is  explicit  in  the  quantitative 
coefficients  of  the  RROS  equation  that 
wildland  fires  in  natural  understory 
fuels  are  largely  propagated  by  the 
combustion  of  loosely  packed  fine  fuels 
at  the  periphery  of  the  fire;  loosely 
packed  fine  fuels  do  not  occur  in  signif¬ 
icant  quantities  in  urban  fuel  beds. 

4.  The  RROS  equation  describes  a  free  burn¬ 
ing  fire  In  which  venting  effects  are 
quantitatively  controlled  by  fuel  pack¬ 
ing.  In  an  urban  fire  venting  effects 
are  generally  controlled  by  the  building 
structure  within  which  the  fire  is  burn¬ 
ing. 

5.  Wind  effects  in  the  RROS  equation  are 

open  loop;  It  has  wind-fire  coupling,  but 
no  fire-wind  coupling. 

6.  The  effects  of  branding  are  not  addressed 
by  the  RROS  equation. 

7.  It  is  not  economically  nor  physically 

practicable  to  experimentally  quantify 

the  detailed  Interparticle  heat  exchange 
mechanisms  for  all  urban  fuels  of  Inter- 
est  In  analogy  with  the  experiments  that 
were  performed  for  natural  fuels. 

8.  There  is  no  mechanism  In  the  RROS  for 

allowing  fire  propagation  across  a  break 
in  the  fuel  bed,  analogous  to  a  street  In 
the  urban  environment. 
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The  RROS  can  be  used  to  describe  fire  growth  In  significant 
areas  of  natural  fuels  within  the  urban  environments  cf  interest.  Even 
for  this  purpose,  however,  modifications  will  be  necessary  to  account  for 
trees  and  the  energy  release  rate  of  unconsumed  fuels  burning  within  the 
fire  perimeter.  These  conclusions  regarding  the  inapplicability  of  the 
RROS  equation  to  the  growth  of  urban  fires  contribute  nothing  positive 
toward  the  solution  of  this  problem.  However,  Section  4  of  this  report 
presents  a  methodology  for  describing  the  growth  of  fires  in  urban  envi¬ 
ronments  of  Interest. 

3.2  THEORETICAL  BASIS 

The  theoretical  basis  of  the  RROS  is  largely  due  to  Fons  (Refer¬ 
ence  20),  Tarifa  and  Torralbo  (Reference  21),  and  Frandsen  (Reference 
22).  Fons  was  the  first  to  describe  fire  spread  by  a  mathematical  model. 
The  Fons  model  embodied  a  contagion  concept  in  which  fire  spread  in  a  fuel 
bed  is  considered  to  consist  of  successive  ignitions  occurring  at  a  rate 
controlled  primarily  by  the  ignition  times  and  distances  between  succes¬ 
sive  fuel  particles.  Tarifa  and  Torralbo  confirmed  Fons'  ideas  by  showing 
that  heating  of  the  fuel  ahead  of  a  progressing  flame  front  is  the  first 
and  most  essential  process  of  the  flame  propagation  mechanism.  Frandsen 
applied  the  principle  of  energy  conservation  to  a  unit  volume  of  fuel 
ahead  of  an  advancing  flame  front  in  a  homogeneous  fuel  bed  to  obtain  the 
following  expression  for  the  quasi-steady  rate  of  fire  spread: 


I  .  +  /  (al  /az)  dx 

Xig  Z  'zc 

‘’be  Qig 


(3-1) 
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where 


R  s  quasi-steady  rate  of  spread  (ft/min) 

Ixig  a  horizontal  heat  flux  absorbed  by  a  unit  volume  of  fuel 
at  the  time  of  ignition  (Btu/ft2-min) 

Pbe  s  effective  bulk  density  (the  amount  of  fuel  per  unit 
volume  of  the  fuel  bed  raised  to  ignition  temperature 
ahead  of  the  advancing  fire)  (lb/ft3) 

(3lz/3z)2c  s  the  gradient  of  the  vertical  intensity  evaluated 
in  a  plane  at  a  constant  depth,  zc,  of  the  fuel  bed 
(Btu/ft-min) 

Qig  -  heat  of  ignition  of  unit  volume  of  fuel  (Btu/lb). 

Equation  3-1  is  a  theoretically  rigorous  description  of  the  rate 
of  spread  of  a  line  fire  through  a  homogeneous  fuel  bed  of  uniform  depth 
in  terms  of  a  ratio  of  the  heat  flux  received  by,  to  the  heat  required  for 
ignition  of,  a  unit  volume  of  fuel  ahead  of  the  advancing  flame  front. 
Equation  3-1  contains  heat  flux  terms  for  which  the  mechanisms  of  heat 
transfer  are  not  known,  and  consequently,  it  cannot  be  solved  analytic¬ 
ally.  Definition  of  an  approximate  form  of  Equation  3-1  that  can  be 
solved  analytically  resulted  from  examination  of  each  term  and  determina¬ 
tion  of  analytical  and  empirical  methods  for  evaluating  the  unknown 
mechanisms  of  heat  transfer  (Reference  18). 

The  form  of  Equation  3-1  was  simplified  by  the  definition  of  the 
following  set  of  terms. 
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3.2.1  Effective  Heating  Number 


The  amount  of  fuel  involved  in  the  ignition  process  Is  specified 
In  the  denominator  of  Equation  3-1  as  the  effective  bulk  density  ptj6. 

The  effective  heating  number  e  is  defined  as  the  dimensionless  ratio  of 
the  effective  bulk  density  to  the  actual  bulk  density.  Hence, 

e  *  pbe/pb  •  (3~2) 

3.2.2  Propagating  Flux 

The  numerator  of  3-1  is  defined  to  be  the  propagating  flux 

Ip  where 

■  >x1,  +  i  (3VHC  d)<  '  <3‘3> 

The  propagating  flux  consists  of  two  terms,  the  horizontal  flux  and  the 
cumulative  gradient  of  the  vertical  flux  up  to  the  fire  front.  The  verti¬ 
cal  flux  is  most  significant  for  wind  driven  and  upslope  fires,  for  which 

there  is  an  effective  flame  tilt  and  an  increased  heat  transfer  to  the 

unignited  fuel  by  radiation,  convection,  and  conduction.  It  is  assumed 
that  the  vertical  flux  is  small  for  no  slope,  no  wind  fires;  and  the  heat 
flux  for  these  conditions  Is  defined  as 

(Vo  *  ‘xig  <3-4> 

3.2.3  Reaction  Intensity 

The  energy  release  rate  of  the  fire  front  produced  by  burning 

gases  resulting  from  pyrolysis  of  the  fuel  bed  ahead  of  the  fire  front  is 

defined  as  the  reaction  intensity  Ir  where 


i 

i 


1  =•  -h(dw/dt) 

r 

and 


dw/dt 


»  mass  loss  rate  per  unit  area 
the  fire  front,  (lb/ft2-min) 


In 


(3-5) 


h  3  heat  content  of  fuel  (Btu/lb). 

The  reaction  intensity  is  the  source  of  the  no  wind  propagating  flux, 
(I  )  and  it  is  assumed  that  the  no  wind  propagating  flux  and  the 
reaction  intensity  Ir  can  be  evaluated  independently  and  then 
correlated.  Thus,  the  coefficient  relating  (Ip)0  and  *r  is 


_  (Vo  (3-6) 

—  .  ■  — —  »  • 


Knowledge  of  this  correlation  will  allow  (Ip)0  to  be  determined  from 
I  which  can  be  determined  from  parameters  that  describe  the  fuel  bed 

complex. 

3.2.4  The  Effects  of  Hind  and  Slope 

The  effects  of  wind  and  slope  are  defined  as  changes  in  the 
propagating  heat  flux,  that  alter  the  heat  received  by  potential  fuels  due 
to  convection  and  radiation.  Based  on  the  assumptions  that  the  con¬ 
tribution  of  the  vertical  intensity  (al/az)  to  (Ip)0  is  small  and  that 
the  effects  of  wind  and  slope  on  the  propagating  flux  are  independent,  the 
total  propagating  flux  is  defined  to  be  of  the  form 

ip  ■  Op)  •  o  ♦  V  •*)  ’  (3'7> 
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where 


=  a  dimensionless  coefficient  dependent  on  wind  speed, 

*s  =■  a  dimensionless  coefficient  dependent  on  slope. 

3.3  EMPIRICAL  DEFINITIONS 

Substitution  of  Equations  3-2  through  3-7  Into  Equation  3-1 
yields  the  Rothermel  rate  of  spread  (RROS)  equation  in  the  form 

(1  +  *  +  A  ) 

R  a  I  £  1- — - L  (ft/min)  .  (3-8) 

r  pb  ^ig 

Each  term  of  Equation  3-7  was  analyzed  and  transformed  into  a  format 
appropriate  for  experimental  quantification.  •  This  resulted  in  a  set  of 
experimentally  determined  relationships  by  which  rate  of  spread  estimates 
can  be  obtained  from  3-8  in  terms  of  parameters  describing  the  fuel  bed  . 
terrain  slope,  and  wind  speed  (Reference  18).  These  experimentally  deter¬ 
mined  relationships  (Reference  7)  are  given  in  Table  3-1. 

The  use  of  Equations  3-8  through  3-26  to  provide  estimation  of 
fire  spread  rate  requires  specification  of  the  parameters  shown  in  Table 
3-2. 

3.4  GENERALIZATION  OF  THE  ROTHERMEL 
RATE  OF  SPREAD  (RROS)  EQUATION 

The  RROS  equation  given  by  Equations  3-8  through  3-26  provides 
local  rate  of  spread  in  a  homogeneous  fuel  bed.  Use  of  the  RROS  equation 
in  describing  a  propagating  fire  front  necessitates  generalization  to 
include  spread  rates  for  heterogeneous  fuels  distributed  over  an  area. 
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Table  3-1.  Defining  equations  for  RROS  equation. 


Equation 

Description 

Equation 

Number 

V  ‘ r  Ws 

Reaction  intensity 
(Btu/ft  -min) 

(3-9) 

vrtiere: 

r'  '  I"max 

Optimum  reaction 
velocity  (min-1) 

(3-10) 

rmax  a  a1-5  (495  +  -0594a1. s}_1 

Maximum  reaction 
velocity  (min-1) 

(3-11) 

$op  »  3.348a-.®189 

Optimum  packing 
ratio 

(3-12) 

A  -  1/(4. 774a.1  -  7.27) 

(3-13) 

2  Mf3 

V  1  * 2-59  if +  »•“  (r 1  - 3-52  [if) 

X  X  X 

Mo i sture 

dampening 

coefficient 

(3-14) 

ns  -  0.174  e^. 19 

Mineral  damping 
coefficient 

(3-15) 

5  -  (192  +  0.2595a)-1 

X  exp [0.792  +  0.681a.5)(8  +  0.1)] 

Propagating 
flux  ratio 

(3-16) 

*  *  cuB(  ) E 

pop 

Wind 

coefficient 

(3-17) 

C  ■  7.47  exp  (-0.133a.55) 

(3-18) 

B  -  0.02526a.54 

(3-19) 

E  *  0.715  exp  (-3.59  x  10“4a) 

(3-20) 
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Table  3-1.  Defining  equations  for  RROS  equations  (concluded). 


w  .  0 

n  l  +  Sj 

Net  fuel  loading 
(lb/ft2) 

(3-21) 

a  5.275  3“  ,3(tan  $)2 

Slope  factor 

(3-22) 

pb  a  V6 

Ovendry  bulk 
density  (lb/ ft3) 

(3-23) 

e  =  exp(-138/o) 

Effective  heating 
number 

(3-24) 

Qig  »  250  +  1,116  Mf 

Heat  of  pre- 
ignition  (Btu/lb) 

(3-25) 

PP 

Packing  ratio 

(3-26) 
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Table  3-2.  RROS  equation  Input  paraaeters . 


A  generalization  of  the  RROS- equation  to  Include  heterogeneous 
fuels  Is  described  in  Reference  18  (pp.  28-33).  This  generalization 
allows  definition  of  an  arbitrary  number  of  uniformly  distributed  fuel 
categories  and  size  classes.  For  natural,  understory  fuels,  two  categor¬ 
ies  (live  and  dead)  and  four  to  five  size  classes  are  used.  Since  neither 
all  size  classes  nor  all  categories  are  equally  effective  in  fire  front 
propagation,  a  weighting  procedure  is  described  in  Reference  18  for  com¬ 
bining  the  descriptions  of  understory  fuel  size  classes  and  categf  es 
into  a  single  spread  rate. 


A  generalization  to  allow  computation  of  spread  rates  over  an 
area  is  accomplished  by  defining  a  computational  mesh  over  the  area  of 
interest  and  specifying  all  of  the  parameters  required  for  the  RROS  equa¬ 
tion  in  each  cell  of  the  mesh.  The  specification  of  mean  parameter  values 
for  the  1-th  fuel  category  and  j-th  fuel  size  class  requires  estimation  of 


each  input  parameter  shown  in  Table  3-3  for  each  cell  of  the  computational 
mesh. 

Two  weighting  parameters  based  on  the  surface  area  of  the  fuel 
within  each  category  and  size  class  are  defined  by 


and 


U'  A. 


Ratio  of  surface  area  of  the  j-th 
size  class  to  total  surface  of  the 
i-th  category  per  unit  fuel  cell 


(3-27) 


f  ■  L 


Ratio  of  surface  area  of  the 
i-th  category  to  total 
surface  area  per  unit  fuel  cell. 


(3-28) 


Table  3-3.  RRQS  equation  Input  parameters  within  the  1-th  fuel 
category  and  j-th  fuel  size  class  of  each  fuel  cell. 


Parameter 

Description 

=  ovendry  loading  (lb/ft2) 

3  surf ace-area-to-volume  ratio  (ft2/ft3) 

s  mineral  content  (lb  minerals/lb  wood) 

<Te>y 

,  .  .(lb  minerals  -  lb  silica), 

=  effective  mineral  content  ( - - - — ) 

=  low  heat  value  (Btu/lb) 

«r>« 

3  moisture  content  (lb  moisture/ lb  wood) 

^Pp^ij 

=  ovendry  particle  density  (lb/ft3) 

Mean  value  within  the  i-th  category: 

<"x>i 

3  moisture  content  of  extinction 
(lb  moisture/ovendry  wood) 

Mean  fuel 

irray  properties: 

I 

3  depth  of  fuel  (ft) 

tan  * 

=  slope  (ft  vertical  rise/ft  horizontal) 

U 

=  wind  velocity  at  midflame  height  (ft/min) 

m 

=  total  number  of  categories 

n 

3  number  of  size  classes  within  the  i-th  category. 

36 


where 


A  * 


(«) 


1j 


Mean  total  surface  area  of  fuel 
of  1-th  category  and  j-th  (3-29) 
size  class  per  unit  fuel  cell 


j=n 
Ai  =  2A1 
j=l 


J 


Mean  total  surface  area  of 
fuel  of  1-th  category  per 
unit  fuel  cell . 


(3-30) 


__  i=m 

At  ■  IA, 
1-1 


Mean  total  surface  area  of  (3-31) 
fuel  per  unit  cell. 


In  terms  of  the  above  weighting  factors  and  the  parameters  defined  in 
Table  3-3,  the  reaction  Intensity  Is  given  by 


v 


iam 

r'i?i  f^n)/i^nd)i(nM)i  * 


(3-32) 


where  the  characteristic  parameters  weighted  by  surface  area  are  as 
follows: 


j-n 

J51f«(w»)u 


Net  loading  of 
i-th  category 


(3-33) 


1  +  <Vij 


Net  loading  of  j-th  class  (3-34) 
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J-n 

h1  ■  J.fuwu 


Low  heat  content  value 
of  1-th  category 


(3-35) 


<\>1  ■  "-IMIS,),-1* 


Mineral  damping  coefficient  (3-36) 
of  i-th  category 


-  j=n 

(S  )  -  l  f  .(S ,)' 
e  1  j=1  1j  e  1j 


Characteristic  effective  • 

mineral  coefficient  of  (3-37) 
1-th  category 


■  1  -  2-S9<r„)1  ♦  5.11(fM)2  -  3.52(fH)| 


(3-38) 


Moisture  damping  co-efficient 
of  i-th  category 


Wi 


(Mf  )i 


ext-*! 


Moisture  ratio  of 
i-th  category. 


(3-39) 


j3n 


Wl  ■“  ^ij^ij 


Moisture  content  of 
i-th  category. 


The  potential  reaction  velocity  (r1)  is  an  experimentally  determined 
function  of  the  fuel  bed  packing  ratio  (e)  and  fuel  particle  size  as 
characterized  by  fuel  particle  surface-to-volume  ratio  (or).  Thus, 


r‘  3  rW0/0op)Aexp[A(1  "  e/0op)l  * 


(3-41) 


where 


F'  a  (S)1-*/(495  +  0.0594  (o)1*5)  ,  (3-42) 

maX  1 

Top  -  3.348  (o)0*8189  (3-43) 

A  ■  (4.774  (a)0,1  -  7.27)“ 1  (3-44) 


and  where 


a 


i=m 

I  V* 

i=l 


Characteristic  surface- area-to-  (3-45) 
volume  ratio  of  the  fuel  complex. 


j=n  ___ 


Characteristic  surface  area-to-  (3-46) 
volume  of  i-th  fuel  ratio  category. 


If  a 


i=m 


J-n 

l 

J-l 


/ 


Mean  packing  ratio  (3-47) 


,  i=m  j=n 

I  I  Poll! 

5  1-1  j=l 


(3-48) 


This  completes  the  calculations  for  the  reaction  Intensity  (Ir)  term  of 
the  RROS  Equation  3-8  for  nonhomogeneous  fuels.  The  remaining  terms  are 
treated  similarly.  The  no-wind  propagating  flux  ratio  is  given  by 
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f; 


£  -  (192  +  0.2595  o)_1  exp [(0.792  +  0.681  o-5)Cp:  +  0.1)]  (3-49) 

The  denominator  of  Equation  (8),  which  contains  the  heat  sink  terms,  Is 
given  by 

1**m  jan  ,,Q  _ 

PbeqigBPb  l  fi  l  f1jt«*P(a^)]  (?ifl)ij  *  (3-50) 

i-l  J-l  °1j 

where 

(Qig)ij  s  250  +  1,116  (Mf)^j  The  heat  of  pre-ignition  (3-51) 

for  j-th  size  class  within 
the  1-th  category. 

The  wind  and  slope  multiplication  factors  are  given  by 


cuB(?/?op)'E 

(3-52) 

5.275  (?)-°-3(tan  $)2 

(3-53) 

where, 

U  a  mean  windspeed  at  midflame  height  (ft/mi n) 

C  *  7.47  exp(-.133  o’55)  (3-54) 

8  *  0.02526  o-slt  (3-55) 

E  *  0.715  exp(-3.59  x  10“4  o)  (3-56) 
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3-5  PROPAGATION  MODEL  EMBODYING  THE  RROS  EQUATION 


The  generalized  RROS  equation  provides  estimates  of  local  rate 
of  fire  spread  In  heterogeneous  fuels  over  an  area  of  Interest.  This  Is, 
however.  Insufficient  to  allow  estimates  of  fire  growth,  which  require  In 
addition  a  fire  front  propagation  algorithm  and  a  local  directional  spread 
rate  algorithm. 


3.5.1  Fire  Front  Propagation  Algorithm  (FFPA) 


The  FFPA  described  here  was  developed  specifically  for  fires  In 
understory  wildland  fuels  (Reference  19)  and  is  probably  not  applicable  to 
urban  fuels.  It  Is  described,  however,  as  an  example  of  the  types  of 
considerations  Involved  in  applying  the  RROS  equation  to  fire  growth 
estimates. 


It  is  assumed  that  all  propagation  interactions  occur  at  the 
perimeter  of  a  wildland  fire  at  the  interface  between  the  flame  front  and 
unignited  fuel.  FFPA  has  these  characteristics: 

1.  The  algorithm  must  accept  required  data  and 
provide  results  of  variable,  user  specified 
resolution. 

2.  Perimeter  resolution  must  be  variable  from 
f ire-to-f ire,  perimeter-to-perimeter,  and 
point-to-point  within  a  perimeter. 

3.  A  perimeter  should  not  cross  an  earlier  perim¬ 
eter  or  itself.  That  is,  a  fire  is  assumed  to 
always  grow  toward  unburned  fuel. 

4.  The  algorithm  should  be  able  to  simulate  fires  ,  ;> 

burning  through  nonhomo geneous,  discontinuous 
wildland  fuels  in  uneven  terrain. 

5.  The  shape  of  the  fire  should  be  a  priori, 
indeterminate,  and  dependent  upon  local  fuel, 
topographical  and  meterological  conditions 
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over  the  region  within  which  the  fire  Is  situ¬ 
ated. 

6.  It  should  be  possible  to  initiate  the  model 
from  an  arbitrary  perimeter  definition  (a 
point,  geometric  shape,  set  of  points,  etc.). 

In  the  FFPA,  each  perimeter  is  described  by  a  set  of  points  con¬ 
nected  by  straight-line  segments.  The  basic  function  of  the  propagation 
model  is  to  estimate  the  location  and  shape  of  the  next  perimeter  by  an 
outward  propagation  of  the  points  that  define  segments  of  the  current 
perimeter.  The  use  of  variable  (i.e.,  segment  length)  resolution  over 
each  perimeter  permits  detailed  description  of  the  perimeter  where  desired 
while  allowing  a  significant  reduction  in  the  total  number  of  calculations 
and  the  amount  of  data  required  to  describe  the  perimeter. 

The  user  specifies  a  minimum  (5s)  and  a  maximum  (as)  allowable 
spacing  between  any  two  adjacent  points  within  a  perimeter.  An  internal 
algorithm  defines  the  spacing  between  each  successive  pair  of  points  such 
that 

1.  The  spacing  falls  within  the  specified  limits, 

2.  The  spacing  is  a  function  of  the  maximum  of  the 
propagation  rates  at  the  two  points,  and 

3.  The  greatest  density  of  points  occurs  where  the 
fire  front  propagation  rates  are  highest. 

Thus,  the  spacing  ds^j  between  two  adjacent  points  p-f  and  pj  is 
given  by 


and 


(3-58) 


AS  <_  ds^j  <_  AS  , 


where 

R.j  s  propagation  rate  at  point  p^, 

Rj  a  propagation  rate  at  point  pj. 

Let  a  perimeter  at  time  T0  be  completely  defined  by  the  set  of 
10  points  shown  in  Figure  3-1.  By  completely  defined,  we  mean  that  no 

further  operations  are  required  upon  the  specified  set  of  points  to 
complete  the  description  of  the  perimeter.  Since  the  function  of  the 

propagation  model  is  to  estimate  the  next  perimeter  in  terms  of  the  set  of 
points  that  define  the  current  perimeter,  it  does  not  matter  how  the 
perimeter  at  T0  was  obtained.  It  could  have  been  calculated  from  an 

earlier  perimeter,  it  could  have  been  estimated,  or  it  could  have  btjn 
obtained  as  input  from  measured  data.  Given  the  perimeter  at  T0,  we 

will  very  briefly  describe  the  sequence  of  operations  required  to  estimate 
the  next  perimeter. 

The  outward  direction  at  each  point  of  the  perimeter  is  defined 
as  the  bisector  of  the  angle  formed  by  the  two  line  segments  intersecting 
at  each  point,  as  shown  In  Figure  3-2.  The  directional  spread  rate  is 
calculated  at  each  point;  the  propagation  direction  is  determined  at  each 
point  (the  outward  direction  is  used  in  this  example);  and  each  point  is 
propagated  in  the  specified  direction,  at  the  specified  rate  over  a  speci¬ 
fied  interval  to  the  next  perimeter  at  T^  as  shown  in  Figure  3-3. 

The  newly  defined  set  of  points  that  describe  the  perimeter  at 
Tj  are  tested  to  determine  if  Equations  3-1  and  3-2  are  satisfied  for 
each  candidate  perimeter  segment,  and 


43 


Figure  3-1.  The  set  of  points  defining  / 

the  perimeter  at  T0.  ' 

Figure  3-2.  The  outward  direction 
at  each  point. 


Figure  3-3.  Preliminary  perimeter 
at  Tj. 


Figure  3-4.  Final  perimeter  at  Tj 


1.  a  new  point  Is  added  at  the  center  of  each 
segment  that  Is  too  long, 

2.  a  point  is  deleted  from  each  segment  that  is 
too  short,  and 

3.  the  process  Is  repeated  until  all  segments 
satisfy  the  above  together  with  Equations  3-1 
and  3-2. 

For  this  example  a  point  was  added  to  each  of  segments  (5,6), 
(7,8),  and  (8,9),  as  shown  in  Figure  3-4. 

The  collection  of  points  is  then  tested  to  determine  If  the 
perimeter  contains  any  loops  (i.e.,  crosses  itself).  If  such  a  condition 
is  detected,  the  points  contained  within  each  minor  loop  are  deleted.  The 
remaining  points  are  enumerated,  and  the  current  perimeter  is  by  defini¬ 
tion  completely  described  by  this  set  of  points,  as  shown  in  Figure  3-4, 

and  the  entire  sequence  of  operations  is  repeated  to  obtain  the  next 

perimeter. 

3.5.2  Local  Directional  Spread  Rate  Algorithm  (LDSRA) 

Local  spread  rate  estimates  are  based  on  the  RROS  Equation  3-8, 

which  is  rewritten  for  convenience  in  the  form  (Reference  5) 

R(q1,q2,./.  .qN,u,$)  =  RQ(q1 ) [1.0  +  ^(u.q^  +  ♦jU.fO]  ,  (3-59) 

where 

R(q1,q9,...qN,u,$)  a  spread  rate  in  direction  of 
1  L  n  local  wind  velocity, 

qN  *  fuel  parameters, 

u  ■  local  wind  speed 

^  *  angle  of  local  slope  measured  from 
the  local  horizontal, 
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R0(qi)  ■  no  wind,  no  slope  spread  rate, 


$w(u,q^)  *  wind  function 

q,.  a  fuel  particle  area-to-volume  ratio, 

$($,&)  s  slope  function, 

B  3  fuel  packing  ratio. 

There  Is  no  explicit  directional  dependence  in  the  RROS  equation.  How¬ 
ever,  development  of  the  RROS  equation  was  heavily  dependent  upon  empiri¬ 
cal  data  obtained  from  wind  tunnel  simulations  of  wildland  fires.  The 
empirically  deduced  effects  of  wind  (4^)  and  slope  ($s)  are,,  thus, 
defined  to  be  for  the  down-wind  and  up-slope  directions,  respectively. 


The  RROS  Equation  3-59  predicts  that,  in  the  absence  of  wind  and 
slope,  a  point  fire  will  spread  circularly  at  the  rate  R^.  The  qualita¬ 
tive  basis  for  a  local  directional  spread  rate  is  the  assumption  that,  in 
the  presence  of  wind  and/or  slope,  a  point  fire  will  spread  elllptlcally 

^  A  ^  A  A  A 

as  a  function  of  RQ,  3  4>w,  a^t  and  4>$  3  4>s  a,,  where  and  ^  are  unit 
vectors  in  the  directions  of  down-wind  and  up-slope,  respectively. 


The  components  obtained  from  the 
3-5  In  a  two-dimensional  reference  frame 
interest.  It  is  assumed  that  the  data  upon 
apply  at  least  to  the  point  of  interest, 
slsts  of  earth-tangent-plane  coordinates  E 
the  angles  a$  and  are  azimuthal  angles, 
direction  toward  which  the  wind  Is  blowing 

A 

logical  convention).  The  unit  vector  ac 


RROS  model  are  shown  in  Figure 

with  origin  at  the  point  of 

which  R„,  and  are  based 
0  w  Ts 

The  reference  frame  used  con- 
(east)  and  N  (north),  so  that 

A 

The  unit  vector  ag  Is  In  the 
(opposite  to  common,  meteoro- 
is  the  negative  horizontal 
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projection  of  the  local  surface  normal,,  and  thus  points  In  the  up-slope 
direction. 


An  ellipse  can  be  completely  defined  in  terms  of  the  location  of 
one  focus,  the  line  of  apsides,  the  apoapsls  distance,  and  the  semi-latus 
rectum.  Define  the  point  of  interest  as  the  location  of  the  occupied 

A 

focus.  Define  the  line  of  apsides  along  a  unit  vector  a^  where 


v.  ■  <»w  +  *S)'K  +  *sl-  <3-60> 

Define  the  semi-1 atus  rectum  by 

%  3  Ro  +  |*w  x  ®mx|  +  |*s  *  *mxl  (3'61) 

Define  the  apoapsls  distance  R^  as  shown  in  Figure  3-5  by 

"mx  “  Ro  +  *  ^nx^  +  *•%  ^nx^  (3-62) 

where, 

amS  a  4  unit  vector  in  the  direction  of  maximum  spread  area, 

Rmx  =»  the  magnitude  of  maximum  spread  rate, 

Rty  *  spread  rate  in  the  direction  normal  to  maximum  spread 
rate. 

The  azimuthal  angle  of  as  shown  in  Figure  3-5  is  given  by 

on,x  *  arctan  [(^  sino^  +  sina$)/(  ^cosa^  +  *s  cos«s)  ].(3-63) 
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The  eccentricity  of  the  ellipse  is  by  definition 

e!l  =  VV  •  '3-«4> 

and  the  directional  spread  rate  in  the  azimuthal  direction  (0)  is  then 
given  by 

R(e)  =  l^(l-e)/(l  +  e  cose)  •  ’  (3-65) 

where  e  a  n  +  a  -  c^x 

An  important  characteristic  of  Equation  3-65  is  that  it  preserves  the 
magnitude  of  spread  rate  in  the  direction  implied  by  Equation  3-59. 

It  was  determined  from  comparing  computational  results  with 
actual  wildland  fire  growth  data  that  the  LDSRA  defined  by  Equation  3-65 
predicts  local  elliptical  spread  rate  patterns  exhibiting  excessive 
eccentricity.  This  was  considered  to  be  a  result  of  the  dependence  of 
Equation  3-59  upon  experimental  data  from  simulated  fires  in  laminar  wind 
fields.  Laminar  wind  flow  is  rare  in  nature,  particularly  near  the  ground 
in  rough  terrain  and  vegetation.  The  effect  of  local  turbulence  is 
considered  to  decrease  the  eccentricity  of  the  local  directional  spread 
rate  pattern  as  a  function  of  local  wind  speed.  The  eccentricity  of  the 
local  spread  rate  ellipse  was,  thus,  modified  in  accordance  with  wildland 
fire  growth  data.  The  resulting  LDSRA  is  given  by 

R(«)  -  Rmx  (1  -  e')/(l  +  «>  cosa)  '  (3-66) 

where 

e*  =  ef(u)  (3-67) 
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and 


f(u)  =  1.0  -  O.Q2u,  0  <  u  <  10  (mph)  (3-68) 

i 

f(u)  -  0.8  -  0.0054(u  -  10),  10  <  u  <  80  (mph)  .  (3-69) 

3.6  ANALYSIS  OF  THE  RROS  EQUATION 

Formulation  of  a  computational  procedure  embodying  the  RROS 
equation  for  estimating  the  growth  of  wildland  fires  revealed  problems 
related  to  input  data,  computation,  and  data  management  requirements. 
These  problems  derive  from  the  use  of  a  grid  structure  over  the  potential 
burn  area  and  the  detailed  fuel  descriptions  specified  for  the  RROS 
equation. 

The  required  cell  size  and  expected  fire  size  define  the  amount 
of  data  to  be  entered,  the  computational  detail  required,  and  the  amounts 
of  data  to  be  accessed  and  stored.  Cell  size  is  determined  by  the  most 
limiting  of  the  rate  of  change  of  terrain  slope,  wind  velocity,  fuel 
diversity  and  desired  firefront  resolution.  Consideration  of  these  param¬ 
eters  indicated  a  minimum  practical  cell  size  of  about  one  acre  for  wild¬ 
land  fires.  The  historical  data  on  wildland  fires  of  interest  implied 
that  a  potential  burn  area  of  about  104  cells  would  not  be  unreasonable. 
Reference  to  Table  3-3  for  two  fuel  categories  and  six  size  classes  shows 
that  98  fuel  parameters  plus  five  additional  parameters  are  required  as 
input  data  for  each  cell.  For  104  cells  this  is  a  very  large  amount  of 
data  and  is  comprised  largely  of  detailed  fuel  characteristics.  For  this 
reason,  the  RROS  equation  was  analyzed  to  determine  the  sensitivity  of 
spread  rate  estimates  to  each  of  the  input  data  elements. 
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An  Initial  analysis  (Reference  23)  examined  the  sensitivity  of 
RROS  equation  results  to  independent  changes  in  each  variable  over  the 
range  of  values  that  each  variable  might  exhibit.  This  analysis  parti¬ 
tioned  the  RROS  equation  variables  into  two  sets,  strong  and  weak  vari¬ 
ables.  Each  strong  variable  resulted  in  a  maximum  rate  of  change  of 
estimated  rate  of  spread  with  respect  to  the  value  of  the  subject  variable 
that  was  greater  than  ten  percent.  Weak  variables  were  below  this  thres¬ 
hold. 


The  strong  variables  were  identified  to  be  u,  o,  6,  Mf,  and 
Mx  and  the  wfeak  variables  were  W0,  h,  pp,  Sf,  Se,  and  tan  $  (see 
Table  3-2  for  variable  definitions!.  A  subsequent  analysis  examined  the 
effects  upon  RROS  equation  results  of  all  independent  variables,  all 
physically  realizable  couplings  between  strong  variables,  all  weak  vari¬ 
ables  combined  and  all  strong  variables  combined.  The  results  of  this 
analysis  are  shown  in  Table  3-4.  Examination  of  the  standard  deviation  of 
the  distribution  of  percentage  change  in  estimated  spread  rate  for  each 
tested  variable  set  reveals  that 

1.  There  are  three  strong  variables  (u,  o,  and  Mf) 

2.  The  weak  variables  Wg,  h,  pp,  Sy,  Sg,  and 
tan  $  can  be  treated  as  constants 

3.  The  two  intermediate  variables  6  and  can  be 
accomodated  with  two  leveis  of  quantization. 

As  a  result  of  these  analyses,  it  was  possible  to  reduce  the  number  of 
data  elements  required  per  cell  from  approximately  100  to  6.  All  of  the 
fuel  characteristics  were  determined  in  terms  of  the  fuel  type  within  a 
cell  and  its  age.  The  remaining  four  data  elements  consisted  of  the  wind 
velocity  components  and  the  slope  components.  This  reduction  in  the 
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amount  of  data  to  be  processed,  coupled  with  the  use  of  the  regmented 
perimeter  propagation  model  described  In  Section  4,  allowed  efficient 
simulation  of  large  wildland  fires. 
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SECTION  4 

URBAN  NASS  FIRE  BEHAVIOR  HODEL  CONCEPT 


4.1  INTRODUCTION 

No  conventional  state-of-the-a-t  exists,  in  either  a  theoretical 
or  empirical  sense,  which  suffices  to  immediately  construct  an  urban  mass 
fire  behavior  model.  Previous  Investigations  provide  a  foundation,  but  In 
themselves  are  insufficient  for  the  purpose.  Consequently,  a  model 
concept  must  have  a  high  judgemental  content  based  upon  assessment  of  the 
most  important  physical  processes  and  factors  together  with  extrapolation 
from  better  understood  situations.  The  present,  or  reasonably  short- 
range,  possibility  for  computing  or  estimating  the  major  physical 
processes  Is  a  necessity.  Simulation  of  a  complex  interactive  situation 
Is,  however,  something  more  than  straightforward  computation.  The  organi¬ 
zational  structure,  which  chooses  and  sequences  the  separate  aspects  and 
their  Interactions,  can  equal  the  physical  calculations  in  Importance. 
The  unique  model  structure  described  in  Section  5  Is  an  essential  part  of 
our  model  concept.  The  present  section  concentrates  on  physical  aspects 
of  the  model. 

Our  concept  alms  at  a  model  description  of  the  behavior  and 
propagation  of  a  mass  urban  fire  In  order  to  assess  Its  ultimate  effect 
and  to  explore  the  requisite  conditions  In  which  various  modes  of  develop¬ 
ment  (propagating  line  fire,  firestorm,  etc.)  occur.  This  emphasis  upon 
fire  development  elevates  propagation  mechanisms  to  primary  Importance. 
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We  consider  mass  fire  propagation  mechanisms,  In  order  of  spread 
potential,  to  be 

1.  Firebrand  transport  (transport  of  burning  material) 

~  0-5  km 

2.  Contiguous  spread  or  pilot  Ignition  (flame  and  hot 
gas  contact,  convection)  ~  0-50  m 

3.  Flame  radiation  ~  0-10  m 

4.  Heat  conduction  ~  cm. 

Heat  conduction  Influences  pyrolysis  rate  (burning  rate)  of  already 
Ignited  material,  but  Is  otherwise  negligible.  Flame  radiation  provides 
an  essential  feedback  of  heat  to  burning  material,  but  propagation  by  this 
mechanism  Is  rather  short  range.  It  becomes  Important  only  when  view 
factors  approach  a  solid  angle  of  4*,  as  in  closed  compartments  or  dense 
fuel  beds.  Fire  spread  between  closely  spaced  buildings  may  occur  by 
flame  radiation,  although,  in  this  case,  It  may  be  obscured  by  convective 
mechanisms.  Only  In  the  absence  of  other  propagation  mechanisms  does 
flame  radiation  become  significant.  The  important  feature  of  the  remain¬ 
ing  major  propagation  mechanisms,  branding  end  convection.  Is  their 
dependence  upon  the  details  of  the  gaseous  hydrodynamics. 

Hydrodynamic  motions  of  ambient  atmosphere  (wind  and  updrafts), 
flame  columns,  end  fire  plumes  are  determined  by  the  background  meteo¬ 
rology  as  well  as  the  state  of  the  fire,  which  Includes  the  spatial 
distribution  and  intensity  of  burning,  both  current  and  past.  In  turn, 
the  burning  state  and  Its  changes  are  determined  by  hydrodynamic  (wind  and 
convection)  effects  upon  burning  rate  and  propagation.  Despite  the  fact 
that  an^  fire  Involves  hydrodynamics,  small  scale  burning  In  the  open 
apparently  does  not  require  explicit  consideration  of  the  mutual,  two-way 
interaction  between  fire  and  wind.  The  small  scale  fire  Is  affected  by 
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the  wind,  but  the  perturbation  of  wind  by  the  fire  Is  negligible.  The 
success  of  the  Rothermel  Spread  Model  (Reference  18)  description  of  fire 
propagation  demonstrates  the  validity  of  neglecting  the  reverse  effect  of 
open  fires  upon  wind,  at  least  for  brush  fires,  which.  In  this  sense,  are 
snail  scale.  On  the  other  hand,  models  of  compartment  fires  In  still 
Intact  buildings  require  a  hydrodynamic  description  of  the  Interaction 
between  air  flow  (ventilation)  and  burning  rate.  In  the  case  of  open 
fires,  as  the  number,  Intensity,  and  areal  distribution  of  burning  regions 
grow,  the  atmospheric  perturbation  Increases  to  the  point  where  fire 
development  becomes  substantially  Influenced  by  the  fire  Itself.  Our 
concept  of  an  urban  mass  fire  model  emphasizes  this  reciprocal  Interactive 
effect,  which  we  expect  to  assume  significance  In  conflagrations  on  the 
scale  to  be  considered.  The  potential  presence  of  a  high-velocity 
fire-induced  wind  as  a  damage  mechanism  Is  also  of  Interest  In  Itself. 

In  addition  to  the  fire-hydrodynamic  reciprocal  Interaction, 
mass  fire  behavior,  of  course,  depends  upon  the  effective  (after  blast) 
Initial  Ignitions  (spatial  and  temporal  distribution)  and  other  factors, 
all  of  which  can  be  considered  fuel  parameters.  Fuel  loadings  of  burnable 
material  a*d  Its  effective  heat  content,  rates  of  combustion,  fuel 
distribution  (building  heights  and  spacing,  streets,  and  open  areas  as 
firebreaks,  vegetation,  etc.)  are  all  essentially  fuel  descriptors.  These 
parameters,  perhaps  drastically  modified  by  air  blast,  clearly  underlie 
the  behavior  of  nuclear-induced  urban  mass  fires,  and  define  the  modeled 
situation.  Recall  that  the  major  research  effort  In  producing  the 
Rothermel  model  was  quantification  of  the  effects  of  fuel  parameters, 
which  are  central  to  it.  In  the  urban  model  context  one  requires  fuel 
parameters  In  sufficient  detail  to  define  the  situation,  but  propagation 
(branding,  pilot  Ignitions)  and  burning  characteristics.  Insofar  as  they 
are  fuel  dependent,  must  be  handled  partly  parametrically  since  detailed 
Information  Is  lacking.  We  advocate  reasonable  functional  behavior  and 
Idealized  physical  models  to  estimate  these  effects.  Any  model  of  urban 
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fire  development,  whether  It  Includes  reciprocal  fire-wind  interaction  or 
not  would,  at  the  present  level  of  knowledge,  suffer  equally  from  these 
fuel-related  uncertainties. 

Once  the  reciprocal  fire-hydrodynamic  (fire-wind,  for  short) 
interaction  is  selected  as  the  central  aspect  of  urban  mass-fire  behavior, 
we  must  consider  methods  for  its  calculation  and  exploitation.  Our  con¬ 
cept  is  essentially  the  choice  of  an  approach  to  this  problem. 

4.2  MODEL  ALTERNATIVES 

Mass  fire  development,  considered  as  an  interactive  fire-wind 
problem,  demands  a  three-dimensional  description  of  the  flow  of  hot  gases 
and  ambient  air.  Although  heat  sources  are  localized  on  or  near  the  sur¬ 
face,  brand  transport  depends  upon  velocities  aloft,  and  winds  near  the 
surface— important  for  pilot  ignltions—are  only  a  special  part  of  the 
unified  total  flow.  Later  we  will  examine  some  evidence  from  the  Flambeau 
tests  which  indicates  the  need  for  three-dimensional  resolution  even  in 
the  near-surface  region  of  heat  production  and  active  flaming. 

4.2.1  Detailed  Hydrodynamic  Calculation 

A  detailed  3-d  hydrodynamic  calculation,  based  upon  fundamental 
equations  of  motion,  could  in  principle  be  used  to  determine  the  flow 
field.  A  major  practical  problem  with  this  approach  Is  the  need  for  high 
spatial  resolution,  sufficient  to  resolve  irregularities  and  discontinui¬ 
ties  in  zones  of  burning  and  convection,  which  are  certainly  present  at 
Ignition  and  immediately  after,  and  which  may  persist  for  extended  per¬ 
iods.  Such  calculations  of  detailed  3-d  flow  fields  would  be  expensive 
and  lengthy  to  perform,  even  if  sufficient  memory  capacity  for  the  com¬ 
puter  calculations  were  available. 
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A  physical  problem  with  this  approach  arises  In  the  need  to  com¬ 
pute  the  effects  of  turbulence,  which.  In  this  case,  Is  largely  produced 
by  buoyancy  of  hot  gases.  It  is  not  clear  whether  there  Is  an  expeditious 
method  for  calculating  the  behavior  of  buoyant,  turbulent  flows  directly 
from  fundamental  principles.  Second-order  closure  schemes  of  turbulence 
modeling  are  far  too  demanding  to  be  practical  in  a  model  which,  in 
addition,  performs  many  other  types  of  calculation  and  which  must  be 
repetitively  run  to  assess  the  effects  of  varying  conditions.  At  this 
level  of  detail,  turbulent  exchange  approximations  based  upon  mixing 
length  hypotheses  do  not  alleviate  the  computational  burden  significantly 
and  of  themselves  introduce  large  uncertainties. 

Numerical  computation  of  detailed  flow  fields  at  high  resolution 
Is  perhaps  a  last  resort,  and  other  avenues  should  be  explored  first. 

4.2.2  Low-Resolution  (Lumped-Parameter)  Hydrodynamic  Calculations 

At  the  other  extreme,  we  might  consider  only  a  global,  compara¬ 
tively  low  resolution  description  of  the  flow  field.  Previous  approaches 
to  the  mass  fire  problem  have  calculated  the  axialiy  symmetric  flow  field 
resulting  from  a  single  large  area  circular  fire  (with  assumed  smoothed 
and  uniform  heat  production  rates  per  unit  surface  area).  By  invoking 
either  mixing  length  assumptions  or  turbulent  entrainment  assumptions,  the 
calculations  have  produced  reasonable  solutions  for  the  single  large  fire 
plume  in  the  asymtotic  region  at  large  altitudes  above  the  source.  These 
calculations,  however,  typically  yield  unrealistic  behavior  In  the  region 
near  the  ground.  One  effect  (Reference  12)  Is  rapid  convergence  of  the 
rising  hot  gases  to  the  center  of  the  burning  area  to  form  a  narrow  plume, 
which  then  expands  normally  In  its  subsequent  rise.  This  exaggerated 
behavior  is  not  necessarily  inherent  in  the  model  type,  being,  we  believe, 
a  consequence  of  oversimple  assumptions.  The  single  plume  formed  at  high 
altitudes  by  coalescence  of  hot  turbulent  gases  from  many  separate  heat 


sources  may  be  adequately  represented  by  this  type  calculation.  The 
coalescence  process  Itself  and  the  detail  of  wind  fields  at  lower  eleva¬ 
tions  resulting  from  discrete  burning  zones  requires  something  In 
addition,  or  perhaps  a  different  approach.  Though  potentially  useful  for 
determining  large  scale  atmospheric  perturbations  and  overall  flow  fields, 
these  low  resolution  hydrodynamic  models  fall  short  of  providing  the  cru¬ 
cial  small  scale  structure  needed  to  model  fire  development. 

4.2.3  Parameterized  Multiple  Plumes 

In  the  face  of  the  limited  understanding  of  mass  fire  phenomena 
and  the  recognized  inadequacy  of  previous  approaches  to  the  problem,  we 
have  evolved  a  concept  for  modeling  fire  development,  which  does  not 
attempt  a  first  principles  calculation  of  the  hydrodynamics.  Instead,  the 
flow  fields  are  deduced  as  a  superposition  of  velocities  resulting  from 
the  actual  distribution  of  segmented  burning  zones  and  their  attendant 
entraining  plumes.  This  approach  requires  that  the  major  features  of 
individual  burning  zones  and  fire  plumes  be  parameterized  in  terms  of  com¬ 
paratively  few  variables:  heat  production  rate  per  unit  area,  areal 
extent,  and  wind  description  being  the  principal  variables.  Heavy 
reliance  is  placed  upon  similarity  solutions  together  with  laboratory  and 
field  experiments,  supplemented  by  physical  arguments  and  scaled  to  mass 
fire  conditions.  In  situations  where  this  procedure  is  considered 
inadequate,  parameterized  functional  relations,  which  summarize  results  of 
off-line  hydrocode  calculations  of  single  plumes,  could  be  employed. 
Essentially  each  burning  zone  (or  subdivision)  is  represented  by  sources 
and  sinks  of  mass  distributed  in  three  dimensions.  From  this  distribu¬ 
tion,  atmospheric  flow  external  to  plumes  and  flames  can  be  calculated. 
The  external  flow  reacts  back  upon  the  plumes  to  establish  the  location  of 
the  source-sink  distribution  in  a  self-consistent  manner.  Once  estab¬ 
lished,  the  wind  field  determines  burning  rates  and  propagation  to  advance 
the  fire  development. 
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The  approach  Is  not  without  difficulties  and  leads  to  several 
problem  areas,  which  require  Investigation.  However,  it  exploits  experi¬ 
mental  measurements  and  qualitative  insights  obtained  In  previous 
research,  avoids  the  inherent  difficulty  of  turbulence  calculations,  and 
can  achieve  detail  and  resolution  without  excessive  computation. 

By  keying  the  description  to  a  simple  demonstration  model,  the 
following  sections  illustrate  the  model  concept  and  clarify  the  procedure 
to  Implement  it.  Problem  areas  are  discussed  as  they  are  encountered. 

4.3  ATMOSPHERIC  MOTION 

First,  we  briefly  sketch  our  conception  of  the  relevant  phenome¬ 
nology  of  fires  burning  in  the  open  as  derived  from  previous  research  and 
observations.  From  the  standpoint  of  urban  mass  fires,  all  detailed  tests 
have  been  small  fires.  Leaving  aside,  for  the  moment,  the  effects  of 
atmospheric  stratification  and  scale  height,  these  small  fires  are 
governed  by  the  same  physical  parameters  as  larger  ones.  The  effect  of 
larger  distance  scales  is  quantitative,  not  qualitative,  and  we  presume 
relations  between  nondlmensioned  groups  of  variables  may  be  used  In  large 
scale  as  well  as  small  scale  fires. 

Figures  4-l(a)  and  4-l(b)  picture  the  features  of  a  small  area 
laboratory  fire.  Vertically,  the  fire  consists  of  a  flame  zone,  a  zone  of 
Intermittent  flaming,  and  a  buoyant  plume,  which  quickly  assumes  the 
asymptotic  behavior  of  a  weakly  buoyant  plume,  since  mixing  with  cool 
ambient  air  is  rapid.  At  small  burning  rates  Q  (J/s),  the  flame  column 
shows  convergence,  but  at  higher  burning  rates  the  flame  column  diverges 
almost  immediately  as  It  rises.  This  distinct  behavior  must  depend  upon 
both  Q  and  the  fire  area  A;  but  as  long  as  the  area  is  small  enough,  the 
burning  takes  place  in  a  simple  flame  column.  The  work  of  McCaffrey 
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Figure  4-1 (a)  and  (b).  Small  area  fire 
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(Inference  24)  and  Cox  and  Chltty  (Reference  25)  has  provided  reliable 
parameterization?-  of  those  single  fires,  but  the  area  effect  awaits  Inves¬ 
tigation. 

When  several  small  area  fires  are  burning  In  proximity,  one 
expects  the  situation  shown  In  Figure  4-2.  Each  separate  fire  and  Its 
plume  entrains  air.  The  effective  "sinks"  of  all  plumes  produce  an  air 
velocity,  which  tips  the  plumes  Inward  and  brings  them  together  at  a 
coalescence  height,  after  which  the  buoyant  gases  form  a  single  plume.  We 
encounter  here  a  fire- Induced  wind,  strongest  at  the  location  of  the  peri¬ 
pheral  fires,  which  may  Influence  their  spread. 

If  all  the  small  area  fires  coalesce,  or  If  a  large  area  is 
burning  without  gaps,  the  situation  approaches  that,  illustrated  in  Figure 
4-3.  In  this  case,  according  to  the  observation  of  Countryman  (Reference 
9)  in  the  Flambeau  tests,  a  single  flame  column  does  not  persist. 
Instead,  the  fire  spontaneously  divides  Into  separate  flame  columns  and 
plumes,  which  maintain  their  identity  through  a  short  transition  zone 
before  combining  Into  a  simple  plume.  Although  detailed  analysis  of  the 
reason  for  this  behavior  In  lacking,  it  clearly  Is  a  manifestation  of  the 
dynamic  instability  of  a  single  flame  column  in  certain  conditions.  It 
may  be  related  to  the  Raleigh-Taylor  buoyancy  instability,  or  to  the 
requirement  of  an  air  supply  for  burning,  which  may  favor  the  large  sur¬ 
face/volume  ratio  of  multiple  flame  columns.  Here  again,  the  fire  charac- 
ter  must  depend  not  only  upon  the  total  Q  but  also  upon  the  fire  area, 
A.  This  behavior  has  a  strong  effect  upon  fire-induced  winds,  since  a 
large  amount  of  entrainment  must  occur  in  the  transition  and  flame  zones 
below  the  base  of  the  combined  plume.  Wind  estimates  depend  upon  resolu¬ 
tion  of  these  features. 

The  configuration  of  Figure  4-3  is  similar  to  that  of  Figure 
4-2,  but  the  separate  columns  have  different  origins.  In  the  one  case, 
they  are  imposed  by  physical  separation  of  the  fires,  ami  in  the  other  are 
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Figure  4-2.  Separated  small  area  fires. 


3.  Large  area  fire  (spontaneous  break -up  into  separate  flame 
columns).  "Countryman  Descriptive  Model." 
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a  result  of  fire  dynamic".  In  either  case,  or  In  a  possible  combination 
of  the  two,  the  atmospheric  motion  and  plume  character  are  interactively 
determined.  Plumes  as  effective  sinks  of  ambient  air  determine  fire 
winds,  but  winds  affect  tUe  plume  (or  sink)  distribution. 


The  illustrations  of  Figures  4-2  through  4-4  depict  some  of  the 
effects  of  atmospheric  fire  plumes  in  the  absence  of  an  ambient  wind. 
Both  ambient  winds  and  atmospheric  stratification  can  modify  the  picture. 
Since  we  have  by  no  means  solved  these  problems,  we  describe  in  the  fol¬ 
lowing  paragraphs  some  partial  features  to  illustrate  the  concept.  Most 
of  these  features  have  been  incorporated  as  elements  in  the  interactive 
demonstration  model. 

4.3.1  Single  Fire  Plumes 

A  single  weakly  buoyant  plume,  in  a  uniform  atmosphere  and  in 
the  absence  of  wind,  approximates  the  form  of  a  right  circular  cone.  Tak¬ 
ing  the  apex  of  the  cone  at  the  ground,  the  similarity  solution  of  the 
buoyant  plume  (References  26  and  27)  requires 


w  =  a  T 


ZB  =  n  r8  , 


(4-1) 


where  w  is  the  mean  upward  velocity  in  the  plume  and  the  ratio  Zg/rg 
between  the  height  and  radius  of  the  plume  has  the  constant  empirical 
value,  n  =  5.51.  The  constant,  a,  is  a  function  of  the  buoyancy  flux  in 
the  plume,  which  is  directly  related  to  the  heat  production  rate  Q  if  the 
plume  is  fire  driven.  By  assuming  3w  s  w,„ax,  the  plume  centerline 
velocity,  the  plume  measurements  of  Cox  and  Chitty  (Reference  25)  permit 
computation  of  the  relation 
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Figure  4-4.  Plu»e  geometry  and  sinks. 
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a  -  3.897(10)-2  Q1/3 


(4-2) 


where  units  are  tries  and  Q  Is  heat  production  In  watts.  An  entrainment 
constant  a  (defined  as  the  constant  of  proportionality  between  the  mean 
upward  plume  velocity  and  the  indraft  velocity  of  ambient  air  through  the 
sides  of  the  plume),  u  *  o  w,  can  be  calculated  from  these  relations,  by 
using  mass  conservation  and  the  assumption  of  Incompressible  flow.  The 
result,  a  *  5/(6n)  3  0.151,  does  not  agree  with  measured  values,  which  are 
near  0.1  .  Vie  use  a  3  0.1. 

Outside  the  plume,  the  ambient  air  motion  may  be  modeled  as  the 
flow  Induced  by  a  distribution  of  sinks  located  along  the  plume  axis.  The 
entrainment  sink  density  X  Is 


x(Z')  3  *  -2waa(Z' )2/3/n 

dZ 

where  V  is  the  vertical  coordinate  of  the  source  point  and  V  (m3/s)  is 
volume  flow  rate.  In  addition,  we  account  for  the  air  supply  necessary 
for  combustion  by  placing  a  point  sink  of  magnitude  Va  on  the  ground  at 
the  cone  apex: 


Va  3  -k  ^(hoPo)  - -5.5(i°)-7  Q 

Here  h0  is  the  heat  of  combustion  of  the  fuel  (incomplete  combustion  of 
wood,  -1(10) 7  0/kg),  k  is  the  mass  of  air  required  to  burn  a  unit  mass  of 
fuel  (~7-10  kg-alr/kg-fuel)  and  ■  1.27  kg/m3,  the  ambient  air 
density.  Q  Is  the  burning  Intensity  In  watts. 

The  plume  geometry  and  sink  distribution  are  shown  In  Figure 
4-4,  where  in  order  to  fit  the  boundary  condition  at  the  ground  (zero 
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vertical  velocity),  a  below-ground  Image  distribution  Is  added.  This  dou¬ 
bles  Va,  the  magnitude  of  the  point  sink  at  the  ground. 

Before  continuing,  we  note  that  this  model  Includes  only  sinks, 
although  one  would  expect  the  effect  of  heat  addition  to  be  gas  expansion 
and,  hence,  a  net  overall  source.  The  source  is  present  at  the  plume  cap 
in  a  developing  plume,  and  must  be  considered  while  the  plume  is  being 
established.  Sudden,  quasi-explosive  increases  of  burning  rate  can  also 
produce  transient  outdrafts  by  this  mechanism.  Here  we  consider  a  steady, 
established  plume,  place  tne  plume  cap  at  infinite  height;  and  neglect  its 
effect. 

The  velocity  field  external  to  the  plume  is  given  by  integration 
over  the  sink  densities.  In  the  cylindrical  coordinates  chosen,  the 
results  take  the  form 


vr(r,e)  ■  Va/[2,rr2(l+e2)3/2]  -  aafp(e)/[2  nr1/3] 


(4-3) 

vz(r,B)  »  Vae/[2irr2(l+02)3/2]  +  oafz(p)/[2  nr1/3] 

where  p  3  Z/r,  and  units  are  mks  throughout.  The  nondimensional  functions 
fr(p)  and  fz(8),  obtained  by  integration  of  the  plume  sinks  between  V 
■  ±  »,  are  shown  in  Figure  4-5.  Values  of  p  are  restricted  to  the  inter¬ 
val,  0  <  p  <  5.51,  the  region  outside  the  plume.  At  fixed  p,  the  atmo¬ 
spheric  velocities  induced  by  plume  entrainment  fall  off  at  the  slow  rate 
of  r-1/3,  a  result  of  the  infinite  plume  length.  The  horizontal  veloci¬ 
ties  directed  radially  inward  always  exceed  vertical  velocities  except 
very  near  the  origin  where  effects  of  the  point  sink  dominate.  Throughout 
most  of  the  external  region,  turbulent  entrainment  by  the  plume  is  a  far 
stronger  driver  of  the  flow  than  is  the  sink  of  combustion  air. 
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This  analysis,  though  very  simple.  Illustrates  the  technique  of 
calculating  Induced  atmospheric  motion  from  plume  parameters.  In  this 
case,  the  sole  free  parameter  Is  Q,  the  power  of  the  fire.  Clearly, 
atmospheric  stratification,  which  not  only  modifies  the  plume  but  also,  In 
stable  conditions,  restricts  vertical  flow  In  the  ambient  atmosphere, 
would  alter  these  results.  Plume  bending  In  a  wind,  either  ambient  or 
induced  by  other  plumes,  would  change  both  the  magnitude  ind  location  of 
sinks.  These  questions  need  to  be  addressed  In  a  complete  model. 

A  simple  algorithm  to  demonstrate  our  conceptual  approach  spe¬ 
cializes  the  above  analysis  to  calculate  ground  level  wind.  Each  separate 
fire  plume  acts  as  a  partial  source  for  the  ground  level  wind  and  contrib¬ 
utes  a  velocity, 

vr(r)  ■  Va/(2irr2)  -  0.862  aa/fnr1/3)  ,  (4-4) 

directed  toward  the  plume  origin.  The  numerical  factor  0.862  Is  half  the 
value  of  fr(p)  at  p  *  0.  This  formula  applies  to  a  developed  plume  of 
infinite  height.  To  account  for  the  plume  formation  transient,  Important 
in  the  initial  stages,  we  extend  the  Integral  of  plume  entrainment  sinks 
only  over  the  actual  vertical  height  of  the  plume.  (The  source  at  the  cap 
Is  still  neglected  for  simplicity.)  The  height  of  the  plume  cap  for  a 
steady  burning  rate  starting  at  ta0  Is  approximately, 

1  -  0.234  qVV/4 
c 

again  in  mks  units.  A  numerical  integration  over  plume  sinks  to  this 
height  introduces  a  function  of  2^/r  *  Y  place  of  the  constant, 

0.862.  This  function  Is  fitted  by  the  approximate  formula. 


(4-5) 


0.862  ♦  g(-y)  -  0.862  y1*55/ (1.328  +  y1*-59) 

which  plays  the  role  of  a  build-up  function  lor  the  plume- Induced  surface 
wind  fields.  Close  to  the  plume,  winds  build  up  to  their  final  steady 
value  more  quickly  than  ut  larger  distances. 

4.3.2  Multiple  Fire  Plumes  ondl  k  ea  Fires 

If  a  mass  urban  fire  were  to  consist  of  separate  point  f l^es  of 
prescribed  burning  rate  whose  plumes  remained  vertical,  separate,  and 
independent,  the  fire- Induced  ground  level  wind  could  be  obtained  as  a 
superposition  of  contributions  from  each  fire  given  by  the  formulas  of  the 
previous  section.  This,  perhaps  zeroth  order,  approximation  Is  In  fact 
the  general  method  used  to  calculate  wind  fields  In  our  demonstration 
model.  However,  the  model  does  not  operate  with  point,  fires,  but  with 
square  grid  cells  (modeling  city  blocks),  and  carries  the  total  burning 
rate  of  each  cell.  A  "fire"  is  defined  as  a  group  of  contiguous  burning 
cells,  separated  from  other  fires  by  unignited  areas.  In  turn,  each  fire 
may  be  divided,  on  a  geometrical  basis,  into  "sub-fires,"  the  two 
horizontal  dimensions  of  which  are  approximately  equal.  This  procedure 
permits  computation  of  the  fire-  induced  wind  by  a  superposition  of  effects 
from  each  "sub-fire."  Outsido  the  sub-fire.  Its  effect  is  that  of  an 
equivalent  plume  at  its  center,  and  the  shape  of  each  sub-fire  may  be 
approximated  as  a  circle  of  radius,  rf. 

In  order  to  perform  the  wind  calculation,  the  number  of  plumes 
or,  what  Is  the  same  thing,  the  number  of  equivalent  point  fires  in  the 
burning  area  of  the  sub-fire  must  be  determined.  This  Is  the  problem 
mentioned  in  Section  4.3  and  illustrated  In  Figures  4-2  and  4-3  relating 
to  the  spontaneous  breakup  of  an  area  fire  Into  separate  flame  columns  and 
plumes.  Its  Importance  stems  from  the  fact  that  air  entrainment  Into  a 
plume  is  governed  by  the  factor  a  defined  In  Equation  4-2,  which  Is 
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proportional  to  Q1/ 3.  A  given  total  Q  of  an  area  fire,  divided  Into  N 
f'lame  plumes  entrains  at  a  rate  proportional  to  NCQ/N)1/3.  Entrainment 
is,  thus,  proportional  to  N2/3.  Many  separate  plumes  produce  more 
atmospheric  perturbation  (fire-induced  winds)  than  does  a  single  fire  of 
the  same  total  power.  This  reasoning  supports  the  contention  that 
resolution  of  separate  fires  and  burning  areas  is  essential  In  a  mass  fire 
model.  On  the  basis  of  these  Ideas,  multiple  separated  ignitions  may  be 
more  likely  to  generate  firestorm  conditions  than  area  burning. 

Our  analysis  of  these  features  has  hardly  begun.'  To  implement 
the  demonstration  model  we  estimate  the  number  of  plumes,  N,  present  in  a 
sub-fire  of  radius  rf  with  a  total  power  Q  as 

N  -  4.94(10) 3  Q-2/3rJ/3  .  (4-6) 

This  expression  has  only  scant  justification,  based  upon  inexact  estimates, 
of  flow  conditions  required  to  minimize  pressure  gradients,  which  we  do 
not  reproduce  here.  It  Is  expected  that  the  number  of  plumes  (or  flame 
columns)  increase  with  fire  dimensions,  and  decrease  with  burning  rate. 
No  obviously  ridiculous  results  occur  within  the  ranges  of  the  variables 
employed. 


Each  sub-fire  Is  assumed  to  consist  of  the  number  of  plumes 
given  by  Equation  4-6,  which  are  uniformly  distributed  over  the  quasi¬ 
circular  area  of  the  sub-fire.  The  effects  of  each  plume,  given  by  terms 
c'  the  form  of  Equation  4-4,  are  superimposed.  The  total  combusion  air 
sink  Va  Is  placed  at  the  fire  center,  and  the  weak  r1/3  dependence  upon 
distance  away  from  the  plume  center  Is  neglected  Inside  the  sub-fire. 
There  results  a  simple  formula,  used  in  the  demonstration  model,  for  the 
contribution  of  each  sub-fire  to  the  total  fire-induced  surface  wind: 
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for  r  >  rf  vf(r)  »  -Va/(2*r2)  -  g< y)Q“  3 

for  r  <  r^  v^r)  -  -Var/(2mr^)  -  g(  YiQ"1/9^0/9*" 1  (4-7) 


where  Q,  Va,  and  r^  are  the  respective  total  power,  combustion  air 
flow,  and  radius  of  the  sub-fire.  The  distance  r  is  the  radial  distance 
from  the  sub-fire  center  and  the  build-up  function  g(y)  is  given  by 
Equation  4-5  computed  on  the  basis  of  the  total  Q  of  the  largest  fire 
(not  sub-fire)  in  the  grid.  The  total  fire-wind  at  any  point  is  the 
vector  sun  of  these  contributions,  each  of  which  Is  directed  toward  the 
center  of  its  sub-fire  source. 

The  fire-induced  wind  at  ground  level  obtained  this  way  is  an 
overestimate,  especially  in  the  region  outside  of  sub-fires,  since,  in 
reality,  separate  plumes  coalesce  into  a  single  plume  whose  entrainment  Is 
less,  (See  Figures  4-2  and  4-3.)  The  far  outside  region  is  influenced 
more  by  the  single  plume  at  high  elevations  than  by  the  multiple  plumes 
lower  down.  The  reverse  effect  of  the  winds  upon  plume  location  and 
entrainment  rate  is  also  neglected.  By  contrast,  very  near  to  and  inside 
of  sub-fires,  the  winds  may  be  more  realistic.  If  an  ambler.t  wind  Is 
present,  its  velocity  Is  simply  added  vectorially  to  the  fire-induced 
wind. 


In  general,  the  demonstration  model  neglects  the  essential 
interaction  between  fire  plumes,  as  well  as  the  effects  of  plume  tilt  and 
atmospheric  stratification.  These  matters  require  further  analysis,  but 
from  the  foregoing,  the  necessity  of  resolving  these  details  and  the  effi¬ 
cacy  of  this  approach  may  be  appreciated. 
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4.4  PROPAGATION  BY  BRANDING 

Winds  are  important  because  of  their  influence  upon  propagation 
mechanisms.  In  this  section,  we  develop  a  simplified  branding  algorithm, 
used  In  the  demonstration  model,  which  Illustrates  the  problen  areas  and 
general  features  of  fire  propagation  by  branding. 

4.4.1  Brand  Transport 

According  to  the  experiments  of  Tarifa  (Reference  28),  the  ver¬ 
tical  component  of  firebrand  motion  can  be  described  by  a  terminal 
velocity,  with  respect  to  surrounding  air,  which  is  a  function  of  time. 
We  approximate  the  families  of  curves  obtained  In  his  experiments  by  a 
simple  linear  function  for  the  terminal  velocity,  wj-, 

w-j  =>  w.|  -  kt 

where  wj-  Is  parameterized  by  w-j,  the  initial  terminal  velocity  of  the 

brand,  which  expresses  the  effect  of  brand  size  and  geometry.  The 

constant,  k  ~  1/12  m/s^,  is  the  same  for  all  brands.  At  tg  *  w^/k, 

the  terminal  velocity  is  zero,  which  is  taken  as  the  time  of  brand 

burnout. 


The  updraft  velocity  of  a  fire  plume  reaches  its  maximum,  w0, 
in  the  region  of  intermittent  flaming  at  a  height,  Zq,  after  which  the 
velocity  decreases  as  w  a  aZ"1/3.  The  values  of  a,  Zq,  and  w0  are 
determined  by  the  power  of  the  fire,  Q,  using  the  measurements  of  Cox  and 
Chitty  (Reference  25).  We  assume  these  relations  to  apply  for  the  verti¬ 
cal  component  of  plume  motion,  even  in  the  presence  of  a  wind.  In  a  wind, 
which  is  needed  for  brand  transport,  both  the  horizontal  plume  and  brand 
velocity  are  taken  as  equal  to  the  wind  velocity.  These  plume  velocities 
are  used  as  the  mechanism  for  brand  lifting.  The  potentially  important- 
brand  lifting  mechanism  of  strong  fire  whirls  is  neglected. 
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In  a  plume  the  vertical  velocity  of  a  brand  Is, 

Z  -  aZ"1/3  -  (v*j  -  kt);  Z  s  ZQ  at  t  «  0  .  (4-8) 

Me  take  the  condition  Z  s  0  as  the  condition  for  the  brand  to  fall  out  of 

the  plume.  An  integration  of  this  equation  would  determine  the  height  and 
time  at  which  a  brand  leaves  the  plume  or  whether  it  burns  out  before 

leaving  the  plume.  Clearly,  small  brands  will  burn  out  inside  the  plume 
before  the  condition  is  met.  An  upper  limit  of  brand  size  is  fixed  by  the 
lofting  condition,  w-j  <  w0.  The  fire  parameters  determine  a  range  of 

relevant  initial  terminal  velocities,  <  w^  <  w0,  of  brands 

that  can  be  both  lofted  and  fall  out  while  still  burning. 

For  expediency,  we  do  not  perform  the  Integration.  Instead,  a 
series  of  linear  approximations  are  employed  to  estimate  the  time  of  fall¬ 
out  and  the  lofting  height  of  the  brands.  We  characterize  a  brand  by  the 

parameter  n  defined  in  terms  of  its  initial  terminal  velocity  and  the 

range  of  relevant  terminal  velocities: 

n  =  (wQ  -  wi)/(w0  -  wm.n)  ;  0  <  n  <  1  . 

A  brand  will  stay  in  the  plume  for  a  time  tp  a  ntyv  where  ty  is  the 

time  of  fallout  of  the  smallest  brand  that  falls  out.  The  height  of 

fallout  Zp  is  given  by  linear  interpolation  of  Z"1/3,  on  the  height 
Zy  to  which  the  smallest  brand  rises, 

Zpi/3  -  (1  -  n)  Z^3  +  nZ^3  . 

When  a  brand  leaves  the  plume,  its  initial  terminal  velocity  for  vertical 
fall  through  the  air  is  «'  s  xj  -  ktp,  and  for  the  period  of  its  fall 


to  the  ground,  its  vertical  velocity  is  Z  ■  «(w'  -  kt).  The  time 

interval  required  to  fall  to  the  ground  from  the  height  Zp  is 

it  •  [i .  a  -  f)  >/»]  , 

where  f  a  2kZp/(w')2.  Provided  f  <  1,  a  real  solution  exists  and  the 
brand  reaches  the  ground  before  burnout.  For  the  approximations  chosen 
below,  this  condition  is  always  met  for  the  range  of  brand  sizes,  which 
fall  out  after  lofting. 

The  total  flight  time  of  a  brand,  tG  *  tp  +  at,  is  less  than 
Its  burning  time,  tB.  Our  simplified  algorithm  assumes  the  horizontal 
transport  of  a  brand  takes  place  in  a  uniform  wind,  V,  the  value  of  the 
surface  wind  field  at  the  point  of  brand  origin.  (The  result  of  this 
approximation  can  be  grossly  unrealistic  in  some  cases.  One  needs  to 
determine  the  brand  trajectory  with  the  complete  3-d  wind  field.)  With 
this  simplification,  the  brand  reaches  the  ground  at  a  point  displaced  a 
distance  D  =  V  tg  from  its  point  of  production. 

One  would  expect  the  probability  for  a  grounded  firebrand  to 
produce  a  new  ignition  to  depend  upon  the  fuel  present  and  the  brand's 
remaining  burning  time,  r  *  tB  -  tG.  This  clearly  stochastic  event  Is 
modeled  with  the  ignition  probability 

Pj  a  Q.3[l  -  exp(r/60)]  , 

where  times  are  in  seconds.  We  have  found  no  data  on  this  process,  so 
this  arbitrary  expression  is  justified  only  by  its  reasonableness.  In 
fact,  it  likely  overestimates  ignition  probability. 
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4.4.2 


Fire  Parameters  and  Brand  Production 


To  Implement  the  transport  calculation  described  above,  the 
values  of  several  quantities  must  be  determined  from  the  fire  characteris¬ 
tics.  The  maximum  updraft  velocity,  w^,,  of  each  flame  column  Is  the 
decisive  quantity  for  brand  lofting.  The  formula.  Equation  4-6,  for  the 
number  of  flame  columns  in  a  area  fire  Is  employed  to  determine  the  burn¬ 
ing  rate  of  each  flame  column  In  a  single  burning  cell.  (We  use  the 
burning  cells  of  the  model,  Instead  of  sub-fires,  to  compute  brand  trans¬ 
port.)  The  plume  parameters  of  each  flame  column  turn  out  to  depend  upon 
the  ratio  q  *  (Qc/rc) l/3,  where  Qc  Is  the  power  of  the  cell 
burning  and  rc  Is  the  radius  of  a  circle  whose  area  is  equal  to  the  area 
of  the  square  cell.  From  the  numerical  values  measured  by  Cox  and  Chitty 
(Reference  25),  we  compute 


wQ  *  9.17(10)-2  q 
ZQ  -  4.20(10)"4  q2 


(4-9) 


in  mks  units.  On  the  basis  of  a  rough  approximation  to  the  solution  of 
Equation  4-8,  the  brand  related  quantities  are 

wmin  “  4‘74™-2  < 

ZT  -  5.82(10)" 3  q  (4-10) 

tj.  -  4.44(10)" 1  q 

again  in  mks  units. 
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The  transport  of  a  brand  with  parameter,  »  produced  In  a 
burning  cell  of  given  q  can  be  determined  with  the  sequence  of  calcula¬ 
tions  outlined  In  the  preceding  section.  The  results  of  the  algebra  are 


tQ  -  q(0.444*y  +  1.100f£  -  0.975tyf£) 
tB  -  q(1.100  -  0.531ry)  , 


(4-11) 


where 

fK  =  4.726(10)“‘*(1  -  0.887^)- 2Cl  -  O^ry)*3 
and 


ffc  ■  1  -  (1  -  fK)1/2 


There  remains  the  critical  question  of  how  many,  and  what  size, 
brands  are  produced  per  unit  time  by  the  burning  cell.  Since  little  Is 
known  of  this,  we  assume  the  number  of  effective  brands  produced  in  the 
time  Interval,  T,  is  proportional  to  T  and  to  the  burning  rate, 

N  -  10- 11  Q  T 
c 

where  the  constant  of  proportionality  Is  frankly  arbitrary,  but  which.  In 
reality,  must  be  fuel  dependent.  A  brand  parameter  n  “  1  corresponds  to 
the  smallest  relevant  brand.  Small  brands  are  assumed  to  be  more  likely 
than  large  ones,  so  a  linear  probability  distribution  of  n  Is  postulated: 

P(n)dn  *  2ndn  ;  0  <  n  <  1 
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During  a  time  Interval  T,  the  r>ode1  selects  N  values  of  from  this 
distribution  for  each  burning  cell  and  computes  the  transport  of  each 
brand.  If  the  brand  grounds  in  unignited  fuel,  Ignition  at  this  location 
either  occurs  or  does  not,  as  determined  on  the  basis  of  the  ignition 
probability.  Should  the  computed  number  of  brands  N  be  less  than  one,  the 
number  Is  Interpreted  as  the  probability  for  producing  one  brand  in  the 
chosen  time  Interval. 

4.4.3  Branding  Discussion 

Stochastic  elements  play  a  large  role  in  fire  propagation  by 
branding.  Their  origin  is  in  part  physical  (random  brand  shape  and  size, 
turbulent  fluctuations),  but  they  chiefly  originate  from  areas  of  igno¬ 
rance.  To  a  good  approximation,  given  adequate  hydrodynamics,  brand 
transport  can  be  modeled  deterministically.  The  simple  algorithms  devel¬ 
oped  above  demonstrate  this.  However,  both  ends  of  the  branding  process 
encounter  more  serious  areas  of  Ignorance— brand  production  estimation  and 
brand  ignition  estimation.  Limits  of  resolution  of  fuel  properties  and 
actual  lack  of  knowledge  will  likely  conspire  to  perpetuate  this  situa¬ 
tion.  Any  branding  model  we  can  imagine  will  have  a  significant 
probabilistic  content. 

These  considerations  make  it  appropriate  to  emphasize  at  this 
point  some  modeling  philosophy,  which  pertains  to  branding,  but  which  also 
applies  generally  to  our  total  model  concept.  In  the  first  place,  proba¬ 
bilistic  elements  of  a  model  do  not  always  Imply  that  Its  overall  behavior 
is  random.  Provided  the  basic  underlying  probability  distributions  are 
narrow  enough  or  provided  choices  from  broad  probability  distributions  are 
made  often  enough,  the  overall  aggregate  results  will  dia.olay  only  a  fin¬ 
ite  and  acceptable  scatter.  In  effect,  the  results  are  deterministic. 


In  the  second  place,  the  bas*c  probability  distributions  govern¬ 
ing  modeling  choices  and  logic  ran  themselves  be  conditioned  and  changed 
by  current  values  of  physical  parameters.  Partial  and  qualitative  know¬ 
ledge,  even  Intuition,  can  In  this  manner  be  Incorporated  with  effect, 
meaning,  and  Interactive  potential.  Absence  of  a  precise,  deterministic 
description  of  a  process  Is  no  license  to  model  It  by  a  rigid  and  Immut¬ 
able  random  choice. 

4.S  CONTIGUOUS  FIRE  SPREAD  (CONTAGION) 

Contiguous  fire  spread— local  spread  from  ignited  fuel  to  adja¬ 
cent  unignited  fuel— Is  a  possible  spread  mechanism  whether  branding 
occurs  or  not.  In  the  context  of  wildland  fires,  the  Rothermel  model 
successfully  models  spread  In  continuous  fuel  beds  on  the  basis  of  local 
fuel  parameters,  together  with  the  local  terrain  slope  and  wind.  The 
physical  framework  of  the  Rothermel  model  Is  a  spread  rate  determined  by 
the  ratio  of  emitted  heat  flux  to  heat  required  to  Ignite  new  fuel,  but 
the  numerical  implementation  of  this  Idea  required  extensive  tests  to 
quantify  fuel  characteristics. 

In  the  urban  mass-fire  context,  fuel  types  and  distributions  are 
probably  less  well  quantified  than  In  wildlands  and  are  certainly  dis¬ 
tinctly  different.  Also,  building  separation  and  the  effect  of  streets  as 
firs, breaks  may  control  spread  (even  after  possible  fuel  scattering  and 
modification  by  blast  effects).  Because  of  these  distinctions,  the  short- 
range  fire  propagation  across  firebreaks  Is  judged  the  most  Important 
mechanism  to  be  modeled  In  the  urban  mass  fire  context.  This  process 
depends  upon  the  local  wind  field  and  nay  be  far  less  dependent  upon 
detailed  fuel  parameters  than  spread  In  continuous  fuel  beds.  These 
statements  outline  our  conception  of  the  proper  modeling  approach  to  con¬ 
tiguous  fire  spread  in  the  urban  mass-fire  model. 


In  the  following,  we  briefly  outline  the  very  approximate  calcu¬ 
lation  of  contiguous  spread  which  we  employ  In  the  demonstration  model. 
Contiguous  spread  Is  specialized  to  spread  between  city  blocks  (the  basic 
grid  cells  of  the  model)  across  streets,  which  are  assumed  devoid  of 
fuel.  No  Inter-building  spread  is  explicitly  modeled,  since  city  blocks 
are  the  basic  units.  Spread  between  separated  buildings  within  a  block 
needs  more  Investigation,  along  with  several  other  "sub-grid"  effects, 
which  must  be  parameterized.  Hlnlmal  cell  size  Is  limited  by  the  large 
area  of  the  simulation.  If  Ignited  anywhere,  an  entire  city  block  In  the 
demonstration  model  eventually  burns.  Multiple  ignitions  of  the  same 
block  are  possible  and  result  In  faster  burning.  The  rate  of  burning  Is 
proportional  to  the  number  of  ignitions  of  the  block. 

4.5.1  Summary  of  Contagion  Concepts 

An  Important  observation  of  the  Flambeau  tests  (Reference  9)  was 
that  no  fires  were  Ignited  outside  the  original  fuel  beds  except  within 
the  horizontal  distance  determined  by  the  vertical  projection  of  the  flam¬ 
ing  volume  onto  the  ground.  Hence,  Ignitions  only  occurred  on  the  lee 
side  of  the  fire  In  a  wind.  We  model  spread  across  streets  by  estimating 
the  configuration  of  tilted  flame  columns  In  the  presence  of  a  wind.  Each 
flame  column  has  the  assumed  shape  of  a  cone  with  a  circular  base  whose 
axis  length  and  tilt  are  functions  of  burning  rate  and  wind  velocity. 
Figure  4-6  shows  some  possible  relations  of  a  flame  column,  burning  at  the 
edge  of  a  block,  to  an  open  street  In  the  downwind  direction.  Structures 
are  suggested  by  the  dotted  outlines,  and  the  flame  column  Is  represented 
by  ground-based  triangles.  If  the  flame  column  does  not  overlap  the 
adjacent  block  and  Its  projected  axis  does  not  intersect  buildings  across 
the  street.  Figure  4-6(a),  no  contagion  Is  possible  and  the  adjacent  block 
remains  unignited.  A  finite  Ignition  probability  per  unit  time  Is  cal¬ 
culated  for  the  two  cases  of  Figure  4-6{b)  and  4-6(c).  Figure  4-6(b) 
provides  a  mechanism  for  propagation  between  tall  buildings  separated  by 
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a.  No  contagion.  No  overlap.  No  axis  Intersection. 


b.  Contagion  probability  finite. 

Flame  axis  intersects  buildings  across  street. 


s' 

s' 


c.  Contagion  probability  finite. 

Flame  column  overlaps  buildings  across  street. 

Figure  4-6.  Contagion  across  streets. 


narrow  streets,  and  makes  explicit  use  of  building  height.  The  configura¬ 
tion  of  Figure  4-6(c)  allows  for  propagation  In  accord  with  the  Flambeau 
observations,  and  an  ignition  probability  proportional  to  the  overlap  a  Is 
calculated. 

The  formulas  and  algorithms  used  by  the  demonstration  model  to 
Implement  these  concepts  are  described  in  detail  In  Section  5.  These 
methods  may  require  substantial  revision,  so  our  summary  description  here 
includes  only  general  features  to  Indicate  our  concept  of  the  physical 
factors  Involved. 

Given  the  power  Q  of  a  fire  and  Its  area,  A  *  *rz0,  where 

r0  Is  the  effective  radius  of  the  cone  base,  the  conical  flame  volume  is 

•  • 

computed  from  the  air  supply  required  ior  combustion  Va  "  Q  and  an  as¬ 
sumed  entrainment  velocity  1.0  (m/s)  through  the  lateral  surface  of  the 
cone.  From  the  volume  and  base  ..rea,  the  axial  slant  height,  h,  of  the 
cone  is  computed.  The  angle  6  of  the  cone  axis  Is  determined  from  the 
ratio  of  the  net  vertical  force  on  the  flame,  its  buoyancy,  to  the  net 
horizontal  force.  The  buoyancy  of  the  flame  equals  Vfg(pa  -  p)  where 
Vf  is  flame  volure,  g  is  the  acceleration  of  gravity,  and  pa  -  p  Is 

the  difference  between  ambient  air  density  and  flame  density.  The  hori¬ 
zontal  force  equals  pa  Vav  where  Va  is  the  total  rate  of  air 

entrainment  and  v  Is  the  local  wind  velocity.  These  rough  estimates 
establish  the  flame  geometry  shown  in  Figure  4-6. 

A  problem  arises  because  the  model  carries  only  the  burning 
rate  Q  of  an  entire  city  block.  If  the  base  area  of  the  single  conical 
flame  equals  the  area  of  the  block,  only  rarely  Is  Q  large  enough  to 
result  In  a  flame  tip  located  outside  the  base  area.  Here  again,  we 

encounter  another  aspect  of  the  "sub-grid"  problem  and  the  need  to  esti¬ 
mate  the  number  of  flame  columns  in  area  burning.  To  resolve  this  problem 


(somewhat  Inconsistently  with  procedures  of  the  wind  calculation)  all  the 
burning  is  assumed  to  take  place  within  an  area  such  that  h  >  r0. 
Since,  in  this  case,  the  location  of  the  burning  sub-area  within  the  block 
Is  not  known,  one  employs  a  probability  that  it  Is  located  on  the  downwind 
periphery  adjacent  to  unignited  cells. 

The  geometry  of  the  wind  direction  and  the  flame  column  in  con¬ 
junction  with  the  position  of  burning  cells  in  relation  to  unignited  ones 
forms  the  basis  for  calculating  spread  probability.  These  probabilities 
are,  in  fact,  probabilities  per  unit  time.  Probabilities  are  computed  at 

m 

time  steps  internally  determined  on  the  basis  of  burning  rate  Q  and  its 
derivative  Q.  Furthermore,  these  probabilities  are  accumulated  so  that 
when  the  time  Integral  of  the  probability  per  unit  time  equals  one,  igni¬ 
tion  is  certain. 

A  more  physical  contagion  algorithm,  based  upon  a  ratio  of  heat 
flux  to  heat  of  ignition  as  in  the  Rothermel  model,  might  be  preferable  to 
the  one  we  have  used  in  the  demonstration  model.  Also,  recall  that  short- 
range  branding  may  play  a  major  role  in  contagion.  In  general,  we  feel 
the  demonstration  model  exaggerates  spread  by  contagion.  However,  it  does 
capture  the  dependence  of  this  propagation  mechanism  upon  wind  speed,  wind 
direction,  and  upon  burning  rate. 

4.6  CONCEPT  SUMMARY 

The  preceding  sections  have  emphasized  the  salient  aspects  of 
our  urban  fire  behavior  concept.  Fire  development  takes  place  by  the 
principal  spread  mechanisms  of  contagion  and  branding,  which  are  governed 
by  atmospheric  motion,  a  major  component  of  which  can  be  generated  by  the 
fire  Itself.  The  demonstration  model,  which  has  guided  the  discussion, 
provides  an  explicit  framework  Illustrating  the  implementation  of  these 
concepts.  The  meaning  of  the  demonstration  model,  we  emphasize,  does  not 
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extend  any  farther.  Narrowly  viewed.  It  only  provides  a  suitable  struc¬ 
ture  and  classifies  the  areas  requiring  further  refinement. 

A  measure  of  uncertainty  must  be  accepted  in  a  model  of  urban 
mass  fires.  Not  only  are  there  limits  of  resolution  of  basic  fuel  and 
structural  data,  but  also  there  are  limits  in  the  understanding  of  basic 
mechanisms  all  the  way  through  the  causal  chain  of  events  that  occur. 
This  circumstance  requires  a  judicious  selection  and  balance  of  model  ele¬ 
ments  and  model  structure  in  order  to  mlnimze  the  uncertainties.  The 
concept  presented  Is  our  approach  to  this  problem. 
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SECTION  5 

DESCRIPTION  OF  DEMONSTRATION  URBAN  FIRE  SIMULATION 


The  demonstration  version  of  the  Urban  Fire  Simulation  embodies 
the  range  of  phenomenology  together  with  the  logic  and  data  management 
structures  needed  to  Implement  the  simulation  concept  described  In  Section 
4.  However,  the  phenomenology  models  are  simplified  versions  of  those 
described  in  Section  4.  Development  of  a  demonstration  version  of  the 
simulation  concept  was  undertaken  to  demonstrate  the  feasibility  of  imple¬ 
menting  the  concept,  to  qualitatively  demonstrate  simple  regimes  of 
behavior  for  the  concept,  and  to  identify  areas  in  which  modeling  problems 
might  exist.  This  section  describes  the  structure  and  operation  of  the 
demonstration  simulation.  Simulation  results  are  presented  in  Section  6, 
and  a  discussion  of  modeling  problems  is  included  in  Section  1. 

5.1  DEMONSTRATION  SIMULATION  STRUCTURE  AND  OPERATION 

The  demonstration  simulation  structure  is  designed  to  provide 
flexibility  and  simplicity  of  maintenance.  The  simulation  is  modular, 
event  sequenced,  and  autolnteractlve.  The  simulation  is  asynchronously 
stepped  forward  in  time  through  execution  of  events  sequentially  selected 
from  an  event  list.  Auto  interaction  is  achieved  by  allowing  each  module 
to  insert  events  in  the  event  list. 

The  overall  structure  of  the  Demonstration  Urban  Fire  Simulation 
is  shown  by  the  block  diagram  of  Figure  5-1.  The  meanings  of  the  blocks 
shown  in  Figure  5-1  are  defined  below: 
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ENVTLST  -  the  list  (queue)  of  events  to  be  executed 

UFMGR  -  the  Urban  Fire  Manager  (i.e.,  the  executive  structure) 

TRANS  -  the  event  module  transfer  list 

IGNITE  -  the  functional  module,  which  ignites  unignited  cells 

8CUP0T  -  the  functional  module,  which  updates  burning  cells 

ACUPDT  -  the  functional  module,  which  updates  unignited  cells 
adjacent  to  ignited  cells 

FUPDT  -  the  functional  module,  which  updates  the  status  of 
fires 

BRNDNG  -  the  functional  module,  which  determines  if  branding 
occurs 

WVMAP  -  the  data  processing  module,  which  produces  wind  veloc¬ 
ity  maps 

CSMAP  -  the  data  processing  module,  which  produces  cell  state 
maps 

QMAP  -  the  data  processing  module,  which  produces  heat  pro¬ 
duction  maps. 

The  Dynamic  Storage  and  Allocation  System  (DSAS)  provides  for 
user  interface  to  and  management  of  the  Data  Base.  The  Data  Base  includes 
input  data,  interim  calculations,  and  final  results.  In  Figure  5-1,  sin¬ 
gle  and  double  headed  arrows  are  used  to  indicate  interfaces  between  the 
blocks.  Single  headed  arrows  denote  control  interfaces  and  double  headed 
arrows  denote  data  interfaces. 

The  simulation  structure  allows  a  simple  executive  procedure, 
which  consists  largely  of  two  nested  sequences  of  operations  as  described 
below.  The  Urban  Fire  Manager  accesses  the  Event  List  and  obtains  the 
event  in  the  list  having  the  earliest  execution  time.  With  this  event, 
the  Urban  Fire  Manager  accesses  the  Transfer  List  and  obtains  the  sequence 
of  modules,  which  must  be  processed  to  execute  the  subject  event.  The 
Urban  Fire  Manager  transfers  control  to  the  first  module  in  the  sequence. 
The  module  accesses  the  Data  Base  through  the  Dynamic  Storage  and  Alloca¬ 
tion  System  to  obtain  its  current  input  data,  performs  the  computations 
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required  by  Its  function,  stores  an  event  In  the  Event  List  as  appropri¬ 
ate,  stores  the  results  of  Its  computations  In  the  Data  Base,  and  returns 
control  to  the  Urban  Fire  Manager.  The  Urban  Fire  Manager  continues  to 
transfer  control  to  the  modules  In  sequence  until  all  the  (nodules  In  the 
sequence  have  been  processed.  When  all  of  the  modules  In  the  sequence 
have  been  processed,  event  execution  Is,  by  definition,  completed.  The 
Urban  Fire  Manager  then  accesses  the  Event  List  to  obtain  the  event  having 
the  earliest  execution  time  and  the  event  execution  process  is  repeated 
until  all  events  in  the  Event  List  have  been  executed.  This  executive 
procedure  Is  both  simple  and  flexible.  It  allows  for  altering  events, 
adding  events,  and  adding  modules  without  any  coding  changes. 

The  demonstration  model  computation  process  Is  based  on  a  two- 
dimensional  grid  structure,  which  is  superimposed  on  the  urban  area  of 
interest.  The  grid  partitions  the  urban  area  into  approximately  city 
block  sized  square  cells  (i.e.,  150  x  150  meters).  As  will  be  seen  In  the 
module  descriptions,  the  grid-cell  structure  provides  an  organizational 
mediuii  for  the  input  data,  the  computational  processes,  and  intermediate 
and  final  results.  A  description  of  the  dynamic  data  base,  which  Is  used 
by  each  module,  is  presented  in  the  next  paragraph.  This  Is  followed  in 
subsequent  paragraphs  by  descriptions  of  each  of  the  modules  in  Figure 
5-1. 


5.2  DYNAMIC  DATA  BASE  DESCRIPTION 

The  dynamic  data  base  Is  keyed  to  the  grid-cell  structure  and 
partitioned  into  data  sets  consistent  with  the  types  of  data  required  or 
produced  by  the  modules.  These  data  sets  are  described  below. 

Grid  Data:  The  grid  data  are  used  to  define  the  grid  dimensions 
for  a  specific  simulation  execution  and  can  be  changed  between  simulation 
executions.  The  grid  data  are  as  follows: 
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1.  The  cell  dimension  (assumes  square  cells) 

2.  The  location  of  the  rectangular  coordinate  system  origin 

3.  The  maximun  number  of  cells  In  each  of  the  two  dimensions. 

Cell  State  Data:  Each  cell  in  the  grid  may,  at  any  time,  be  in 
one  of  the  following  five  states: 

1.  Unignlted 

2.  Unignited,  but  adjacent  to  a  burning  cell  (Adjacent) 

3.  Ignited  and  on  a  fire  periphery  (Peripheral) 

4.  Ignited  and  In  a  fire  interior  (Interior) 

5.  Burned  out. 

The  cell  state  data  are  an  array  containing  the  state  of  each  cell  in  the 
grid.  This  array  is  updated  each  time  a  cell  changes  state. 

Fuel,  Data:  The  fuel  within  each  cell  is  assumed  to  be  uniform 
and  the  fuel  type  within  each  cell  is  identified  in  terms  of  an  occupancy 
class.  An  arbitrary  number  of  occupancy  classes  can  be  defined,  ranging 
from  heavy  industrial  through  single  family  residential  to  parks  and  unim¬ 
proved  areas.  The  fuel  data  specified  for  each  occupancy  class  are: 

1.  Fuel  loading(  kg/cell) 

2.  Fuel  heat  content  (J/kg) 

3.  Fuel  height  (m) 

4.  No  wind  combustion  rate  (m/s) 

5.  Scale  factor  for  time  of  maximun  heat  production. 

Event  List:  The  event  list  is  a  queue  of  all  the  events  that 
have  been  scheduled  for  execution  and  Is  updated  each  time  an  event  is 
taken  from  or  added  to  the  queue. 
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Burning  Cell  Data:  The  burning  cell  data  consist  of  a  sequence 
of  three  data  sets  for  each  ignited  ceil.  Each  cell  is  assigned  a  unique 
Identification  number  (ID)  based  upon  its  row-column  location  in  the 
grid.  The  cell  ID  is  used  to  access  the  cell  state  data,  and  the  cell 
state  together  with  cell  ID  are  used  to  access  data  for  both  ignited  and 
unignited  cells.  The  following  burning  cell  data  are  maintained  for  each 


ignited  cell : 

1. 

Entry  ID  (1) 

2. 

Occupancy  class  (fuel  type) 

3. 

Cell  topographic  height  (m) 

4. 

Time  of  cell  Ignition  (s) 

5. 

Number  of  times  cell  has  been  ignited 

6. 

Cell  ID 

7. 

Entry  ID  (2  or  3) 

8. 

Current  cell  heat  production  rate  (J/s) 

9. 

Maximum  cell  heat  production  rate  (0/s) 

10. 

Current  cell  burn  time  (s) 

11. 

Time  of  cell  maximum  rate  of  heat  production  (s) 

12. 

Time  at  which  these  data  are  valid. 

Data  elements  2  through  6  do  not  change,  and  entry  ID  =  1  denotes  this 
subset  of  a  burning  cell  data  set.  Data  elements  8  through  12  may  change 
at  each  calculation  time,  and  this  is  denoted  by  entry  ID  *  2  or  3.  Sub¬ 
sets  of  calculated  burning  cell  data  are  maintained  at  two  times,  which 
always  subtend  the  simulation  time.  Data  required  at  any  other  times  are 
obtained  by  Interpolating  between  the  data  stored  at  these  two  times. 
Entry  ID  *  2  denotes  the  earlier,  and  entry  ID  3  denotes  the  later  of 
these  two  data  subsets. 
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Adjacent  Cell  Data:  Adjacent  cell  data  consist  of  a  sequence  of 
data  sets  for  each  unignited  cell  that  Is  adjacent  to  a  burning  cell.  The 
following  Information  Is  maintained  for  each  adjacent  cell: 

1.  Entry  ID  (1) 

2.  Occupancy  class  (fuel  type) 

3.  Cell  topographic  height  (m) 

4.  Time  of  last  cell  update 

5.  Cumulative  probability  of  cell  Ignition 

6.  Cell  ID 

Data  for  entry  ID  *  1  are  included  in  this  case  because  the  data  for  an 

adjacent  cell  will  be  used  to  form  the  first  subset  of  data  for  a  burning 

cell.  When  the  cell  is  ignited,  the  entry  ID  -  1  subset  will  be  updated, 
and  subsets  for  entry  levels  2  and  3  will  be  generated. 

Unignited  Cell  Data:  Unignited  cell  data  consist  of  a  sequence 
of  data  sets  for  each  unignited  cell  In  the  grid  that  is  not  adjacent  to  a 
burning  cell.  The  following  quantities  are  maintained  for  each  unignited 
cel  I : 

1.  Entry  ID  (1) 

2.  Occupancy  class  (fuel  type) 

3.  Cell  topographic  height  (m). 

Data  for  entry  ID  =  1  are  included  in  this  case  because  the  data  for  an 

unignited  cell  will  be  used  to  form  the  entry  ID  *  1  subset  of  data  for  an 

adjacent  cell  when  the  subject  cell  becomes  adjacent  to  a  burning  cell. 

Fire  Data:  A  fire  is  defined  to  consist  of  a  set  of  contiguous 
burning  cells.  Depending  upon  its  geometry,  a  fire  may  be  partitioned 
into  subfires  that  are  reasonably  symmetrical.  This  is  done  because  some 


of  the  phenomenology  models  Include  assumptions  regarding  geometric  sym¬ 
metry  of  the  fires  they  describe.  The  following  data  are  maintained  for 
each  defined  fire: 

1.  Fire  ID 

2.  Location  of  the  effective  center  of  the  fire  (x,y  coordin¬ 
ates  -  m) 

3.  Effective  radius  of  the  fire  (m) 

4.  Fire  heat  production  rate  (J/s) 

•  5.  Fire  area  (m2) 

6.  Volume  rate  of  surface  air  inflow  to  fire  (m3/s) 

7.  Time  at  which  these  data  were  calculated  (s) 

8.  Number  of  cells  in  fire  at  the  time  given  in  (7) 

9.  List  of  cells  currently  in  fire. 

The  use  of  an  effective  center  and  radius  for  each  fire  recognizes  that 
the  rate  of  heat  production  over  a  fire  will  generally  not  be  uniform. 
Each  time  a  cell  is  added  to  (ignited)  or  deleted  from  (burned  out)  the 
fire  data,  element  9  is  updated.  The  other  eight  data  elements  may  be 
updated  much  less  frequently  for  reasons  of  computational  efficiency. 

The  data  sets  described  above  con  titute  most  of  the  dynamic 
data  base,  which  is  accessible  to  and  updated  by  the  simulation  models. 
The  remaining  data  consist  of  flags,  constants,  and  other  numerical  quan¬ 
tities  that  are  not  generally  part  of  the  dynamic  data  base.  The  follow¬ 
ing  paragraphs  describe  the  logic  flow  and  data  interfaces  within  the 
modules  shown  in  Figure  5-1. 

5.3  IGNITE  MODULE 

The  IGNITE  module  initializes  burning  in  cells  that  have  been 
determined  to  have  been  ignited.  As  shown  in  Figure  5-2,  the  event  data 
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CANDIDATE  CELL  ID 
NO.  OF  IGNITIONS 
SIMTIM 


IS  CELL  IGNITED  OR  BURNED  OUT? 


CALL  8CINIT 


INITIALIZE  BURNING  CELL 


Figure  5-2.  Logic  flow  of  I6NITE  nodule. 


passed  to  IGNITE  are  the  Identity  of  the  candidate  cell,  the  number  of 
Ignitions  and  the  current  simulation  time  (SIMTIM).  The  candidate  cell 
status  is  tested  to  determine  If  the  cell  can  be  Ignited.  If  the  cell  Is 
already  ignited  or  burned  out,  It  cannot  be  Ignited  and  IGNITE  returns 
control  to  UFMGR.  If  the  candidate  cell  Is  unignited,  the  Burning  Cell 
Initialization  Subroutine  is  called  and  then  control  is  returned  to  UFMGR. 

The  Burning  Cell  Initialization  Subroutine  (BCINIT)  performs  the 
tasks  associated  with  changing  the  status  of  a  cell  from  unignited  to 
ignited  as  shown  In  Figure  5-3.  BCINIT  Ignites  the  candidate  cell  at 
SIMTIM  and  changes  the  cell  state  from  unignited  to  Ignited.  While  the 
cell  was  unignited,  the  entry  ID  3  2  and  3  cell  data  subsets  did  not 
exist.  The  entry  ID  3  1  cell  data  subset  is  initialized  at  SIMTIM  with 
all  zero  values.  The  wind  velocity  in  the  just  ignited  cell  is  calculated 
at  SIMTIM  +  360  seconds  by  the  Wind  Velocity  Subroutine  (WINDV).  The 
entry  ID  3  2  data  subset  is  initialized  at  SIMTIM  +  360  seconds  with 
values  calculated  for  that  time,  and  a  burning  cell  update  event,  is  placed 
in  the  event  list  for  execution  time  at  SIMTIM  +  360  seconds.  Any  sub¬ 
routine  requiring  data  on  this  cell  between  the  current  time  and  current 
time  plus  360  seconds  will  obtain  it  by  interpolating  between  these  two 
data  subsets.  At  current  time  plus  360  seconds,  execution  of  the  burning 
cell  update  event  will  reset  data  subset  ID  =  2  to  ID  *  1,  generate  a  new 
data  subset  ID  =  2  at  a  future  time  dependent  upon  cell  heat  production 
rate,  and  place  an  update  event  in  the  event  list  for  execution  at  the 
same  future  time.  In  this  way,  two  data  subsets  that  always  subtend 
SIMTIM  are  maintained  for  each  burning  cell  . 

The  cell  state  and  the  fire  to  which  it  belongs  (if  burning)  are 
identified  for  each  of  the  eight  cells  adjacent  to  the  just  Ignited  cell. 
Also,  the  oldest  fire  represented  among  the  eight  adjacent  cells  is 
identified.  Fires  are  assigned  unique  ID  numbers  in  sequential  order,  so 
that  older  fires  have  lower  ID  numbers.  When  fires  are  merged,  younger 
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CANDIDATE  CELL  ID 
NO.  OF  IGNITIONS 
SIMTIM 


Figure  5-3.  Logic  flow  of  burning  cell  Initialization  module. 
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ADO  IGNITED  CELL  TO  OLDEST  FIRE  REPRESENTED 
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Figure  5-3.  Logic  flow  of  burning  cell  Initialization  module 
(continued) . 
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IDENTIFY  8  CELLS  ADJACENT  TO  THIS 
ADJACENT  CELL  AND  IDENTIFY  THE 
STATE  OF  EACH  CELL 


ARE  ALL  3  OF  THESE  CELLS  IGNITED? 


CHANGE  THE  STATE  OF  THIS  CELL  FROM 
PERIPHERAL  TO  INTERIOR 


HAVE  ALL  8  CELLS  ADJACENT  TO  IGNITED 
CELL  BEEN  EXAMINED? 


<? 

STORE  CELL  DATA  AND  FIRE  DATA 


C 


RETURN 


Logic  flow  of  burning  cell  initialization  module 
(concluded). 


fires  are  always  merged  Into  older  fires  to  maintain  continuity  of  (older) 
fire  Identities. 


If  at  least  one  cell  adjacent  to  the  just  ignited  cell  is  burn- 
ing,  BCINIT  transfers  to  label  (A),  otherwise  transfer  is  to  le.be  1  (B). 
At  label  (A),  the  fire  Identities  of  all  of  the  ignited  adjacent  cells  are 
examined.  If  the  fire  10  of  an  ignited  cell  is  different  from  the  ID  of 
the  oldest  fire  represented  among  the  adjacent  cells,  transfer  is  to  label 
(D).  At  label  (D),  the  fire  that  the  subject  burning  cell,  adjacent  to 
the  just  ignited  cell  belongs  to,  is  merged  with  the  oldest  fire  repre¬ 
sented  among  the  eight  adjacent  cells.  Merging  the  two  fires  involves 
combining  the  burning  and  adjacent  unignited  cell  lists  of  the  younger 
fire  with  the  burning  and  adjacent  unignited  cell  lists  of  the  older  fire, 
updating  the  fire  geometry  partitioning,  and  updating  the  fire  physical 
data.  After  all  of  the  required  mergers  have  been  performed,  the  just 
ignited  cell  is  added  to  the  oldest  fire  represented  among  its  eight 
adjacent  cells  and  a  transfer  is  made  to  label  (E). 

At  label  (B),  none  of  the  eight  cells  adjacent  to  the  just 
ignited  cell  were  found  to  be  ignited.  This  is  the  definition  of  an  iso¬ 
lated  fire.  Thus,  a  new  fire  is  created  for  the  just  ignited  cell.  This 
involves  assigning  to  this  cell  the  next  fire  ID  number,  creating  burning 
and  adjacent  unignited  cell  lists,  and  initializing  the  fire  physical 
data. 


The  Wind  Velocity  Subroutine  (WINDV)  calculates  the  wind  veloc¬ 
ity  at  a  specified  point  due  to  all  fires  in  the  grid  area  and  due  to  the 
ambient  wind.  The  wind  velocity  calculation  assumes  a  flux  sink  at  the 
center  of  each  circularly  symmetrical  subfire.  As  mentioned  before,  fires 
consist  of  sets  of  contiguous  burning  cells,  which  are  partitioned  into 
subfires  that  are  approximately  symmetrical.  Thus,  each  fire,  even  if  not 
partitioned,  includes  one  subfire  (i.e.,  the  whole).  As  shown  in  the 
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logic  flowchart  of  Figure  5-4,  WINDV  loops  over  all  fires  and  subfires  In 
the  grid  and  calculates  the  contribution  of  each  subfire  to  t.he  fire 
Induced  wind  velocity  at  a  specified  field  point.  Three  cases  are  consid¬ 
ered  In  this  calculation:  (1)  the  field  point  Is  located  at  the  center  of 
the  subfire  under  consideration;  (2)  the  field  point  Is  Inside  the  subfire 
under  consideration;  and  (3)  the  field  point  Is  outside  of  the  subfire 
under  consideration.  Upon  completion  of  the  fire  Induced  wind  calcula¬ 
tions,  the  ambient  wind  velocity  is  added  to  the  fire  Induced  wind  veloc¬ 
ity  to  define  the  total  wind  velocity  at  the  specified  point. 

Calculation  of  the  fire  induced  wind  velocity  is  based  on  the 
following  equations. 

w  ■  X  w  *  A  w  <5‘1> 

where 

Vp(rp)  a  the  fire  induced  wind  velocity  at  the  grid  point  rp  , 

f*l,  N  are  the  subfires  that  do  not  contain  grid  point  rp  , 

f ' =1,  M  are  the  subfires  that  do  contain  grid  point  r  . 

The  contribution  to  the  fire  induced  wind  at  point  ?p  due  to  subflro  f, 
which  does  not  contain  grid  point  rp,  is  given  by 
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Figure  5-4.  Logic  flow  of  wind  velocity  subroutine 
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Figure  5-4.  Logic  flow  of  wind  velocity  subroutine  (concluded) 
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The  contribution  to  the  fire  induced  wind  at  point  rp  due  to  subfire  f, 
which  does  contain  grid  point  rp,  is  given  by 


W 


V„ 


2  ir  r  2 , 


+  g(w)h  (lilt)  q-iMf9 

M  i  l  l 


(5-3) 


where 


Vfi  =  volume  rate  ^f  air  inflow  to  subfire  f ' , 
h'  =  ?cf'-  ?p  » 

rcf.  -  location  of  the  effective  center  of  subfire  f ' , 

• 

Qfi  =  total  energy  release  rate  of  subfire  f', 
r^i  =  effective  radius  of  subfire  f ' . 

The  function  g(y)  is  given  by 

g(y)  =  0.5298’f1 ' 5863/(l  .3285  +  y1’5863)  (5-4) 


where 


Y  ■  Zc/|Sf|  if  |Sf|  >  rf  or 


(5-5) 
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(5-6) 


Zc/rf  if  ) af )  <  rf  , 


Zc  »  0.23376  qj/4  t3/4 


Q^x  3  total  energy  release  rate  of  the  largest  fire  in 
the  grid 

t  =  simulation  time. 

The  function  h(|afj/rf)  is  given  by 


2  cos" 


(5-7) 


5.4  BURNING  CELL  UPDATE  (BCUPDT)  HODULE 

The  Burning  Cell  Update  Module  maintains  dynamic  cell  data  that 
provide  the  characteristics  of  a  burning  cell.  Two  dynamic  cell  data 
lists  are  maintained  that  always  subtend  SIMTIM.  Characteristics  at  any 
SIMTIM  are  determined  by  interpolating  with  respect  to  time  between  these 
two  data  lists. 


The  energy  production  rate  of  a  burning  cell  is  defined  by 


(J/s) 


(5-8) 
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where 


i 


qc  a  energy  production  rate  of  cell  c  (0/s) 

Qc  -  total  energy  content  of  cell  (J) 
t^  -  cell  time  of  maximum  energy  production  (s) 
t  *  cell  burn  time  (s). 

The  total  cell  energy  content  is  given  by 

Qc  -  HfLfAc  ,  (j)  (5-9) 

where 

Hf  *  fuel  heat  content  (J/kg) 

Lp  3  fuel  loading  (kg/m2) 

Ac  -  cell  area  (m2). 

The  cell  heat  production  rate  given  by  Equation  5-8  rises  Lo  a  peak  value 
at  t  *  rmx  and  then  decays  as  shown  In  Figure  5-5. 

As  shown  In  the  logic  flowchart  of  Figure  5-6,  the  information 
transferred  to  BCUPDT  by  the  update  event  is  the  ID  of  the  burning  cell 
and  the  current  SIMTIM. 

BCUPDT  Identifies  the  fire  the  burning  cell  belongs  to  by  refer¬ 
ring  to  the  cell  state  map.  In  this  map,  a  burning  cell  is  Identified  by 
the  ID  of  the  fire  to  which  It  belongs. 

The  fraction  of  the  cell  that  has  been  burned  at  current  SIMTIM 
Is  calculated.  As  shown  In  Equation  5-8  and  Figure  5-5,  the  cell  combus¬ 
tion  function  Is  characterized  by  a  maximum  heat  production  rate 
(qc)mx,  ^ch  occurs  at  time  ^x.  Cell  burnout  Is  defined  to  occur 
at  t  ■  4  trax,  at  which  time  the  cell  is  approximately  90  percent 
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Figure  5-6.  Logic  flow  of  burning  cell  update  module  (continued) 
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Figure  5-6.  Logic  flow  of  burning  cell  update  module  (concluded), 
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consumed,  according  to  Equation  5-8.  Thus,  a  measure  of  the  fraction  of 
the  cell  that  has  been  burned  is  given  by  t/4^x. 

The  current  energy  production  rate  of  the  cell  Is  calculated 
from  the  dynamic  cell  data.  This  cell  production  rate  Is  used  to  update 
the  energy  production  rate  of  the  fire  to  which  the  cell  belongs.  In  the 
case  of  fires  (as  opposed  to  cells),  only  one  set  of  data  Is  maintained  to 
describe  the  fire  status.  In  between  fire  update  events,  the  fire  state 
is  assumed  not  to  change,  except  for  fire  heat  production  rate,  which  Is 
•adjusted  each  time  one  of  the  fire's  constituent  cells  is  updated. 

A  burning  cell  update  event  Is  always  scheduled  to  occur  when 
SIMTIM  is  equal  to  the  time  of  the  most  recent  dynamic  set  of  cell  data. 
A  significant  function  of  the  burning  cell  update  is  to  determine  the  next 
update  time,  project  dynamic  cell  data  to  that  time,  and  arrange  for  a 
cell  update  at  that  time. 

The  time  to  be  associated  with  the  new  data  set  and  with  the 
next  update  of  the  burning  cell  is  based  on  the  current  rate  of  change  of 
cell  heat  production  rate.  The  time  increment  to  the  next  update  is  given 

by 

dT  *  0,5  Tmx  ^  “  q/qo*  <s)»  (5-10) 

where  t^,  ■  d/(t*0).  This  provides  for  updates  at  approximately  2000 
second  intervals  when  t  >  tfflx  and  provides  more  frequent  updates  when  t 
<  tmx.  The  rate  of  change  of  cell  heat  production  rate  with  respect  to 
time  (qc)  Is  from  Equation  5-8 

«c  *  <V4>d  -  e't/,mx  (5-11) 
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A  call  is  made  to  WINDV  to  calculate  the  wind  velocity  in  the 
burning  cell  at  t'  s  t+dt,  the  time  of  the  next  cell  update. 

At  label  (A)  of  Figure  5-6,  the  time  of  maximum  cell  heat 
production  is  calculated  from 

V  k  Lc/<N19Ms  +  Sw»  <s>>  <5-12> 

where 

k  *  fuel  constant 
Lc  a  cell  dimension  (m) 

N 3  number  of  ignitions  in  cell 
hf  3  fuel  height  (m) 

a  =  zero-wind  spread  rate  coefficient  (m/s) 

Sw  =  total  wind  speed  in  ce1!  (m/s). 

The  effects  of  increasing  wind  speed  upon  cell  combustion  are,  thus,  to 
increase  the  maximum  cell  heat  production  rate  and  to  cause  its  occurrence 
at  an  earlier  time,  Maximun  cell  heat  production  rate  occurs  at  t  = 
xmx,  which  from  Equation  5-8  yields 

*  °-368  VV<  •  <5-131 

where  Qc  is  defined  by  Equation  5-9  and  tmx  by  Equation  5-12. 

The  time  of  maximun  heat  production  Is  a  function  of  total  wind 
speed  in  the  cell,  which  Is  a  function  of  the  geometry  and  Intensity  of 
all  the  fires  In  the  grid,  together  with  the  ambient  wind.  The  wind  In 
the  cell  may,  thus,  fluctuate  and  cause  both  the  time  and  level  of  maximum 
heat  production  rate  to  fluctuate.  These  fluctuations  are  acceptable  and 
represent  an  important  part  of  the  fire  wind  interaction. 


A  change  In  wind  speed  affects  Tfflx,  (4c)mx»  <lc  »nd  burnout  time, 
together  with  other  cell  and  fire  variables.  Given  a  change  In  we 
must,  to  be  physically  consistent,  require  the  cumulative  heat  production 
to  be  unchanged.  The  cumulative  heat  production  rate  (Q  )  Is,  from  Equa- 

V 

tion  5-8,  given  by 

0.(t)  ■  /T  5  dt  -  Q[1  -  (1  +  -S-)  e-T/^]  .  (5-14) 

c  o  mx 

which  we  write  In  the  form 


Qc(f)  *  Q[1  -  (1  +  f)  e'f]  ,  (5-15) 

where  f  *  x/xm  is  a  measure  of  the  fraction  of  the  cell  that  has  been 
burned.  The  cumulative  cell  heat  production  is  held  constant  under  vary¬ 
ing  tmx  by  also  varying  t  so  that  f  remains  constant.  The  fraction 
burned  at  the  next  cell  update  time  (t')  with  the  current  value  of  t^x 
is 


f«  -  T',xm  ■  (5-16) 

The  fraction  burned  at  the  next  cell  update  time  with  the  value  of  T*mx 
at  that  time  Is 

f?-  s  T'(^x/Tmx)  ‘  (5‘17) 

Calculate  the  new  value  of  cell  time  at  the  next  update  (t")  of  this  cell 
as 


Xmx7xmx^ 


(5-18) 
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by  requiring  that  fj  ■  f2.  The  manipulation  at  cell  time  Is  possible 
without  affecting  simulation  time  because  the  cell  time  Is  coupled  to 
SIMTIM  only  by 

SIMTIM  -  t1  <  t  <  4^x  -  SIMTIM  .  (5-19) 

Within  the  limitations  Imposed  by  Equation  5-19,  t  may  advance  or  regress, 

as  required  by  Equation  5-18,  without  affecting  SIMTIM. 

Given  the  adjusted  time  of  the  next  cell  update  (t"),  the  test 

t"  >  4Tmx  Is  performed.  If  the  up  i  time  Is  after  cell  burnout,  the 

update  time  Is  adjusted  at  label  (B)  to  coincide  with  cell  burnout. 

If  t"  <  4tmx,  transfer  Is  made  to  label  (C).  If  t"  a  4^, 
transfer  Is  to  label  (D),  otherwise  transfer  is  to  label  (J). 

At  Label  (D),  a  burned  out  cell  has  been  identified.  The  state 
of  the  subject  cell  is  changed  from  burning  to  burned  out,  and  adjustments 
are  made  to  the  related  fire  data. 

At  label  (E)  the  subject  cell  was  the  only  burning  cell  In  the 
related  fire;  thus,  the  fire  Is  deleted  by  deleting  its  burning  cell  and 
adjacent  cell  lists.  Processing  of  this  cell  is  then  complete  and  BCUPDT 
returns  control  to  UFMGR. 

At  label  (F)  the  subject  cell  was  not  the  only  burning  cell  In 
the  related  fire;  thus,  the  subject  cell  is  deleted  from  the  burning  cell 
list.  The  8  cells  adjacent  to  the  subject  cell  are  Identified  to  deter¬ 
mine  If  their  state  has  been  changed  by  the  state  change  of  the  subject 
cell.  All  burning  cells  are  either  interior  cells  or  peripheral  cells. 


An  Interior  cell  Is  one  that  has  8  burning  adjacent  cells  and  a  peripheral 
cell  Is  one  that  has  at  least  one  nonburning  (un Ignited  or  burned  out) 
adjacent  cell.  Thus,  the  only  state  change  possible  for  a  cell  adjacent 
to  a  newly  burned  out  cell  Is  a  change  from  an  Interior  cell  to  a  periph¬ 
eral  cell.  In  labels  (G)  through  (I)  all  adjacent  cells  of  the  subject 
cell  that  were  Interior  cells  are  changed  to  peripheral  cells. 

At  label  (J)  all  fire  data  manipulations  have  been  completed. 
The  burning  cell  heat  production  rate  at  the  adjusted  time  of  next  cell 
update  Is  calculated.  A  burning  cell  update  event  for  the  subject  cell  Is 
placed  in  the  event  list  for  execution  at  the  SIMTIM  corresponding  to  the 
adjusted  cell  update  time.  At  label  (K)  the  new  set  of  dynamic  cell  data 
are  stored  in  the  dynamic  data  base  In  place  of  the  older  of  the  two 
existing  sets  of  dynamic  cell  data  for  this  cell.  The  data  stored  are 
Entry  ID,  Heat  Production  Rate,  Maximum  Heat  Production  Rate,  Cell  Tlmetof 
These  Data,  Cell  Time  of  Maximum  Heat  Production  Rate,  and  SIMTIM  of  these 
data.  After  storing  these  data  in  the  dynamic  data  base,  BCUPDT  returns 
control  to  UFMGR. 

5.5  ADJACENT  CELL  UPDATE  (ACUPDT)  MODULE 

As  shown  in  Figure  5-7,  the  information  passed  to  the  ACUPDT 
module  by  the  update  event  Is  the  current  SIMTIM  and  the  ID  of  the  cell  to 
be  updated.  The  cell  to  be  updated  is  called  an  adjacent  cell  because  It 
is  adjacent  to  a  fire.  We  will  also  be  concerned  with  th2  eight  cells 
adjacent  to  the  cell  to  be  updated.  To  prevent  confusion,  the  cell  to  be 
updated  is  referred  to  here  and  in  Figure  5-7  as  the  subject  cell. 

An  unignited  cell  adjacent  to  a  fire  may  experience  a  state 
change  by  becoming  Ignited.  Ignition  may  result  from  branding  or  from 
contagion.  Branding  is  defined  as  Ignition  due  to  the  physical  transport 
of  burning  material  and  will  be  discussed  in  a  later  section.  Contagion 
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Figure  5-7.  Adjacent  cell  update  wodule  (continued) 
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Figure  5-7.  Adjacent  cell  update  module  (contfnued) 
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Is  defined  as  Ignition  due  to  the  transport  of  thermal  energy  by  conduc¬ 
tion,  convection,  or  radiation.  The  purpose  of  ACUPDT  Is  to  determine  If 
and  when  the  subject  will  be  Ignited  by  contagion. 

The  subject  cell  ID  provides  access  to  the  appropriate  cell 
data,  which  contains  the  cell  fuel  type  ID.  The  fuel  type  ID  provides 
access  to  the  fuel  characteristics  associated  with  that  fuel  type.  The 
eight  cells  adjacent  to  the  subject  cell  are  Identified. 

Label  (A)  Is  the  beginning  of  a  loop  over  the  eight  cells  adjac¬ 
ent  to  the  subject  cell  In  which  only  burning  cells  will  be  examined. 
Observe  that  since  the  subject  cell  is  adjacent  to  a  fire,  at  least  one  of 
its  adjacent  cells  must  be  Ignited.  Each  adjacent  cell  Is  tested  (via  the 
cell  state  matrix)  to  determine  If  it  is  burning.  If  it  Is  not,  the  next 
cell  Is  examined  via  label  (L).  If  the  cell  Is  burning,  transfer  is  to 
label  (B). 


At  label  (B)  the  unit  vector  is  calculated  in  the  direction  from 
the  cent'-r  of  the  burning  adjacent  cell  to  the  center  of  the  subject 
cell.  At  label  (C)  a  call  is  made  to  WINDV  to  calculate  the  wind  velocity 
in  the  burning  adjacent  cell.  Examination  of  the  scalar  product  of  the 
un  vector  and  wind  velocity  reveals  if  a  wind  component  exists  In  the 
direction  from  the  burning  adjacent  cell  to  the  subject  cell.  If  there  Is 
no  wind  component  In  the  direction  of  the  subject  cell,  the  next  adjacent 
cell  Is  examined  via  label  (L).  If  there  is  a  wind  component  in  the 
direction  of  the  subject  cell,  transfer  Is  to  label  (0). 

At  label  (D)  the  dynamic  cell  data  are  accessed  for  the  burning 
cell  adjacent  to,  and  upwind  of,  the  subject  cell.  If  this  cell  has  just 
been  Ignited  (at  the  current  SIHTIM),  its  heat  production  rate  Is  still 
zero  and  the  next  adjacent  cell  is  examined  via  label  (L).  If  the  cell 
was  Ignited  prior  to  the  current  SIMTIM,  transfer  Is  to  label  (E). 


At  label  (E)  all  of  the  preliminary  conditions  have  been  satis¬ 
fied  for  contagion  to  occur.  An  Ignited  cell  with  nonzero  heat  production 
has  been  Identified  adjacent  to,  and  upwind  of,  the  subject  cell.  A  test 
Is  performed  to  determine  If  a  plume  can  be  defined  for  the  entire  burning 
cell.  The  plume  test  assumes  a  uniformly  burning  circular  area  of  radius 
r0,  a  constant  entrainment  velocity  ve,  and  a  conical  plume  of  base 
radius  r0  and  height  h.  The  volume  rate  of  air  entrained  by  the  plume 
(V6)  is 


We  ’ 


where  the  cone  area  Ap  is  given  by 


Ap  =  *ro(ro2  +  *2)1/2  ‘ 


(5-20) 


(5-21) 


Combine  Equations  5-20  and  5-21  and  solve  for  the  plume  height  h  to  obtain 


°r°[Mr)  _11 


(5-22) 


The  total  heat  production  Qc  of  the  cell  Is 


m. 


5c  -  (*£)  • 


(5-23) 


where 


q0  a  heat  content  of  fuel  (0/kg) 
pa  ■  air  density  (kg/m3) 
k  ■  mass  of  air  per  mass  of  fuel. 


The  average  heat  production  rate  over  the  cell  (c[)  is 

q  -  Qc/Ag  (5-24) 

where  a  base  area  of  conical  plume.  From  Equations  5-23  and  5-24, 
Equation  5-22  can  be  written  In  the  form 


h 


(5-25) 


where 


q  ■  average  heat  production  rate  per  unit  cell  area. 

The  condition  for  a  cell  to  have  a  definable  plume  Is  that  the  plume 
height  (h)  as  defined  by  Equation  5-25  be  real  valued.  For  h  to  be 
real,  we  must  have 


(5-26) 


If  Equation  5-26  is  satisfied,  transfer  is  to  label  (F),  else  transfer  is 
to  label  (G). 

At  label  (F)  the  base  area,  base  radius,  and  height  of  the 
conical  plume  are  defined  from  Equation  5-25  and  the  cell  area,  and 
transfer  is  to  label  (H). 

At  label  (G)  a  plume  cannot  be  defined  for  the  entire  cell.  It 
is  assumed  that,  since  the  cell  is  burning,  a  plume  exists  and  can  be 
defined  for  some  fraction  of  the  cell.  The  appropriate  cell  fraction  is 
defined  by  requiring  that  the  minimum  plume  height  h^  be  equal  to  the 
radius  of  the  burning  cell  fraction  responsible  for  its  generation.  In 
terms  of  a  fraction  f  of  cell  area  ^  =  L2C,  Equation  5-25  can  be 
written 


(5-27) 


where  the  effective  cell  fraction  in  plume  generation  is  obtained  by 
partitioning  cf  and  assuming  that  q  is  uniform  over  the  cell.  The 
requirement  for  h^  s  rmn  yields,  from  Equation  5-27,  a  definition  of 
the  cell  fraction  f,  which  is 


f  -  (k^)/(  2  ve) 


(5-28) 


The  plume  base  area,  base  radius  and  height  are  defined  as 


(5-30) 


ro  “  (Ap/w^/2 

h  3  rp  ,  (5-31) 

respectively,  and  transfer  is  to  label  (H) . 

At  label  (H)  a  plume  has  been  defined,  either  for  the  entire 
cell  or  for  a  fraction  of  the  cell,  and  we  desire  to  calculate  the  plume 
tilt  resulting  from  wind  velocity  in  the  burning  cell  adjacent  to  the 
subject  cell.  The  horizontal  ambient  wind  force  Fx  is  given  by 


Vave 


(5-32) 


and  the  vertical  buoyancy  force  F2  is  given  by 

FZ=  V(pa-  p1}  »  (5-33) 

where 

VG  3  volume  rate  of  air  entrained  (m3/s) 
ve  3  entrainment  velocity  (m/s) 
pa  3  ambient  air  density  (kg/m3) 

Vp  3  plume  volume  (m3) 
g  3  gravitational  acceleration  (m/s2) 

Pi  3  air  density  inside  plume  (kg/m3) 

The  inclination  e  of  the  force  vector  measured  from  the  horizontal  is 
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(5-34) 


tan  a  *  Fz/Fx  =  (0.8  Vpg)/(v’Gve) 

where 


(Pa  “  Pi)/pa  *  °-8 

The  velocity  (Vp)  of  a  conical  plume  tilted  (*/2  -  a)  radians  from  the 
vertical  is 

vp  =  (1/3)hAb  sin  0  (5-35) 


where 

h  a  cone  height  (m) 

A[j  -  cone  base  area  (m) . 


Combine  Equations  5-34  and  5-35  to  obtain 

1  tan  0  0.8 

— t —  «*  -  -  ziz  .  (5-36) 

cos  e  sin  e  3  u  v 

6  e 

Substitution  of  Equations  5-23,  5-24,  apd  5-25  into  5-36  and  simplifying 
yields 


1  .  0.8 
CS'5'15  ”  ~2~ 


(5-37) 


The  angle  e  given  by  Equation  5-37  is  defined  to  be  the  plume  tilt  angle. 
It  can  be  seen  that  Equation  5-37  may  be  imaginary  and  that  the  condition 
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for  Equation  5-37  to  be  real  valued  Is  the  same  as  the  condition  for 

Equation  5-25  to  be  real  valued.  Thus,  the  plume  definition  tests,  which 

define  h  to  be  real  valued,  also  assure  that  0  Is  real  valued. 

The  cells  In  the  grid  are  considered  to  be  separated  by  streets, 
which  are  lines  of  specified  width  that  contain  no  fuel.  The  ACUPDT 
module  considers  two  methods  by  which  heat  from  a  burning  cell  may  Ignite 
an  adjacent  unignited  cell  across  a  street.  The  method  considered  depends 
upon  whether  the  plume  centerline  lies  above  (label  I)  or  below  (label  J) 
the  top  of  the  fuel  In  the  un Ignited  adjacent  cell. 

At  label  (I)  the  plume  is  tilted  sufficiently  far  from  the 

vertical  that  the  plume  centerline  lies  below  the  top  of  the  fuel  In  the 
adjacent  cell.  In  this  case,  it  is  assumed  that  the  plume  Is  blowing 
directly  on  the  unignited  fuel  in  the  adjacent  cell,  and  the  probability 
of  Ignition  of  the  subject  cell  is  calculated  as 


P 


ig 


min(T/r 


mx’ 


(5-38) 


where 


t/tmx  =*  ratio  of  cell  burntime  to  time  of  maximum  heat 
production 

A 

u  a  unit  vector  from  the  burning  cell  center  to  the 
unignited  adjacent  cell  center 

vq  ■  wind  velocity  In  the  burning  cell 
e  a  plume  tilt  angle 

emx  *  angle  to  the  top  of  fuel  In  adjacent  unignited  cell. 


Transfer  Is  then  made  to  label  (L). 
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At  label  (J)  the  plume  is  not  tilted  enough  for  the  plume 
centerline  to  lie  below  the  top  of  the  fuel  in  the  subject  cell.  In  this 
case  the  plume  is  required  to  extend  across  the  Intervening  street  and 
overlap  the  fuel  in  the  adjacent  cell  for  ignition  to  be  considered.  If 
overlap  occurs,  transfer  is  to  label  K,  else  transfer  is  to  label  L. 

At  label  (K)  the  probability  of  Ignition  of  the  subject  cell  due 
to  radiation  from  the  overlapping  plume  of  the  burning  cell  is  calculated 
as 


(5-39) 


where 


Ap/Ac  *  fraction  of  cell  (f)  for  which  plume  is  defined. 

Transfer  Is  then  made  to  label  (L). 

At  label  (I)  as  at  label  (J),  prior  to  the  transfer  to  label  L, 
the  probability  of  subject  cell  ignition  is  calculated  by  Equations  5-38 
and  5-39,  respectively  and  accumulated  for  each  appropriate  burning  cell 
adjacent  to  the  subject  cell.  This  accumulation,  when  performed  for  all 
cells  adjacent  to  the  subject  cell,  represents  the  spatial  accumulation  of 
probability  of  ignition  for  all  burning  cells  surrounding  the  subject 
cel  1 . 


At  label  (L)  a  test  is  made  to  determine  if  all  eight  cells 
adjacent  to  the  subject  cell  have  been  examined.  If  all  adjacent  cells 
have  not  been  examined,  transfer  is  to  label  (A)  where  the  next  cell  is 
selected,  else  transfer  ts  to  label  (M). 
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At  label  (M)  the  just  calculated  cumulative  probability  of  Igni¬ 
tion  over  the  cells  surrounding  the  subject  cell  is  accumulated  to  the 
probability  of  Ignition  calculated  at  the  last  subject  cell  update.  If 
the  minimum  value  of  x/tmx  for  the  eight  adjacent  cells  (T/%,x)mn 
Is  greater  than  four,  transfer  Is  to  label  (N),  else  transfer  is  to  label 
(0). 

At  label  (N),  (T/%x)mn  *s  greater  than  four,  which  means 
that  all  burning  cells  adjacent  to  the  subject  cell  are  beyond  the  time  of 
burnout.  Under  these  conditions,  no  Increase  will  occur  In  the  cumulative 
probability  of  subject  cell  ignition.  For  these  reasons,  the  time  inter¬ 
val  to  the  next  examination  of  the  subject  cell  by  ACUPDT  is  set  at  the 
maximum  allowable  value. 

At  label  (0),  (t/%x)  is  less  than  or  equal  to  four.  Thus,  at 
least  one  of  the. cells  adjacent  to  the  subject  cell  is  currently  burning. 
The  time  interval  to  the  next  examination  o'  the  subject  cell  by  ACUPDT  is 
scheduled  in  accordance  with  the  maximum  rate  of  change  of  the  heat  pro¬ 
duction  rates  of  the  cells  surrounding  the  subject  cell  by  the  relation 

DT  =  (1/2) (1  -  q  )(x  )  ,  (5-40) 

'  mv  mx'mn 

where 

V-[U-(t/^)m)]exp!,/,mx)m  .  (5-4!) 

Control  is  then  transferred  to  label  (W). 

At  label  (P)  the  just  calculated  probability  of  cell  ignition 
is  tested.  If  Pjg  is  zero,  transfer  is  to  label  (Q),  else  trans¬ 
fer  is  to  label  (R) . 


I 


At  label  (R)  the  just  calculated  cumulative  probability  of  cell 
Ignition  (Pclinl)  Is  tested.  If  Pcum  Is  equal  to  or  greater  than  unity, 
transfer  Is  to  label  (S),  else  transfer  Is  to  label  (T). 


At  label  (S),  Ignition  Is  certain  and  the  subject  cell  Is 
Ignited  by  a  call  to  BCINIT.  The  cumulative  probability  of  Ignition  Is 
allowed  to  exceed  unity  In  cases  where  Ignition  conditions  are  appropri¬ 
ately  severe.  The  value  of  Pcun  is  passes  to  BCINIT  In  such  cases  as 
the  number  (real  valued)  of  simultaneous  Ignitions,  thus  allowing  for 
multiple  (fractional)  ignitions.  The  effect  of  multiple  Ignitions  Is  to 
cause  the  cell  to  burn  hotter  and  faster  than  it  otherwise  would.  These 
effects  are  implemented  as  shown  In  Equations  5-12  by  making  ^  an 
Inverse  function  of  the  number  of  ignitions  N^g  in  the  cell  and  In 
Equation  5-8  by  making  qc  a  function  of  1/%X* 

At  label  T  the  cumulative  probability  of  ignition  for  the  sub¬ 
ject  cell  Is  greater  than  zero  but  not  certain.  The  determination  of  what 
to  do  in  this  case  is  based  on  a  definition  of  the  rate  of  change  of  cumu¬ 
lative  probability  of  cell  Ignition  per  unit  time.  The  cumulative  proba¬ 
bility  of  cell  Ignition  Pcum  Is 


cum 


N 

J-l 


<Vj 


(5-42) 


where  (Pjg)j  Is  the  probability  of  ignition  associated  with  the  j-th 
cell  update  time  interval  Atj  and  where  the  time  interval  ATcutl  over 
which  probability  of  ignition  has  been  accumulated  is 


AT 

cum 


N 

w 

l 

j-l 


(5-43) 
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The  increment  of  cumulative  probability  remaining  until  cell  ignition  is 
certain  ( aP ' )  is  given  by 


AP'  -  1 


cum 


(5-44) 


Denote  the  time  interval  until  ignition  is  certain  by  AT1  and  define  the 
average  rate  of  change  of  cumulative  probability  of  cell  ignition 
dPCuii)/dt  over  the  interval  aT'  by 


AP' 

AT1 


(5-45) 


Define  the  current  rate  of  change  in  cumulative  probability  of  cell  igni¬ 
tion  in  terms  of  the  most  recently  calculated  values  as 


d(pcj  .  ">ia)N 
dt  AT.. 


(5-46) 


At  the  current  rate  of  change,  the  time  interval  until  cell  ignition  is 
certain  is  given  by  Equations  5-45  and  5-46  as 

AT'  -  AT  ,  (5-47) 

ig  N 

which  is  called  the  projected  time  to  ignition.  The  interval  until  the 
next  scheduled  cell  update  AT^+j  has  already  been  calculated  at  either 
label  (N)  nr  label  (0).  If  aT'  is  less  than  aTn+j  transfer  is  to  label 
(U),  else  transfer  is  to  label  (V). 


At  label  (U)  the  projected  time  to  Ignition  Is  less  than  the 

time  to  the  next  scheduled  update.  The  cell  is,  thus,  scheduled  for 

Ignition  at  the  projected  ignition  time. 

At  label  (V)  the  time  to  the  next  cell  update  Is  less  than  the 
projected  time  to  Ignition.  The  cell  Is,  thus,  scheduled  to  be  updated  at 
the  calculated  update  time. 

At  label  (W)  all  events  that  have  been  scheduled  by  ACUPDT  are 
placed  in  the  Event  List  and  all  new  cell  and  fire  data  calculated  by 

ACUPDT  are  stored  In  the  Dynamic  Data  Base.  Control  is  then  returned  by 

ACUPDT  to  the  UFMGfL 

5.6  FIRE  UPDATE  (FUPDT)  NODULE 

The  logical  organization  of  the  Fire  Update  Module  (FUPDT)  Is 
shown  in  Figure  5-8.  The  FUPDT  module  deals  primarily  with  the  problem  of 
scheduling  the  next  update  of  the  subject  fire.  The  mechanics  of  the 
update  process  are  performed  by  the  Fire  Geometry  (FGEOM)  and  Fire  Physi¬ 
cal  Data  (FPHYSD)  subroutines,  which  will  be  described  next. 

During  initialization,  all  ignitions  need  not  occur  simultane¬ 
ously,  but  can  be  spread  over  a  specified  interval.  This  can  Initially 
cause  many  small  Independent  fires  with  rapid’y  changing  sizes  and  heat 
production  rates  which  can  result  in  the  necessity  for  many  fire  updates 
early  In  the  simulation.  These  updates  occur  over  time  intervals  that  are 
significant  to  each  fire  but  are  not  significant  over  all  of  the  fires. 
To  satisfy  the  tradeoff  between  frequency  and  significance,  all  fire  up¬ 
dates  are  constrained  ay  the  FUPDT  module  to  occur  at  Intervals  that  are 
consistent  with  significant  changes  In  global  fire  characteristics. 
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GET  DATA  OH  SUBJECT  FIRE  FROM  DYNAMIC  DATA  BASE 
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Fire  data  are  treated  In  two  parts,  physical  data  and  geometri¬ 
cal  data.  The  fire  physical  data  consist  of  the  burning  characteristics 
of  the  aggregate  of  cells  that  make  up  the  fire.  The  geometrical  data 
consist  of  a  partitioning  of  the  fire  into  a  set  of  suitably  symmetric 
subfires.  To  maintain  physically  significant  fire  data,  the  following 
fire  update  criterion  are  employed: 

1.  A  fire  physical  update  occurs  whenever  the  fire  heat 
production  rate  has  changed  by  ten  percent  relative 
to  the  value  at  the  last  update  time. 

2.  A  fire  geometry  update  occurs  whenever  the  fire  area 
has  changed  by  ten  percent  relative  to  the  area  at 
the  last  update  time. 

3.  A  fire  physical  data  update  can  occur  without  an 
accompanying  geometry  update. 

4.  A  fire  geometry  data  update  must  be  accompanied  by 
a  physical  data  update. 

5.  No  fire  may  be  updated  more  frequently  than  once 
each  360  seconds. 

The  period  of  360  seconds  was  chosen  as  the  standard  update  interval 
because  it  represents  10  percent  of  the  minimum  time  in  which  an  ignited 
cell  can  reach  maximum  rate  of  heat  production. 

Fire  updates,  thus,  occur  whenever  the  heat  production  rate  or 
the  size  of  the  fire  are  changed  by  ten  percent  relative  to  the  respective 
values  at  the  last  update.  The  Burning  Cell  Update  (BCUPDT)  module  is 
responsible  for  maintenance  of  heat  production  rate  values  for  burning 

/ 
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cells.  It,  thus.  Incrementally  changes  the  heat  production  rates  of  fires 
and  determines  if  a  fire  update  is  required  under  the  change  in  fire  heat 
production-rate  criterion.  In  the  process  of  updating  burning  cells,  the 
BCUPOT  module  also  identifies  cells  that  burn  out.  It,  thus,  decrements 
fire  size  as  cells  burn  out,  and  determines  if  a  fire  update  Is  required 
under  the  change-ln-f ire-area  criterion.  The  Ignition  (IGNITE)  module  Is 
responsible  for  igniting  cells.  It,  thus,  increments  fire  area  and  deter¬ 
mines  if  a  fire  update  is  required  under  the  change-ln-flre-area 
criterion. 


Upon  entry  to  the  FUPDT  module,  shown  in  Figure  5-8,  fire  data 
for  the  subject  fire  are  accessed.  If  the  fire  heat  production  rate  is 
greater  than  zero,  transfer  is  to  label  (A),  else  transfer  is  to  label 
(C). 


At  label  (A)  the  most  recent  update  of  the  subject  fire  is  com¬ 
pared  with  the  last  standard  update  time.  If  the  most  recent  update  was 
later  than  the  last  standard  update  time,  transfer  Is  to  label  (B),  else 
transfer  Is  to  label  (C) . 

At  label  (B)  the  fire  has  been  updated  since  the  last  standard 
update  time.  The  currant  update  is,  thus,  rescheduled  to  occur  at  the 
time  of  the  last  update  plus  the  standard  fire  update  interval.  Control 
is  transferred  to  label  (E). 

At  label  (C)  either  the  fire  heat  production  rate  is  zero  or  the 
last  update  was  at  or  before  the  last  standard  update.  In  the  first  case, 
the  initial  update  is  required  and  in  the  second  case,  the  next  update  is 
required.  In  either  case,  an  update  Is  currently  required,  and  the  next 
update  is  scheduled  to  occur  at  SIMTIM  plus  the  standard  fire  update 
interval . 
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At  label  (D),  the  geometry  and  physical  data  for  the  fire  are 
calculated  by  the  Fire  Geometry  subroutine. 

At  label  (E),  all  events  scheduled  by  the  FUPDT  module  are 
stored  in  the  Event  List,  all  data  are  stored  in  the  dynamic  data  base, 
and  the  FUPDT  module  returns  control  to  the  UFMGR  module. 

5.6.1  Fire  Geometry  (FGEOM)  Subroutine 

The  algorithm  used  for  fire  induced  wind  velocity  calculations 
assumes  that  the  fires  involved  are  circular.  In  the  processes  of  ini¬ 
tialization  and  propagation,  high  aspect  (length-width)  ratio  fires  can  be 
generated,  which  violate  this  assumption.  The  fire  geometry  subroutine 
employs  a  simple  pattern  recognition  algorithm  to  identify  long  thin  fires 
and  partitions  them  into  sets  of  approximately  square  subfires,  as  a  suit¬ 
able  representation  of  the  circular  fires  expected  by  the  fire  induced 
wind  model . 

Upon  entry  to  the  FGEOM  subroutine  shown  in  Figure  5-9,  the 
existing  physical  and  geometrical  fire  data  for  the  subject  fire  are 
deleted.  The  overall  fire  physical  data  are  independent  of  the  fire  geom¬ 
etry.  Thus,  the  Fire  Physical  Data  (FPHYSD)  subroutine  is  entered  to 
calculate  physical  data  for  the  entire  fire.  Upon  return  from  the  FPHYSD 
subroutine,  the  number  of  cells  in  the  fire  is  tested.  If  the  fire  con¬ 
tains  four  or  more  cells,  transfer  is  to  label  (A),  else  transfer  is  to 
label  (J). 


At  label  (A),  the  fire  has  been  found  to  have  four  or  more 
cells.  Fires  smaller  than  four  cells  are  not  examined  for  elongation. 
This  is  consistent  with  the  attempt  to  maintain  the  aspect  ratio  of 
involved  fires  to  be  no  greater  than  four.  The  number  of  candidate  sub¬ 
fires  is  set  to  zero  (i.e.,  the  entire  fire),  and  all  cells  in  the  fire 
are  assigned  to  this  subfire. 


FIRE  ID 
SIMTIM 


DELETE  EXISTING  FIRE  PHYSICAL  AND  GEOMETRY  DATA 


CALL  PPHYSD 


UPDATE  FIRE  PHYSICAL  DATA 


DOES  FIRE  CONTAIN  FOUR  OR  MORE  CELLS? 


SET  NUMBER  OF  CANDIDATE  SUBFIRES  TO  ZERO  AND  ASSIGN  ALL 
BURNING  CELLS  TO  THE  ZEROTH  SUBFIRE  (THE  WHOLE  FIRE) 


Figure  5-9.  Logic  flow  of  fire  geometry  subroutine. 


1 


Figure  5-9.  Logic  flow  of  fire  geometry  subroutine  (continued) 
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YES 


0 


THE  ASPECT  RATIO  OF  THE  BOUNDING  ROWS  AND 
COLUMNS  GREATER  THAN  FOUR? 


Figure  5-9.  Logic  flow  of  fire  geometry  subroutine  (continued) 
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HAVE  ALL  CANDIDATE  SUBHRES  BEEN  EXAMINED? 


Figure  5-9.  Logic  flow  of  fire  geometry  subroutine  (concluded). 


At  label  (B)  the  cells  for  the  next  candidate  subfire  are  acces¬ 
sed.  At  the  first  entry  to  label  (B),  the  only  candidate  is  the  whole 
fire  (i.e.,  zeroth  subfire),  and  all  of  the  cells  in  the  fire  are  acces¬ 
sed.  A  subfire  consists  of  a  set  of  contiguous  burning  cells,  cells  are 
identified  by  their  row-column  Indices,  and  the  rectilinear  boundary  of  a 
subfire  is  defined  by  the  maximum  and  minimum  row  and  column  indices 
exhibited  by  the  set  of  cells  in  the  candidate  subfire.  These  Indices 
identify  the  rows  and  columns  that  bound  the  candidate  subfire. 

At  label  (C)  the  area  within  the  bounding  rows  and  columns  Is 
calculated.  If  the  burning  cells  constitute  less  than  half  the  bounded 
area,  transfer  is  to  label  (D),  else  transfer  is  to  label  (F). 

At  label  (0),  if  the  total  number  of  cells  within  the  bounding 
rows  and  columns  is  greater  than  ten,  transfer  is  to  label  (E),  else 
transfer  is  to  label  (H). 

At  label  (E)  the  bounded  area  is  less  than  half  filled  by  the 
burning  cells  and  the  total  number  of  burning  cells  Is  greater  than  ten. 
The  intent  here  is  to  identify  subfires  represented  as  sparsely  filled 
bounded  areas.  The  long  dimension  of  these  candidate  subfires  is  approxi¬ 
mately  halved  along  cell  boundrles  to  create  two  potentially  less  sparsely 
filled  subfires.  A  simple  example  of  such  a  fire  is  shown  in  Figure  5-10, 
where  it  can  be  seen  in  (a)  that  the  burning  cells  (shaded)  constitute 
only  30  percent  of  the  bounded  cells.  The  largest  dimension  of  the  fire 
is,  thus,  partitioned  as  shown  in  Figure  5-10(b),  creating  candidate  sub¬ 
fires  1  and  2. 

The  partitioning  of  the  example  fire  shown  In  Figure  5-10  will 
be  described  here,  although  some  of  the  logic  will  not  be  encountered 
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Figure  5-10.  Example  of  fire-subfire  partitioning. 
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until  later.  One  of  the  reasons  for  re-examining  each  newly  created  can¬ 
didate  subfire  by  a  subsequent  entry  to  label  (B)  is  clearly  shown  in 
Figure  5-10(c)  and  is  as  follows.  As  shown  In  Figure  5-10(b),  the  burning 
cells  of  candidate  subfire  1  only  occupy  25  percent  of  its  portion  of  the 
bounded  area.  However,  when  the  boundary  for  candidate  subfire  1  is  rede¬ 
fined  at  label  (C),  the  burning  cells  constitute  100  percent  of  the 
bounded  area  as  shown  in  Figure  5-10(c),  Candidate  subfire  1  will  be 
identified  as  an  elongated  subfire  and  appropriately  subpartitioned.  Can¬ 
didate  subfire  2  of  Figure  5-I0(b)  has  only  30  percent  burning  occupancy 
and  is  subpartitioned  into  candidate  subfires  2A  and  2B.  In  a  subsequent 
entry  to  label  (B),  candidate  subfire  2A  will  be  found  acceptable  with 
respect  to  burning  cell  occupancy  and  elongation,  and  will  be  made  into  a 
subfire.  Candidate  subt'ire  2B  will  be  found  unacceptable  with  respect  to 
occupancy  (33  percent),  but  too  small  (less  than  10  cells)  to  subdivide, 
and  will  also  be  made  into  a  subfire. 

The  candidate  subfire  under  examination  has  been  subpartitioned, 
the  new  candidate  subfires,  thus,  created  have  been  assigned  IDs,  and  the 
total  number  of  subfires  Included  in  this  fire  has  been  appropriately 
incremented.  At  this  point,  transfer  is  back  to  label  (B)  where  the  next 
subfire  will  be  accessed  ard  the  above  described  tests  will  be  repeated. 

At  label  (F)  the  bounded  area  of  the  candidate  subfire  is  50 
percent  or  more  filled  with  burning  cells.  The  ratio  of  the  longer  to  the 
shorter  rectilinear  boundary  lengths  (aspect  ratio)  is  calculated  for  the 
candidate  subfire.  If  the  aspect  ratio  is  greater  than  four,  transfer  is 
to  label  (G),  else  transfer  is  to  label  (H). 

At  label  (G)  an  elongated  candidate  subfire  has  been  identified 
and  the  degree  of  elongation  quantified.  The  candidate  subfire  is  parti¬ 
tioned  into  subfires  each  having  an  aspect  ratio  no  greater  than  two.  An 
example  elongated  subfire  is  shown  in  Figure  5-10(c)  as  subfire  1.  It  can 
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be  seen  that  subfire  1  has  an  aspect  ratio  of  six  and  has  been  partitioned 
Into  three  subfires  (1A,  ID,  and  1C)  each  having  aspect  ratios  of  two. 
After  the  new  subfires  have  been  defined  and  the  number  of  subfires  appro¬ 
priately  incremented,  transfer  is  back  to  label  (B)  where  the  next  subfire 
is  accessed. 

At  label  (H)  the  candidate  subfire  has  failed  the  occupancy  test 
at  label  (C),  failed  the  size  test  at  labef  (0),  or  failed  the  elongation 
test  at  label  (F).  In  any  of  these  cases,  the  candidate  requires  no  fur¬ 
ther  modifications  to  become  a  subfire.  The  cells  are  assigned  to  a 
subfire  and  physical  fire  data  are  calculated  by  the  FPHYSD  subroutine  for 
this  subfire.  Transfer  Is  then  to  label  (I). 

At  label  (I),  if  all  candidate  subfires  that  have  been  identi¬ 
fied  have  also  been  examined,  transfer  is  to  label  (J)  where  control  Is 
returned  to  the  FUPDT  module,  else  transfer  is  to  label  (B)  where  the  next 
subfire  is  accessed. 

5.3.2  Fire  Physical  Data  (FPHYSD)  Subroutine 

The  Fire  Physical  Data  subroutine  calculates  current  values  for 
the  specified  physical  parameters  of  a  fire  or  subfire.  The  fire  physical 
(and  geometrical)  data  are  an  exception  to  the  rule  of  maintaining  two 
sets  of  data  subtending  current  time  and  Interpolating  between  them  to 
obtain  current  values  of  specified  variables.  It  is  not  possible,  in  this 
case,  because  changes  in  the  fire  geometry  may  not  allow  correlation 
between  the  subfires  existing  at  different  times.  Fire  physical  data  are, 
thus,  updated  more  frequently  than  Interpol aed  data,  but  are  held  constant 
between  updates. 

The  logic  flow  of  the  FPHYSD' subroutine  Is  shown  In  Figure  5-11, 
where  it  can  be  seen  to  be  devoid  of  logic  branches  and  to  consist  simply 
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Figure  5-11.  Logic  flow  of  fire  physical  data  subroutine. 
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of  a  sequence  of  calculations.  Upon  entry  to  the  FPHYSD  subroutine,  the 
list  of  cells  assigned  to  the  subject  fire  (or  subfire)  Is  accessed  and 
the  total  heat  production  rate  Is  calculated  from 

Nk 

L  •  l  q,  (0/s)  (5-48) 

K  1*1  1 

where 

• 

*  total  heat  production  rate  of  fire  k 
N|<  *  number  of  cells  In  fire  k 

q-j  •  heat  production  rate  of  the  1-th  cell  In  fire  k. 

The  location  of  the  effective  fire  center  Is  calculates  In  terms 
of  the  fire  extent  weighted  by  local  heat  production  rate,  and  Is  used  In 
the  fire  Induced  wind  velocity  calculations.  The  coordinates  of  the  ef¬ 
fective  center  of  fire  k,  (TT^,  7^),  are  given  by 

\  \ 

CV>k) '  h  Vi  •  Vi  (■)  (5-49> 

where  x^.y^  a  coordinates  of  the  center  of  the  i-th  cell  of  fire  k. 

The  effective  radius  (r^  of  fire  k,  also  used  In  fire 
Induced  wind  velocity  calculations.  Is  calculated  from 

(  Nk  ) 

(r„>i<  "  {jp-  £  [(*i  •  "n!2  +  <*i  -  lijl/2  (")•  (5-50> 

The  volume  rate  of  air  Inflow  (V^)  to  fire  k  Is  calculated  from 


\  *  aQ^/bc  (m3/sec)  (5-51) 

where 

a  ■  mass  of  air  required  to  burn  a  mass  of  fuel, 

(kg  air/kg  fuel) 

b  -  heat  produced  per  weight  of  fuel  (J/kg) 
c  »  air  density  (kg/m3). 

The  area  [\)  of  fire  k  Is  given  by 

A,  -  N.A  ,  (5-52) 

k  k  c 

where  Ac  is  the  specified  cell  area.  The  cell  area  cannot  be  changed 
during  an  execution,  but  may  be  changed  between  executions.  The  time 
these  data  were  calculated  Is  denoted  as  SIMTIM,  and  control  Is  returned 
to  the  calling  routines. 

5.7  BRANDING  (BRANDG)  NODULE 

The  execution  of  a  burning  cell  update  event  consists  of 
sequentially  p^oceslng  the  BCUPDT  module  (Subsection  5.4)  and  the  BRANDG 
module.  Thus,  upon  entry  to  the  BRANDG  model,  current  BCUPDT  data  are 
available. 

The  logic  flow  diagram  of  the  BRANDG  module  Is  showr  In  Figure 
5-12,  where  the  event  data  are  shown  to  Include  the  (burning)  cell  ID, 
SIMTIM,  cell  heat  production  rate,  and  the  time  of  the  next  cell  update. 
Only  burning  cells  can  launch  brands,  and  this  Is  assured  by  making  the 
test  for  branding  a  part  of  each  burning  cell  update  event.  The  time  of 
the  next  cell  update  Is  calculated  by  the  BCUPDT  module  and  passed  to  the 
BRANDG  module. 
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Figure  5-12.  Logic  flow  for  branding  nodule  (continued) 
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Figure  5-12.  Logic  flow  for  branding  Module  (concluded). 
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Upon  entry  to  the  BRANDG  module,  SIMTIM  is  tested,  and  If  SIMTIM 
Is  greater  than  1000  seconds,  transfer  is  to  label  (A),  else  transfer  is 
to  label  (J)  where  control  is  returned  to  the  UFMGR.  Branding  is  not 
allowed  for  the  first  1000  seconds  (■  17  min)  of  a  simulation  for  both 
physical  and  computational  reasons.  In  the  cases  that  have  been  run, 
initialization  (i.e.,  initial  Ignitions)  have  been  allowed  to  occur  over 
the  first  1000  seconds  , of  simulation  time.  Thus,  cells  are  being  ignited 
throughout  this  Interval.  In  general,  large  numbers  of  cells  are 
initialized,  and  all  have  very  low  probabilities  of  branding  over  the 
first  1000  seconds.  Significant  computational  resources  are  saved  by  not 
examining  these  low  probability  occurrences. 

At  label  (A),  If  the  subject  cell  is  a  peripheral  cell  of  the 
fire  to  which  it  belongs,  transfer  is  to  label  (B),  else  transfer  is  to 
label  (J) .  Only  peripheral  cells  are  tested  for  branding.  It  is  assumed 
that  brands  launched  by  interior  cell  plumes  will,  if  they  fall  out,  be 
caught  up  in,  and  consumed  by  the  plume  of  an  adjacent  burning  cell. 

At  label  (8)  the  time  Interval  until  the  next  cell  update  At  is 
calculated,  and  the  number  of  brands  that  will  be  produced  In  the  Interval 
until  the  next  cell  update  Njj  is  estimated  as 

Nb  *  q  At/ (1.63  x  1011)  ,  (5-53) 


where  q  ■  cell  heat  production  rate.  The  constant  Includes  the  nature  of 
the  fuel  in  the  cell,  which  can  be  made  explicit.  The  effective  heat 
production  rate  per  unit  cell  area  is  calculated  from 
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where  Lc  B  the  cell  linear  dimension.  From  the  cell  ID,  the  coordinates 
of  the  cell  center  are  calculated. 

At  label  (C)  the  wind  velocity  (calculated  in  the  BCUPDT  module) 
is  accessed.  A  loop  is  formed  over  the  estimated  number  of  brands.  Each 
brand  is  characterized  by  its  terminal  velocity,  which  is  obtained  by  a 
random  draw  from  a  linear  distribution.  A  linear  distribution  fits 
empirical  data  and  produces  many  more  smaller,  than  larger,  brands.  From 
the  terminal  velocity  of  the  k-th  brand  the  brand  flight  time  \ 
is  calculated  as 

Tk  -  q(0 .444  nk  +  1.10  f*  -  0.975  r^f')  (s)  ,  (5-55) 

where 

fk  -  1.0  -  (1  -  fk)l/2  (5-56) 

and 

fk  -  4.726  x  10-4(1  -  0.887  t^)“2(1  -  0.482  r^)**3  .  (5-57) 

The  time  at  which  the  brand  will  burn  out  is  calculated  as 

Tbo  -  q( 1.1000  -  0.5314  r^)  (s)  ,  (5-58) 

and  the  amount  of  time  the  brand  will  burn  after  impact  is  given  by 
tBI  3  ^BO  "  Tk  (s)*  The  distance  the  brand  is  transported  from  the 
cell  center  is  given  by 
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(5-59) 


Ar  «*  vGTk  (m)  , 

where  vg  *  surface  wind  velocity  at  cell  center. 


At  label  (£),  if  the  the  brand  transport  distance,  jar|,  is  such 
that  impact  is  within  the  launch  cell,  transfer  is  to  label  (I),  else 
transfer  is  to  label  (F). 

At  label  (F)  the  Impact  coordinates  anc!  the  ID  of  the  impact 
cell  are  calculated.  If  the  impact  cell  Is  already  ignited  or  burned  out, 
transfer  is  to  label  (I),  else  transfer  is  to  label  (G). 

At  label  (G)  the  probability  that  the  brand  will  ignite  the 
impact  cell  is  calculated  by 

P.  »  1  -  exp(-TDT/60)  .  (5-60) 

ig  til 

If  the  cell  is  to  be  ignited,  transfer  is  to  label  (H),  else  transfer  is 
to  label  (I). 

At  label  (H)  an  event  is  added  to  the  Event  List,  which  will 
ignite  the  impact  cell  at  SIMTIM  +  .  At  label  (I),  if  all  (Ng) 
brands  have  been  tested,  transfer  is  to  label  (J),  else  transfer  is  to 
label  (D)  where  the  terminal  velocity  is  drawn  for  the  next  brand.  At 
label  (J)  events  generated  by  the  BRANDG  module  are  stored  on  the  Event 
List  and  control  is  returned  to  the  UKMGR. 
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Figure  5-13.  Logic  flow  of  typical  HAP  module. 
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5-8 


OUTPUT  NODULES  (WVMAP,  CSNAP,  )MAP) 


The  output  modules  present  the  results  of  the  simulation  in  the 
form  of  maps  of  wind  velocity  (WVMAP),  cell  state  (CSMAP),  and  heat  pro¬ 
duction  rate  (QMAP)  over  the  specified  urban  area.  The  logic  flow 
diagrams  of  these  modules  are  very  similar  and  are,  thus,  represented  in 
Figure  5-13  by  a  logic  flow  diagram  typical  of  the  three  modules. 

At  label  (A)  each  cell  of  the  specified  urban  area  (matrix)  is 
identified  by  row  and  column.  The  map  quantity  (i.e.,  wind  velocity,  cell 
state,  or  heat  production  rate)  is  calculated  for  the  identified  cell  and 
written  to  a  plot  file.  This  Is  continued  until  a  value  of  the  map  quan¬ 
tity  has  been  calculated  and  stored  for  every  cell  In  the  urban  area. 

At  label  (B)  specified  maximum  and/or  minimum  values  of  the  map 
quantity  are  written  to  the  plot  file.  The  map  module  then  stores  its 
next  update  event  in  the  Event  List.  The  update  event  is  timed  according 
to  a  specified  schedule,  which  may  be  periodic  or  aperiodic. 
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SECTION  6 

EXAMPLE  MODEL  RESULTS 


Results  from  the  model  for  eleven  example  cases  are  presented  in 
this  section.  The  results  for  each  case  presented  consist  of  the  follow¬ 
ing: 

1.  A  description  of  the  initial  conditions  for  the  fire 

2.  Sets  of  maps  depicting  these  results: 

a)  the  state  of  each  cell  in  the  area 
bi  heat  production  rate  contours  over  the  area 

c)  heat  production  rates  of  each  cell  in  the  area 

d)  wind  velocity. 

Map  results  are  presented  at  times  of  1000,  1750,  2750, 
4750,  8750,  16750,  and  32750  seconds  after  simulation 

initialization. 

Included  with  the  results  for  each  case  is  a  description  of  the  ini¬ 
tialization  conditions  and  brief  comments  on  the  case  results. 

Initialization  is  an  important  aspect  of  uncontrolled  urban 
fires  in  the  cases  of  interest.  For  this  reason,  and  because  of  its 

uncertainty,  the  model  has  been  designed  to  accept  arbitrary  Initializa¬ 
tion  processes  For  the  example  cases  presented,  the  initialization 
consisted  of  selecting  a  candidate  set  of  cells  for  ignition  and  Igniting 
them  over  a  1000  second  Interval  with  a  uniform  random  distribution  in 

time  anL  a  Poisfon  distribution  of  from  zero  to  four  ignitions  per  candi¬ 

date  cell  for  a  specified  expected  value  of  ignition. 


All  example  cases  are  simulated  in  an  urban  area  that  Is  50  x  50 
city  blocks  (about  56.25  km2).  Cellulosic  fuels  loaded  at  40  kg/m2  are 
uniformly  distributed  over  the  area.  The  model  will  allow  different  fuel 
types  and  loadings  In  each  cell;  however,  for  the  example  cases,  it  was 
desired  to  have  the  results  reflect  phenomenology  variations  Instead  of 
fuel  variations. 

An  objective  of  developing  the  demonstration  model  was  to  demon¬ 
strate  a  model  structure  capable  of  integrating  the  range  of  variables  and 
variable  interactions  required  to  describe  the  onset  and  growth  of  mass 
fires  in  an  urban  setting.  Even  for  demonstration  purposes,  the  model 
must  contain  physically  based  phenomenology  models  that  include  most  of 
the  variables  and  variable  interactions  of  interest.  In  this  sense  the 
phenomenology  models  included  in  this  demonstration  urban  fire  model  are 
considered  j  be  qualitatively  appropriate,  although  it  is  known  that  they 
are,  in  many  cases,  quantitatively  in  error. 

In  the  interest  of  aiding  in  model  development  and  in  mitigating 
model  maintenance  problems,  the  model  structure  has  been  designed  to  read¬ 
ily  allow  the  introduction  of  new  models  and/or  processes  and  for  the 
alteration  of  existing  models  and/or  processes. 

6.1  LARGE  AREA  FIRE  WITH  NO  AMBIENT  HIND 

As  shown  in  Figure  6-1,  the  central  35  x  35  cells  of  the  50  x  50 
cell  urban  area  are  chosen  as  candidates  for  ignition,  and  ignited  with  a 
probability  of  ignition  of  0.1.  There  is  initially  no  ambient  wind  veloc¬ 
ity;  however,  very  strong  fire-induced  winds  are  generated,  as  shown  in 
Figures  6-1(1)  and  6-l(m).  It  is  noted  that  the  wind  velocities  are  con¬ 
sidered  to  be  excessive  in  all  of  the  cases  presented. 
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Figures  6-1  through  6-1 (g) ,  together  with  Figures  6-l(j)  and 
6-1(1),  show  that  the  fire,  once  ignited,  tends  to  propagate  inwards,  fol¬ 
lowing  the  inward  fire-generated  wind  flow.  These  figures  also  show  that, 
as  the  fire  grows  inward,  it  generates  radial  "fingers"  along  the  wind 
velocity  stream  lines.  This  is  consistent  with  wildland  fire  experience, 
which  reveals  similar  combustion  patterns  under  high  wind  velocity  condi¬ 
tions. 


Figures  6-1  and  6-1 (c)  show  a  quantitative  phenomenology  that 
occurs  in  other  example  cases  as  well.  Looking  at  the  upper  right  hand 
corner  of  Figure  6-1,  it  can  be  seen  that  no  cells  are  ignited  above  row 

44.  In  Figure  6-l(c),  it  can  be  seen  that  a  cell  has  been  ignited  in  ro.,» 

50.  Because  this  cell  is  disjoint  from  any  other  Ignited  cell,  ignition 
had  to  have  occurred  by  branding.  However,  as  shown  by  Figures  6-1(1)  and 
6-1  (m),  the  winds  are  all  directed  inward  toward  the  center  of  the  urban 
area.  In  view  of  this,  the  only  way  the  cell  in  row  50  could  have  been 

ignited  by  branding  Is  for  a  brand  to  have  been  thrown  from  the  lower  left 

corner  of  the  fire  all  of  the  way  across  the  burning  area.  This  is 
clearly  erroneous  and  is  an  artifact  of  the  simplicity  of  the  brand  trans¬ 
port  model  that  is  used. 
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Figure  6-1.  Cell  state  Map  for  area  fire. 
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Figure  6~l(a).  Ceil  state  nap  for  area  fire 
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Figure  6-1 (b) .  Cell  state  nap  for  area  fire. 
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Figure  6-lfa).  Cell  state  «ap  for  area  fire 
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Figure  6-l(d).  Cell  state  *ap  for  area  fire. 
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Figure  6-1 (e) .  Cell  state  »ap  for  area  fire. 
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Figure  6-l(f).  Cell  state  Map  for  area  fire. 
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Figure  6-1 (g) .  Cell  state  Map  for  area  fire 
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Figure  6-1(1).  Heat  production  rate  contour  for  area  fire. 
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Figure  6-l(j).  Cell  heat  production  rate  Hap  for  area  fire. 
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Figure  6-1 (k) .  Cell  heat  production  rate  asp  for  area  fire. 
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Figure  6-1 (a) 


Wind  velocity  aap  for  are*  fire 


6.2 


LARGE  AREA  FIRE  WITH  10  M/SEC  AMBIENT  WIND  VELOCITY 

The  Initial  fire  geometry  In  this  case  Is  exactly  the  same  as  in 
Case  6.1.  The  ambient  wind  is  10  m/sec,  and  directed  from  the  lower  left 
corner  to  the  upper  right  corner.  Figures  6-2  through  6-2(g)  show  that 
the  fire  burns  very  similarly  to  that  of  Case  6.1.  The  spread  is  slightly 
more  rapid  and  tends  to  follow  the  wind  direction  more  than  in  Case  6.1. 
The  reason  the  fires  are  so  similar  is  because  the  fire  Induced  winds  are 
a  factor  of  from  two  to  four  times  as  large  as  the  ambient  wind.  Propaga¬ 
tion  effects,  which  are  strongly  related  to  wind  velocity,  are  thus 
dominated  by  the  fire  Induced  winds.  Figures  6-2(1)  and  6-2(m)  show  the 
interactions  between  the  ambient  and  fire  Induced  winds. 
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Figure  6-2.  Cell  state  aap  for  area  fire  irith  10  Vs  wblent  wind 
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Figure  6-2 (a) .  Cell  state  aap  for  area  fire  with  .10  m/s  aablent  wind 
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Figure  6-2(b).  Cell  state  wap  for  area  fire  with  .10  m/s  awblent  wind. 
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Figure  6-2(c).  Cell  state  Map  for  area  fire  with  .10  m/s  art  lent  wind 
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Figure  6-2(f).  Cell  state  nap  for  area  fire  with  .10  m/s  wbient  wind. 
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Figure  6-2(j).  I  11  heat  production  rate  map  for  area  fire  with  .10  m/s 
amolent  wind. 


189 


iMWtw  wiWti!  cifrwMrro* 

UMMI  MM  HiML 
r*  tiouawTiM  ■  4w*  mcmm 

M*M  I  M  H,  CfcKM  I  II  M. 
mi  mh  utttlM  imun  «rv  tr  mi  m.n% 

—  mi  nwi  mu  «u  mt  m  mm  u  iwtr?*  mm 
wvr  nu  is  LOMMirwnrc  mm  10  i  vi.w  1 


4* 

4  3  4  4  4 

+  44-44 


;  3  4- 
14 


4  4  4  4  4-  3  4  4  3  3  3  4  3 

4444.3  43433  3  4 
444444543333332 

544  4  444443+  0  3  3  J 
444+44444443444 
4444444444444- 5  4  4 

+  43  444  +  444444  + -3  45 
44454444444  4  4454  3 
44+ +4 44444 +44+3 534 

y  liiMtmh 

+4 5  5  44  4+4  544 
43+  5  443544 

5 


imm  NteMi 

uwn  p  sm  moil 


MM  l  II  M,  MM  l  t*  M; 

•Mt  m i  imw  1  tv  1  in  ■  mm  v  wr  mm 

—  Ml  MM*  mw  NU  *M  M  «M«M  |l  PlinU  «■— 

r  nu  (9  UMMirHmc  mr  io  t  o-ur*  j 


:|342f4 


mn 


4 


4IU4 

n4 

44  .  4 


43 

tt++44++++44 
4444.44  3  43  + 
44444I44+43+ 
444444+441+33 
444  +  A444  +  444  .  ,  .  ,  . 
4+4444444444+44444 
4 


,114+4+4+4444+3+ 

■  44  ■ .  , 

444444  +  444  44  4  +  44-4  + 


4-  +  +  +  5  +  4  +  +  +  4  +  +  +  45+4 

4  +  +  4'H  +  +-  +  4‘H  +  44S4''- 

*hnlhhlhtlnt 

34+4+44+544+44444 
4444+  4+  +4445+4 
4+4  5  444544 

+ 


1357913579135791357913579 


1357913579135791357913579 


*tt*ic*  hcmmch  cnMMTiM 

UfttftM  MM  H0DCL 


IIMIM  MMMMM  CMfMWTTM 

UMM  MM  HMML 


i  ti 


ftov  iTUin m  tim  - 
mm  in  u.  m 
M*l  MI  UTTMt  v  «v  V  m  •  I 
— «  '.Ml  mu  «TW  (umm  mm  <«  n.mw  — ■ 
wtr  khl*  i3  uMMinvttc  m  to  i  o-i.tv  t-i.iu  » 


+  3 

4+444443444 


4 

4 


ntwuiimn  *t\h 

4  +  +  +  +  + +  +  +  4+  +j  +  +  +f 

444444444334  44444 

444444444444  4+44 

444444444444434 
444  +  4+  +  +  +  ++  +  +  +  +  4  +  4 
A 


4444+444444+4+44443 
4444444444444+4+4+4 
4  +++4+44+44+4444+4 
++++++++++++++  4+ 
++++++++++++++4++ 
4444444444+444+44 
4 .144  44  44444  4 

44  - 


I  10  <  0-1.  t*  9*1.  Ill  J 

t  +  44  4 
£  4  +4  +  444 
o  4+4+44 

4444+44 


vnor  want  i*  itomitwhc  i 

444+44443 
344444+4444' 

+  .144  +  44+  4' 

4+444++444- 
4  +  +  4  444  +  4- 

444444  44- 
4+444444444444+44 

444444444+4444+444 

4+4+4+  4+444444+444 

+4+++++4+444+44++44 
++++4+44+4444  4444 
44+4+++  ++++4  4+4445 

++++++++4+++4++4  +4+4 

4  +444444444444+4+4+ 
4+4+  444+4+4+4+  4+4 
+  44  +  44.444.4441444  4  4 

4  4  4444444444+4444 

+444+  +44444  4++ 

44  414444  444+44444 
444  4  44444444  44+ 

444+1+  4  44+  44 


13579 135791 3575 1 3579 13579' 


.1357913579135791357913579 


Figure  6-2(k).  Cell  heat  production  rate  for  area  fire  with  .10  m/s 
ambient  wind. 
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Figure  6-2(1).  Wind  velocity  nap  for  area  fire  with  .10  */s  anbient  wind. 
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Figure  6-2(»).  Hind  velocity  nap  for  area  fire  with  10  n/s  anbient  wind. 
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6.3 


LINE  FIRE  WITH  NO  AMBIENT  WIND 


As  shown  in  Figure  6-3,  the  central  3  x  35  cells  were  ignited 
with  a  probability  of  ignition  of  0.5.  Fire  induced  winds  are  generated 
ranging  from  7  to  18  m/s  as  shown  in  Figures  6-3(1)  and  6-3(m).  This  is 
not  sufficient  wind  speed  to  cause  a  significant  amount  of  fire  propaga¬ 
tion.  The  ignited  area  does,  however,  partially  fill  in,  due  primarily  to 
contagion.  The  "dumbbell"  appearance  shown  in  Figures  6-3(f)  and  6-3(g) 
is  an  artifact  of  brands  being  tossed  completely  across  the  fire  as  des¬ 
cribed  in  Case  6.1. 
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Figure  6-3.  Cell  state  map  for  line  fire 
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Figure  6-3(a).  Cell  state  map  for  line  fire. 
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Figure  6-3(b).  Cell  state  map  for  line  fire 
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Figure  6-3(c).  Cell  state  map  for  line  fire. 
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Figure  6-3(d).  Cell  state  Map  for  line  fire. 
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.  Cell  state  map  for  line  fire. 


199 


Figure  6-3(e) 
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Figure  6-3(f).  Cell  state  Map  for  line  fire 
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Figure  6-3{g).  Cell  state  aap  for  line  fire. 
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Figure  6-3(h).  Heat  production  rate  contour  for  line  fire. 
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Figure  6-3(1).  Heat  production  rate  contour  for  line  fire. 
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Figure  6-3(4).  Cell  heat  production  rate  aap  for  line  fire. 
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Figure  6-3(k).  Cell  heat  production  rate  map  for  line  fire. 
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Figure  6-3(1).  Wind  velocity  map  for  line  fire. 
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Figure  6-3(»).  Mind  velocity  map  for  line  fire. 
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LINE  FIRE  WITH  10  M/SEC  AMBIENT  UIND 


The  ignition  geometry  is  the  same  as  for  case  6.3;  however,  the 
location  of  the  candidate  cells  was  moved  to  the  lower  left  corner.  This 
was  done  so  that  the  fire  would  have  room  to  propagate  along  the  wind 
direction  from  the  lower  left  corner  to  the  upper  right  corner.  Figures 
6-4  to  6-4(g)  show  that  the  fire  propagates  much  more  rapidly  than  in  Case 
6.3  and  that  the  propagation  is  along  the  wind  stream  lines.  Comparison 
of  the  wind  maps  for  Cases  6.3  and  6.4  shows  that  the  ratio  of  the  fire 
induced  to  the  ambient  winds  is  such  that  the  10  m/s  anbient  winds  would 
be  expected  to  cause  the  fire  to  burn  significantly  differently  than  in 
the  case  with  no  ambient  wind. 
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Figure  6-4.  Cell  state  nap  for  line  fire  with  .10  n/s  ambient  wind. 
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Figure  6-4 (a).  Cell  state  Map  for  line  fire  with  .10  «/s  anblent  wind. 
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Figure  6-4(b).  Cell  state  wap  for  line  fire  with  .10  «/s  anient  wind. 
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Figure  6-4(c).  Cell  state  map  for  line  fire  with  .10  a/s  aiblent  wind. 
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Figure  6-4(d).  Cell  state  nap  for  line  five  with  .10  m/s  wbient  wind 
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Figure  6-4(e) .  Cell  state  nap  for  line  fire  with  .10  m/s  anblent  wind. 
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Figure  6-4(f).  Cell  state  nap  for  line  fire  with  .10  n/%  anblent  wind 
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Figure  6-4(g).  Cell  state  Map  for  line  fire  with  .10  m/s  aoblent  wind 
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Figure  6-4(h).  Heat  production  rate  contour  aap  for  line  fire  with 
.10  n/s  abient  wind. 
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Figure  6-4(1).  Heat  production  rate  contour  wap  for  line  fire  with 
.10  m/s  afalent  wind. 
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Figure  6-4(1).  Wind  velocity  for  line  fire  with  .10  m/s  aMbient  wind. 
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Figure  6-4(*).  Hind  velocity  for  line  fire  with  .10  m/s  ant lent  wind. 
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6.5 


RING  FIRE  WITH  NO  AMBIENT  WIND 


As  shown  in  Figure  6-5,  a  four-cell  wide  annulus,  42  cells  In 
diameter,  was  Initially  ignited  with  a  probability  of  ignition  of  0.5. 
There  was  no  initial  wind  velocity;  however,  strong,  inward-directed  fire 
Induced  winds  resulted  in  heavy  branding  and  contagion  In  the  interior  of 
the  annulus.  It  is  clear  in  all  of  the  cases  that  fire  spread  is  much  too 
rapid.  However,  this  is  graphically  shown  by  the  marked  increase  in  the 
number  of  ignited  cells  In  the  1000  seconds  (from  1750  to  2750  seconds), 
as  showi  in  Figures  6-5 ( b)  and  6-5(c).  The  "tails"  shown  in  Figures 
6-5(d)  through  6-5(g)  are  artifacts  of  the  erroneous  branding  process 
described  in  Case  6.1. 
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Figure  6-5.  Cell  state  arp  for  ring  fire 


PUUT  81wUL*  T  ION  Tint  •  l7o,.,  »u  StCJi.ua 

flowft  .1.10...  a*..  .cuuJ«Nft _ u..tu _ 

S  On  Si f 

HOW  40  . . _  .  .. 

flow  44 

ROW  4;  . 

HQw  40 

HOW  _„.49 _ _ 

HOW  44 

flO»  43  ....  . . . . 

HO-  49  •« 


HOW  4  1 

H0«  44 

now  30 

ROW - 3ft 

ROW  37 

ROW 
ROW 

ROw 

ROW 

flaw.- 

flOw 
HOw 
ROw 
HOW 
ROW 
HOW _ 

ROW 
ROw 
ROw 
ROw 
ROw 

fow. 

ROW 
ROW 
NOW 
RQw 


29 
29 
2? 

.2ft- 
28 
.24  ... 
23 

22  .  ... 

21 

..29 _ 

19 

18 . 

1? 

1*. 


*•_ 


t  .t 


•  •ft 
-At, 

.  ft.  -  . 

ft  ft 
ft*** 
ft* 


. -  •  -  •  ft**..**.. 

*  I  ••  *  «*• 

ft*  *  -  ••  _ _ _  **-  A-  4. 

•*  *♦  ft*  ♦****#  *«« 

,»■ .  ft  ft.*.* _ ftft _ - _ JL. 


ft*  •  *•  t  ♦  •*  **  ftftftft* 

.  t  ..  -ft  ft*-.  -  _  «•  ft.  •• ft  ft*. - ftft.A - 

ft*  ♦  #•••«  *  ft  *  I 

..  .  * -  .  ..  ft*.  _  t-  -A  .  •  *.«**...  .  *  - 

•  ••  ft  ftft  ft 

— t  ft - «LJ _ « - ft  ft.  ft  *AJ*.  *_* - 

ft  ft  ftftftft 


... _ .ft _ 


_ •  _ _  *J f _ t _ 

ft  ft  ft  ft*  ft 

♦  ...ft _ _ ft* _ 


ROw 

19 

ft  ft  * 

***  • 

F.QW. 

14  . 

.  ft  *  -  ft  .  . 

•  **•*-  *  -  .  .  .  _ 

ROw 

13 

•  ft  • 

• 

*  ft  ftft** 

RQW 

12  .. 

. . ft  ftft  * 

ft 

•  ft 

•  **ft  ft . .  _ _  .....  .  . 

ROW 

11 

♦ 

* 

ft  ft 

ROW 

1* 

ft* 

ft. 

•  *• 

.  .  ft  ft* . . .  .  . . __ . 

flQw 

0 

♦  • 

«*ft 

* 

JO". 

1... 

♦ 

JL 

# 

ROW 

7 

ROw 

0 

ROw 

9 

ROW 

4 

ROW 

ROW 

3 

_  2 

ROW 

1 

COLUMNS  1234967(99 1 
NUWftIR  0»  UNTCnITED, 


234967994 1234867899 1234(6 7(99 1234867(99 
BURNING  AWO  auRNEO-OUT  CELLS  *  2194 


39s 


Figure  6-5(b).  Cell  state  map  for  ring  fire. 
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Figure  6-5(d).  Cell  state  nap  for  ring  fire. 
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Figure  6-5(f)„  Cell  state  Map  for  ring  fire 
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Figure  6-5{g).  Cell  state  map  for  ring  fire. 
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Figure  6-5(1).  Heat  production  rate  contour  «ap  for  ring  fire 
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Figure  6-5(J).  Cell  heat  production  rate  sap  for  ring  fire. 
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Figure  6-5(k).  Cell  heat  production  rate  Map  for  ring  fire. 


235 


ttmm  n hkm  (m m 

SSSaanSc  nmi 


-sarpai®r 


*  ~j^v£»ssr£^an^ 


varCT.VJSwrvftsifires" 


ounwft  NMI  KKr 


RUM  *  *a.Mtr 


m  A**ir4r.VAAA4-A3^AAAAA haaah  *  ] 
nnn  *  Aft-wA-VA.-. 4ii*^i4niuiu  k  1 

a  x  >n  -*•*•*<  >A-Atffr***tf*  it  k  n  kQ 

3*  x  xi  AAA-i-fc-to-fcAAfrAfrfrAAAfr*  K  at  at  it? 
AAAAi aJfliff-JkrA*  +  ♦  *  *  *  I* 

aa  aaaaIaa*  +  ♦•*»*»  ****&!£  KKt  3 

A'A’tf^’jSA  . . .  *  h  AAtf 5 

AA-AAAA-*  . . *  •  *  HAAfrtfKKftftg 

AAAAAh  . . ••■•  +  +  h 

-A-A'A'AA  *  •  •  •  - . ■»**♦« 

-A-  ^  »aa  . . .  if  ^rXjrAij 

•  ■  . *■< 

*  »  •  •  . . .  ■  h<. 

■  •  . 'Afr^  . 

. . «  «frfrfrfr^ 

ftftft35ft3**  . . .  «RfW  _ 

-  *  * . .  Hft.RK1SJ.tf.frC] 

ftftA3535AA.ftfr  *  ***••*«  fttfRRKtffrfrfr 

*  ft  75  71  35  71 75  75  A  H  HtHM-vftftPrKRRKRfrfrl 

H  )5  71  35  35  7575  35-frfr  A  K  K  K  K  fr! 

ft  ft  71  71  35  75  75  75  AAA  A  A' 

*  *  75  7*  7»AAAAAAA' 

91 

1 


Wff«  K  35  55  fr 

ck*** 


a  a  a  aa-V^-VA1 

A  A  ft  Aft  A-V'V1' 
AAAAAAi7A. 


AfrfrfrAfrtftfKftl 
_ '  frfrfrfrlLKttWK! 


A  AA  A  JflaflilfiU J»^ArA<A7A frfrfr  A  frfrfr  K  6f  AC  At? 
A  A  7®  A  A  Mil  iff  At  A  *  7  lift*  HAA&AKKACKKl 
AAAA'AAMft  «♦*..***  *  AAtltKAfifAfAri 

aaaaa  A  a  •* . . 

aaaaa*-* . WwjrtS 

-AA-ft-AA  **♦••■ . •  m  V§%AA'P 

A~H>A  +  H+  ** . *  4'4f'^ 

»A+*'»  *  •  •  •  . . *•»+•■< 

-T-fr-5  4  +  ♦••♦• . •  •  *  A< 

._.?***♦«•♦ . *  *  A< 

JjJpft  if  . . .  tfft 


it  tf,  frfr.fr  fr‘ 
"  fr  ■ 


»  ... 

*♦•••••  ******  Afrfrfrfr . 

*  .........  *  AKRRtf.tf.tfJ 

».«...*..»»  RKKRfrtf  - 
Aft  ft  ft  ft  75  ATS  *  *  *  *  *  •  *»  «»KTO(tlW1 
AftftftATSATSAAHH  *  *  A  A  ftUCKfSJKKRRfr 
A  ft  ft  A  7J  73  75R' A' A' A' A' A-A  ft  AftKR  K 

A  35  71  35  71  tflTftA  ATPT  P5  K  R  tf  FJ 

35  35  35  3"  35  75  A  A  Afr^T  TTyA*^*1-  *f  ft  ft  ?5  ft  ft  It 


iBmmm 


IMM  riM  Ml 


MUMtlll 

M  MM 


-jw&T’jsftarj'jEra,*-- 

—  Ml  MM  VMM  MIMM  mum  N  VUV1W  — 
I  .MM  W.l TT  IVMIVM  -  W.  Mi T 


A  A  A  A  *  M  A  **<VAr4r-trA  ■ AAA  AA  AAA  A  fr  *1 
A  A  A  A  A  A  A  AAA  A-VAA-f  fr  A  fr  frfr  frfr  tCUK 
A  A  A  A  A  A  A  A  AAAfrfrfrfrfrfrfrfrfr  frfr  fr  fr  n 
A A A  A A AAAAAAfrfrfrfrfr-frfrfrfrtffrfrfrfr 
A  A  A  A  A  A  AAA A  AAAAfrfr.frfrfr.frtf.frlK  tf  A 
A  A»  A  A  A  A  A  A  A  A  *  A  A  AfrfrfrfrfrtftKIKA if: 
AAAAAAAAAAM  A  A  *  *  *  Hfrtf  tfKtK&KA 
aaaaaaaaammv  a  a  *H-*A*tti6£ic4rtr4r6 

AAAAAAAAMM  mw  w  *  wAAHAK  && A^Afg 
A  A  A  A  Ai*  A  ammmwwwv**h  HA4cr'4s'AAP 
AAA  AAAAAMMM****ww**4r 
AAAAA  A  A  AMM3IW......AH' 

AAAAAAAA  MM.....  • 

AAAAAAAA  *»«»•*■• 

AAAAAAAa  . . .  H 

^AAAa  »  *  . . *  *  tffffffj, 

^AAA  »  . .  5-AA'?-’f4-4r 

■  ••  *  HAtfAtf-tf-tf-tfA 

.  .  .  Hfrtf.tf.tf.tf.tf.tf.fr- 

J9^J?7J.^7575^75  A  H  *  ♦  V  H  A  ft RRtf.tf.fr  tf.fr  tfj" 

7J  7}  7571777  75  RRRRfrt 

75  75  75  ft  75  73  R*' frtf5fr^  R  ft  r  r  r  g 


“js'.'aft'jsaar.rxTEftyu- 


AAAAAJrJf^rV^»-'>’WAAAAtfAtfAKK. 

A  AA  AAAAA^>^-V^vAAAA^5*A«A  A  *  ■ 
A  A  A  A  A  AAA  n  rAAAAAtf  RbLU  tf  **  - 

A  A  A  A  A  A  H  A  iliit  isW.  M  p^AAltftSC  R IA  A  t*  itfl 
AAA  73  ^.A.5t5Ht6£5t  6fK«4fe 

A  A  A  A  A  A  A  A  Jo  Jfriswv  r>  rd*5M&6£.  tf.  tf  K  #£  Iff  Ilf  4r[ 
A  A  A  A  A  A  A  Jd  fi£  CS  &  «C  K  *r  *t"C 


AA^AiiAAAijvtf-^-W'+^oiteCbllaCte'fa'trfri 

. 

-^A^A^W^AA  A  •V  +  MAKKi.,,  ,  ,  ,  ,. 
AAAAAAA'AAA  *»**»«  *f  *,4r-<Nr4'4A-4- 1 
.....  .  »..«*•  *r^r4"%^4^4"^"  3 

; ; ; ; ;  JtjUUIti 

35  fr  H  H  tf  ft  (5  R  tf.tf.tf.frfr.fr fr  ] 

5J73  A  A -F  A  ft  ft  ftgtf  tf  frfrfrfrfr  3 
^^^J979^7?^73AAfrtAt^£«Ktf  tf.frtf.frfr  3 
ft  ft  jj  3? 35  ft  73  35RRfrfrA4^M'^ft-tf  R  R K  fr  fr  frg 
ft  ft  75  ft  ft  ft  75  »  RfrJafr'M-JftJft/ftft'ft  K  K  K  R  fr  fr 7 
7!ft757575R  75  39tflfrfrfr- 
ft  ft  75  73  35  35  75  Rfrfrfr  * 
ft  35  75  35  35  75  A  fr  AAA' 

35  ft  3ft  35  7?  AAA  AAA 

J§§?g?ms 


■A-ftftKK  KRfrfr? 
ft  er  r  k  r  r  r 
555505(15  K  RF 
■tArftftprftiftfrtft 

?T?S?9T^7§£ 


Figure  6-5(1).  Hind  velocity  nap  for  ring  fire. 
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Figure  6-5(m).  Wind  velocity  naji  for  ring  fire. 
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6.6 


THREE  SMALL  FIRES  WITH  NO  AMBIENT  WIND 


As  shown  In  Figure  6-6,  three  centrally  located  5x5  candidate 
sets  of  cells  were  Ignited  with  probability  of  Ignition  of  0.5.  Figures 
6-6  through  6-6(9)  show  that  none  of  the  fires  spread  outwardly,  although 
each  fire  filled  In  some.  The  wind  velocity  maps  of  Figure  6-6(1)  and 
6-6 (m)  show  the  reason  for  this.  The  fire  induced  wind  speeds  generated 
by  these  fires  were  high  enough  to  cause  only  slight  contagion  and  insuf¬ 
ficient  to  cause  branding.  Since  these  winds  were  essentially  directed 
inward  toward  the  center  of  the  three  fires,  no  outward  spread  vas  pos¬ 
sible.  Only  Inward  spread  and  filling  in  could  have  occurred. 
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Figure  6-6.  Cell  state  map  for  three  clump  fires. 
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Figure  6-6{a).  Cell  state  map  for  three  clump  fires. 


240 


PLOT  SIMULATION  T I Mt  »  1 .nu  SECONDS 

ROMS  -  .  ..  1  TO..  SB,...  COLUMNS _ 1_IQ _ Sit 

ROW  9b 

RQw  49  . . 

ROM  4k 

ROM  47  .  . .  . . 

ROM  4k 

.  ROw _ 49 _ ; _ 

ROn  44 

RO"  43  _  .  . 

RQm  49 


ROw 

HOW 

ROw 

..ROrt _ 

41 

48 

39 

36 

. . -  ....  .  -  -  * . - 

...  »•  ft  .  ..... 

ROM 

37 

** 

•  * 

RO"  . 

3ft.  . 

- _  *  **•  .  .. 

.  -*•** _ . _ 

ROM 

3ft 

***** 

ROM  . 

34  . 

***  _  _  . 

ROW 

33 

-J1QM _ Jl. 


ROM 

31 

ROM 

3* 

RQm 

29 

ROm 

-26 

RQm 

27 

_R.QW_ 

_ 2ft. 

ROw 

29 

ROw 

24 

ROM 

21 

HOm 

22 

RQw 

21 

— BOw. 

..  29.. 

RQw 

19 

ROw  IS  ..  . . . . •••** _ _ _ 

ROM  17 

ROM  Ik  .  _ _ _  ..  .. _ *  • _  _ ... _ _ 

RO»  1ft  f  •* 

—ROM. _ .14 _ 

ROw  11 

.Rom  .  ia  .  . .  . . . . . .  . . . 

ROM  U 

ROM  ..  IB  ....  ....  .  _ . _ . .  . .  . 

ROM  9 

_.RQ"  _ ft _ _ 

row  r  . . . . 

RO"  ft 

RO"  9 

ROw  4 

RO"  1 

ROm  1 

COLUMNS  12349379 98 1234S67S94l234So789912349ft780Ol2345ft7S0B 

NUMBER  OR  i.inI ONl TED#  BURNING  AND  BijRNEO-OUT  CELLS  •  2448  92  » 


Figure  6-6(b).  Cell  state  Map  for  three  cluMp  fires. 
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Figure  6-t»(d).  Cell  state  eap  for  three  cluvp  fires. 
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Figure  6-6{e).  Cell  state  Map  for  three  clusip  fires 
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Figure  6- 6{f ) .  Cell  state  map  for  three  cluap  fires. 
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Figure  6-6(g).  Cell  state  i#an  for  three  cluep  fires 
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Figure  6-6(h).  Heat  production  rate  contour  for  three  cluap  f*res 
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Figure  6-6(1).  Heat  production  rate  contour  for  three  cluq>  fires. 
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Figure  6-6(j). 


Cell  heat  production  rate  nap  for  three  clunp  fires. 
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Figure  6-6(1).  Wind  velocity  Map  for  three  clump  fires. 
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Figure  6-6(«) ,  Mind  velocity  aap  for  three  eluwy  fires 
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SMALL  AREA  FIRE  WITH  AMBIENT  HINDS 
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As  shown  in  Figure  6-7,  a  10  x  10  cell  candidate  set  of  cells  in 
the  lower  left  corner  of  the  urban  area  was  ignited  with  probability  of 
Ignition  of  0.1.  This  location  was  chosen  because  the  effects  of  wind 
velocities  directed  from  the  lower  left  to  the  upper  right  at  speeds  of 
10,  20,  30,  40,  and  50  m/s,  were  examined  and  are  shown  In  Figures  6-7, 
6-8,  6-9,  6-10,  and  6-11,  respectively. 

These  cases  depict  the  sensitivity  of  propagation  to.  local  wind 
velocity.  There  are  also  the  first  cases  In  which  significant  collections 
of  cells  ere  seen  to  reach  burnout.  The  effects  of  this  are  clearly  shown 
in  all  of  the  maps  at  32,750  seconds,  but  are  best  seen  in  the  heat 
production  rate  contour  maps. 


253 


7LQT  situation  U*t 


1949,99 

40*9  1  TO  5<S,  COLUMNS  1  TO  50. 

40*  «Q 

40*  40  ...  ... 

40*  47 

40»  -  . -  - -  •  . - 

40*  49 

oo«  43 

»0«  •  *t  ...  -  -  _  .  - -  -  - - 

40*  41 

40*-  4»  ...  -  -  . . . -  -  -  -  .  .  - . . . . .  — 

40*  39 

40*  ..  341 - - - — ■ - 

40*  J7 

40*  36 .  — . .  - - - - 

40*  39 

40*  -  34  .  ...  -  — . . .  - .  -  •  ■ 

40*  33 

no*  ■  -38 —  - 

40*  31 

40*  34  -  -  . . 

40*  39 

4o*  a«.  ■  -  - . . . 

40*  27 

40*  29 

40*  24 . ~.  - -  - -  —  •  - . 

40*  23 

40*  22  -  — . . 

40*  21 

40*.  -  24  - 

40*  |9 

40"  19  . . . . .... 

40*  17 

40*  19 _ - 

40*  13  •• 

RUSt - 4.4 _ . _ 

40*  13 

40*  12  •  .  *  .  .  .  . .  .  - - 

40*  11  * 

40*  19  ..  .. - - _  „  . .  .  .. 

40*  9  ft  •  • 

rqr _ a _ _ _ 

40*  7 

40*  6  ..  .  .  *  . . . . — _ _ _  . 

40*  9 

40*  4  .  . . . . . 

40*  3 

40*  1 

COLUMNS  ) 2345(1 7994 1234  967999 1 234867499 1 234  567949 1234967994 
*UMBER  07  U*t6H ITEfti  bu9h ING  »*0  BORMID-OUT  CELLS  ■  2491 


SECONDS 
1  TO 


Figure  6-7.  Cell  state  nap  for  snail  area  fire  with  anbient  wind 
■  10  n/s. 


254 


*1.01  SI'HJi.*  I  Ion  T  1  "*t  a  t2!>*..,l- 

H0"3  1  TU  . .  52.  CULUflNl.  - 

»0w  5K 

(tow  49 

90*  48 


dtCOMUS 
-l-  TO - 50*. 


ROw 

47 

ROW 

46 

now 

48 

non 

44. 

ROW 

43 

ROW 

42 

NOW 

4  t 

ROW 

44 

ROM 

35 

now 

38— 

ROu 

3? 

ROW 

36 

ROW 

38 

ROW 

34 

now 

33 

now.  .  .. 

32— 

now 

31 

now 

3H 

ROW 

20 

RON 

26 

ROW 

27 

-RON..—  . 

..36 _ 

ROW 

20 

-ROW. 

24  _ 

ROW 

23 

RON 

22 

RON 

21 

RON  ..  . 

aa- 

ROW 

10 

ROw 

i« 

ROW 

\7 

RON 

16 

ROw 

18 

ROW.  .. 

IK- 

ROw 

lS 

ROW  . 

12 

ROW 

11 

ROW 

10 

now 

9 

.AON.... 

m 

ROW 

7 

ROW 

6 

ROw 

9 

ROw 

4 

ROw 

3 

ROw 

.1. 

ROw 

1 

* 

•  • 


COLUMNS  1234587*9* 1 3349678911 1234887«0*iaj4567*9«t 33*887890 


NU*1B£H  Of  UMGNlTEO#  BURNlWS  4 NO  8URn£D-0UT  CELLS  * 


2451 


Figure  6-7(a).  Cell  state  sap  for  saall  area  fire  with  awblent  wind 
*  10  m/s. 


255 


*»UUT  31*UL.t  1 1  UN  TJ“tfc  i  l/Sy.’.io  StCO'^OS 

HOWS  I  TO  3H,  -CQLUXHS- . .  t  TO....  Sa,.~ . ....... 

R0»  94 

HOW  49  .  ..  .... 

HOW  <•« 

HOW  4 7  ..  . 

HQN  40 

POw  44 

(iOw  ..49  .  . .  - _  - .  -  . 

HQW  42 


ROw 
HOW 
ROw 
HOW  - 
HO" 

ROW 
ROW 
ROW 
RO" 

ROW  . 

ROW 
ROW 
ROw 
RO" 

RO" 

ROW  .  . 

ROW 
ROW 
ROW 
ROW 
ROW 
ROw  -. 

ROW 
ROW 
ROW 
ROM 
ROw 

ROW - - 

ROw 
ROW 
RO" 

ROW 
ROw 

ROW - 4. - - - 

ROW  7  • 

ROW  6  -  -  -  *  ...  ....  . . 

ROw  9 

ROw  « 

ROt  3 

ROW  | 

COLUMNS  1234S67R941 2349679914 1 354667 9QD13349S79 90 1  234 & 

"OWBCn  07  UNICNITEO*  RURnIng  *no  0uRn£D«QOT  cells  • 


41 

414 

30 

3« - - 

37 

36 

30 
34 
33 

33  — - 

31 
3<* 

39 

38 

37 

36 — 

39 

34  . 

33 

32  . . 

31 

3« - 

19 

It  -  —. 

17 

1* 

19 

14  - _ 

13  i 

12 * 

15  * 

li  -  * 

9  «  * 


W789U 

3467 


13 
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Figure  6-7(1).  Wind  velocity  map  for  nail  area  fire  with  nblent  wind 
*  10  m/%. 
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Figure  6-7 (w). 


Wind  velocity  map  for  small  area  fire  with  ambient  wind 
»  10  m/s. 
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Figure  $-8.  Cell  state  map  for  saall  area  fire,  subient  wind  *  20  m/s. 
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Figure  6-8(a). 


Cell  state  aap  for  saall  area  fire,  aablent  wind 
-  20  ml s. 
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Figure  6-8{b). 


Cell  state  Map  for  saall  area  fire,  aMblent  wind 
*  20  m/s. 
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*5ite  ■**»  for  9M11  area  fire,  wblent  Mind 
■  20  m/%. 
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Figure  6-8(g). 


Cell  state  aap  for  area  swell  fire,  awblent  wind 
-  20  w/s. 
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Figure  6-8(h).  Heat  production  rate  contour  for  seall  area  fire, 
Mb  lent  wind  -  20  ■/$. 
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Figure  6-8(1). 


Heat  production  rate  contour  for  saall  area  fire 
artlent  wind  -  20  «/s. 
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Figure  6-8(4).  Cell  production  rate  map  for  small  area  fire, 
ambient  wind  -  20  m/s. 
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Figure  6-8(k).  Cell  production  rate  nap  for  snail  area  fire, 
mb  lent  wind  ■  20  n/s. 


.s-hlV': ; 


r. jf=c^oa^^aF^i~ccT'-^<=c<x&3^-+;^T CGcroi-t-joco  r\x=&xxcj\iPovcatt\x^a 


MMt«M  nmmu 

UftOflM  MM  M9ML 


■u  mh  m*tm  md  r*t*N  ti  hwwwi  ^  —» 
MMM«U«  M  M-MPM**  N  m  ■«*  <«« 

<—  *MT  fW  IW  <iMA  Mi  *»  MW  m  fWW  - 

(cuiwfl  tHMi  mam  raniTuoi  •  u.  m« 


*******9A**vs*i4*'VVi7Ai'iiii4i*i!',^i4i!ka 

**54**54*14-5414  94*54  5494**14  9454  5494*94540 

*94******i4**i454  ***54***9494  54**4 
********lA*****54***9!Si494141S^1V 
14 1414 14 14 14 14 14 14 14  Vi  14  VI  <4 14  Vi  14 14 IV 14 1*1414 14 14(1 
**isis*i4**5454'**i4*l4i4*i4*54****Vg 
**************54*******i4i4*b 
**<4***********VS**i4VI***is  54*4 
***9494******il**54**54**54****2 
Vl*-414ii94****l4****5A****iS54***r] 
*54-**********************  *9 
14 14 14 14 14  V4  *54  Is  V4 14 14 1*  14 1*1454 1*14 1414  V4 1414 14146 
14 14 14  54  **  14 1*  1*  *9i  Vi  14 1*  14  V 1*  Vi  1* 14 14*14  5*  14*54  V 144 
14 141454  54 14 1*1*  Vi  14 1*1*  Vi  14  *54 1*14141*  IS 14 14 74 lAli? 
1414VSS  14  *  A  *  *******54  *  54iS******h 
VilA9474*3  AA74********Vi*l4******H 
•?i  *^4  V\  54  54  54  54  Vi  14  *54  Vi  <4  <4  *54  54*14 14 14 1*  *4 14 14  54 14  ViQ 
**VS**V\<S****'4**'4**54,4*****54il 

^■Aft******* /V4*4*A***Vl****54*2 
"4  *  *i  -  V '*«*••*•■  *  *'4  -*  *  *4  '*  *  <\  -* ' 14  **  **  * *  V.  ‘4**  A  Vi  94  *h 
•:.V*\-1 '*'*•*■'  M*  *  '*41V7i  7 55  Si*******549\*g 
*,  •’*  -',*,**,  *54  *-i  ■*  *54  *  '*54*  *  IS  *9;  1***  •*  1",  A  A  54  54  54 
•“  •■',  *‘\  “i  A  -1  *  4  A  <\  54  -4  -a  is  *  94  IS  "1  •*  *  A  *  *  *4  *  5*ll 
■*■**<,•*•* '\  •*  is  *5S-41S  49s  is  is  1*-*<s  Vi  1*  ii  I*'*-3 

I  357  31  ?6 7-3  136331  3579 l^S^Sr 


w  tm  won. 

MMamii  -  m.a  msam 

■M  IN  M<  MUM  L  N  Mi 

u  *“-a,aT,as,as*„'jrsa,Wi— * 

—  — .1  MM*  M  IUMM  MM  M  flffW  * 

van  HIM  MUXITT  MMinjM  -  n.»ll 


94  7»  *  74  *)4  *  5*  54  14  5*  V .* 1*  *  *  1*  1*  *  *  5*  *  *  *  1*  *Q 
********************* ****fi 
74 * *14 *  *  ?*  *** *)4  ***)4  *14 **74 ******  *K 
***********************54*4 

**************  *54  ********  *? 
*************************0 
****************** *******8 
74  94 1*  7*  1*  V  V  5*  7*  1*  7*  V  V  1*  14 1*  V 1*  V  V  J4 1*  V 1*  *h 
*54*****  *******  ****54*****  *4 
*************************2 
74*54*  ***************  ******] 

*************************3 

*************************3 

**********54  *************  *q 
5454*****  *  ****************  *2 

*****54*******************] 

74 *54 54  A  *  <  ******************  j 
*54*14*  H  ***-,' 4***********  ****9 
*54 **********************  *Q 

*  14  *  5?i  *  *  *  ******  -A  ******  *4  *  *  *  *  2 

■*9VSS*i'4*';i<5.*-A*ViV'* ** *  ******* *0 
VHAlAVi******* liVilA*-**  *  *** IS  *  *  *0 

AA* A  a a a A a a  a  a a A A A  A*  a  *  *  *  *  *  *B 
A  il  i,  A  A  A  A  A  '  4  54 1 4  A  ’A  A  V  ’A  54  A  A  A  A  *  *  *  *il 
ii  VSIS*  54*54  54  Vi  *54 14 14  54  54  VS  *  54  <4  54  ii  VS*  *  *2 

1 36?3 1 3  5  ?  9 1 36591367913579 


*t***CH  CMOMrttM 

UWM  rm  no  DC  l 


p^r  iiununM  -  »*«  m 

MM  IN  *1.  tliMM  4  Tf  M, 

mi  4(M  mmJ*  «r.  mm*  nnn  nwifflifw, 
MM  4tN  14  HMNHIMI  l«  «M  tNN 
■M.  Ml  Cum  «TM*  «uMk*l  (««  tl  rv«ITM  • 

MVlMl*  M1M0  VW.0CITT  MCNPUi*  -  «.373D 


MNM  HUHMCM  CPNNmP 

vmm  nm  howl 

furaubsimuia  own  i 


_ _ _  j  mm  nun  »  imnm  w  mm. 

Mi  •  (M  M  MMtIM  II  TIM  IM  M 
-MU  Mr  trm  <«MM  MU  WM  «MM  l«  FltTTH  — 
MKTMJh  Ml  MI  'rO.OCITT  MCHITUW  •  W.IW 


************* 

n*r\*r\n*xxx*xx 

************* 
x************ 
************* 
*************** 
***»********* 
************* 
************* 
n  *  *  *  ******** 

*7\**X  XX****** 
JirV'i****  ****** 

nrwi*  *  *  *  ****** 
rt  x  ******* 

■^71  j?-i  -*'*-*  **•*•* 
~1  ~  71 H  ^  *h  •*  *  * '« -A 


*********** 

************Q 

************& 

************a 

*********** 

*  * *** *******& 
************& 
****** ******& 
*  * *  +  ******  +  +? 
******  ******1 

*********** *g 
************a 
*** *********p 

*  *  **********n 

*********** 
************a 
*********** *p 

*  *  *  *  *  *  <  ****  v'.h 


r:  *  *f  •- *  *  *  *  *******  *  *  *  *  *r. 

-‘\'*’?\^*-****‘r*  *******  *  ******{ 

******  *  ***********  *[ 1 

’T'.  ******  ********  *  *P 

i  35/5  1  6b>3  l  3679  1  33761  3579 


j»NNan*»*i«i»i>**»*i*aaa***A**i 

nununummnmmmmnmnnm******* 

n»»m4»»»»*»»»'»*9999999  9 

4S|  4  *  4  S  »  »»«***•»**•******♦** 

4  4  4  4  4  4  i  H»****»**A***A*14*-f 

4ras»a»»  •»*»*♦■»**»******** 

JJ  44444*i»»*»55*****»**1»* 
*^4|3**l*»  •  ♦»**1k******'»*»*1»'* 
Tin*******'***********  ***** 

V3  *  7  *  ^54- + *  **  *  ************ 

53  91 71  *54*4,  -A-4  *  A  *  *  *  *  *  *  *  *  *  -k  *  *  -s.  -tf 

*  vivsvivi**  vi-5*v>-M4  ******  -r-  *  *  *  *  *£, 

V4l474i4Vi5'S*  *  **A‘ A  *  *  *  *  *********  *1 
^ **  *54*  ****  ■*■ *************  *1 
14^4 54 is**********  *  *******  *  **P 

13533136731 35791 35791 3579 


Figure  6-8(1) .  Wind  velocity  sap  for  snail  area  fire, 
mb  lent  wind  ■  20  n/s. 
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Figure  S-Sta).  Mind  velocity  nap  for  snail  area  fire, 
onblent  wind  »  20  n/s. 
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Figure  6-9.  Cell  state  aaap  for  swall  area  fire,  aablent  wind 
*  30  m/s. 
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Figure  6-9{a).  Cell  state  Rap  for  snail  area  fire,  art  lent  wind 
*  30  r/s. 
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Figure  6~9(b) 


Cell  state  Map  for  spall  area  fire,  ambient  wind 
-  30  m/%. 
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Figure  6-9(c). 


Cell  state  Map  for  SMall  area  fire,  aMblent  Mind 
-  30  m/s 
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Figure  6-9(f). 


Cell  state  map  for  small  area  fire,  ambient  mind 
■  30  m/s. 
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Figure  6-9(g).  Cell  state  Map  for  snail  area  fire,  oablent  wind 
*  30  n/s. 
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Figure  6-9(h).  Heat  production  rate  contour  for  snail  area  fire, 
«b lent  wind  -  30  m/s. 
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Figure  6-9(j).  Cell  heat  production  rate  map  for  saall  area  fire, 
ambient  wind  -  30  m/s. 
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Figure  6-9(k).  Cell  heat  production  rate  map  for  aaall  area  fire, 
ambient  Mind  *  30  */s. 
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Figure  6-9(1).  Hind  velocity  nap  for  seall  area  fire,  anblent  wind 
a  30  «/s. 
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Figure  6-9(w).  Hind  velocity  nap  for  smI!  area  fire,  awbient  wind 
*  30  m/Sm 
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Figure  6-10(a) .  Cell  state  aap  for  seall  area  fire,  art lent  Mind 
-  40  «/s. 
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.  Cell  stete  up  for  small  area  fire,  ambient  wind 
-  40  m/s. 
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Figure  6-1 0(b) 
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Figure  6-10(c).  Cell  state  nap  for  snail  area  fire,  ambient  wind 
-  40  m/s. 
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Figure  6-10{e).  Cell  state  nap  for  snail  area  fire,  anblent  wind 
■  40  n/s. 
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Figure  6-10(f).  Cell  state  Map  for  saall  area  fire,  sapient  wind 
*  40  a/s. 
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Figure  6-l0(g).  Cell  state  map  for  swall  area  fire,  aoblent  wind 
-  40  ii/s. 
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Figure  6-10(h).  Heat  production  rate  contour  nap  for  snail  area  fire, 
ambient  wind  ■  40  n/s. 
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Figure  6-10(1).  Heat  production  rate  contour  Map  for  SMall  area  fire, 
aablent  wind  -  40  m/s. 
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Figure  6-10(j). 


Cell  heat  production  rate  «ap  for  snail  area  fire, 
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Figure  6-10{k). 


Ceil  heat  production  rate  map  for  small  area  fire, 
ambient  wind  *  40  m/s. 
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Figure  6-10(1) . 


Hind  velocity  mp  for  saall  area  fire,  Mblent  wind 
-  40  mfs. 
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Figure  6-10(m).  Hind  velocity  nap  for  snail  area  fire,  anblent  wind 
*  40  n/s. 
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Figure  6-11.  Cell  state  map  for  snail  area  fire*  anbient  wind  *  50  n/s. 
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Figure  6-ll(a).  Cell  state  map  for  small  area  fire,  ambient  wind 
-  SO  m/s. 
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Figure  6-ll(c).  Cell  state  map  for  small  area  fire,  ambient  wind 
-  50  m/s. 
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Figure  6-ll(e). 


Cell  state  aap  for  small  area  fire,  ambient  wind 
■  5  m/s. 
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Figure  6-ll(f).  Cell  state  aap  for  aall 
■  50  m/s. 
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Figure  6-ll(g).  Cell  state  nap  for  snail  area  fire, 
■  50  m/s. 
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Figure  6-ll(h) 


Heat  production  rate  contour  Map  for  snail  area  fire, 
Mb  lent  wind  -  SO  m/s. 
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Figure  6-ll(j).  Cell  heat  production  rate  wap  for  suall  area  fire, 
aeblent  wind  ■  50  m/s. 
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Figure  6-ll(k).  Cell  heat  production  rate  nap  for  snail  area  fire* 
aablent  wind  3  50  n/s. 
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Figure  6-11(1). 


Hind  velocity  «ap  for  snail  area  fire,  ambient  wind 
-  SO  n/s. 
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Figure  6-ll(n)i  Wind  velocity  aap  for  smal 1  area  fire,  ambient  wind 
-  50  w/s. 
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